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Abstract 



We review the foundations as well as a number of important applications of 
light-cone dynamics. Beginning with the peculiarities of relativistic particle dy- 
namics we discuss the choice of a time parameter as the gauge fixing within 
reparametrization invariant dynamical systems. Including Poincare invariance, 
we are naturally led to Dirac's forms of relativistic dynamics. Among these, the 
front form is our main focus as it is the basis for light-cone dynamics. We explain 
the peculiar features of the light-cone formulation such as boost and Galilei invari- 
ance or separation of relative and center-of-mass motion. Combining light-cone 
dynamics and field quantization leads to the introduction of light-cone quantum 
field theory. We show how the positivity of the kinematical longitudinal mo- 
mentum implies the triviality of the light-cone vacuum. We point out that its 
special features make the light-cone formulation a unique framework to deal with 
bound states as few-body systems based on quantum field theory. In a first ap- 
plication, we analyze spontaneous symmetry breaking for scalar field theory in 
1+1 dimensions. The importance of modes with vanishing longitudinal momen- 
tum is elucidated. For fermionic field theories, we suggest to reconstruct vacuum 
properties, like chiral condensates, from the particle spectrum. The latter can be 
obtained by solving the light-cone Schrodinger equation as we explicitly demon- 
strate for the 't Hooft and Schwinger models. Finally, we make contact with 
phenomenology by calculating the pion wave function within the Nambu and 
Jona-Lasinio model. We are thus able to predict a number of observables like 
the pion charge and core radius, the r.m.s. transverse momentum and the pion 
distribution amplitude. The latter turns out being not very different from the 
asymptotic one. 
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Chapter 1 

Introduction: Strong Interactions 



1.1 Hadrons 

The word 'hadron' originates from Greek where aSpoa means 'big' or 'strong'. It 
refers to all particles participating in strong interactions. One can therefore say 
that hadrons are the physical degrees of freedom in strong interaction physics. Ev- 
erything we know empirically in this area has been obtained by studying hadronic 
processes and properties. A (very) incomplete list of hadrons is provided by Ta- 
ble [TI]. 

There are more than a hundred additional hadrons which, together with their 
properties can be found in the reviews of the particle data group [[!]]. One dis- 
tinguishes between baryons having baryon number B = 1 and half-integer spin 
(fermions), and mesons with B = and integer spin (bosons). This great variety 
of hadrons has been classified in the sixties within the quark model of Gell-Mann 
and Zweig |2|, |[ which states that hadrons are grouped into multiplets of the group 
SU(3)f, where the 'F' stands for the quantum number 'flavor'. (Anti-)Quarks 
are postulated as the realizations of the fundamental (anti-)triplet of SU(3)p. 
This implies that quarks come in three flavors, 'up' u, 'down' d and 'strange' s 
(the anti-quarks being u, d and s), which combine to give the quantum num- 
bers of the hadrons. The quarks were postulated to be fermions carrying charges 
which are multiples of a third of the electron charge. Nowadays one knows a total 
number of Np = 6 quarks with the additional heavy flavors 'charm' c, 'bottom' 
b and 'top' t. In the quark scheme, mesons consist of a quark and an antiquark, 
(anti-)baryons of three ( ant i-) quarks. Thus, for example, one finds the following 
octet of pseudo-scalar mesons with spin-parity J n = CT, 

7T° = ^/T/2(uu-dd) 

7T~ = du . , ... ( i r 

7T + = 11(1 



K° 


= ds 


K+ 


= us 


K~ 


= su 


K° 


= sd 



7] = 1/1/3(10* + dd + ss) 
This octet structure is reminiscent in the name 'the eightfold way', the title of a 
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Table 1.1: Masses and quantum numbers of a few hadrons. The upper part 
consists of baryons (B = 1 ), the lower part of mesons (B = 0). Masses are given 
in MeV. 



particle 




mass 


isospin 


spin 


B 


proton 


V 


938 


1/2 


1/2 


1 


neutron 


n 


940 


1/2 


1/2 


1 


delta 


A 


1236 


3/2 


3/2 


1 


lambda 


A 


1116 





1/2 


1 


pion 


7T 


138 


1 








kaon 


K 


495 


1/2 








rho 


P 


770 


1 


1 






famous book by Gell-Mann and Ne'eman ||. 

It was soon realized thereafter, that the quark-flavor model is in conflict with 
the Pauli principle. The baryon resonance A ++ = uuu, with all quark spins 
pointing upwards, cannot have an anti-symmetric wave function as it should as a 
three-fermion state. The loophole is to introduce an additional quantum number 
called color, giving rise to another symmetry group, SU(3)c- Upon postulating 
that hadrons should always appear as color singlets the Pauli principle is satisfied. 

Using flavor and color symmetry, the wave function of mesons and baryons 
can thus be written as 

|M> = -^M.M^)), (1-2) 

Vive 

\B) = -^=B ijk e ahc \qUU) ■ ( L3 ) 

Here, Nc = 3 is the number of colors (indices a, b, c), My and Bijk are tensors in 
flavor space. It is obvious that the baryon state is anti-symmetric in color which 
solves the A ++ -puzzle. 

If we include spin degrees of freedom and the spatial part of the hadron wave 
function, the total hadronic state can be written symbolically as 

|hadron) = \space) (g) \spin) ® \ flavor) £g> \color) , (1.4) 

where the space, spin and flavor part is always symmetric, the color part guaran- 
teeing overall anti-symmetry. The non-spatial part, \spin) ® \ flav or) <g> \color), is 
governed by the symmetry group SU (2)5 ® SU (3) p ® SU (3)c. For the positively 
charged pion 7r + one finds, for instance, 

V Z1SI C 
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where the spatial part has been suppressed. The latter will, however, become 
relevant in later chapters. 

It should be pointed out, that the flavor symmetry is not realized in nature 
as an exact symmetry, because the states in the SU(3)p multiplets have different 
masses. The color symmetry, on the other hand, is unbroken. 

If one calculates the magnetic moments (fii = ej/2Mj) of the hadrons i = p, n, 
A by adding those of their quark constituents, one can determine the constituent 
masses of the quarks 

M u = M d = 340 MeV , M s = 510 MeV . (1.6) 

The constituent quark masses thus roughly sum up to the total bound state 
(hadron) mass. 

At this point, the question arises whether quarks can really be 'seen' exper- 
imentally. The answer is both 'yes' and 'no'. On the one hand, quarks have 
never been observed as free particles outside a hadron, a property called quark 
confinement. On the other hand, the deep-inelastic scattering (DIS) experiments 
beginning in the late sixties have revealed a rich substructure of hadrons, in par- 
ticular the nucleon (Nobel prize 1990 for Friedman, Kendall, and Taylor). This 
substructure consists of quasi-free, point-like constituents ('partons') which carry 
the quantum number of quarks. 

This large variety of facts about hadrons has to be explained by the dynamical 
theory of strong interactions which is the subject of the next section. 

1.2 Quantum Chromodynamics 

Twenty-five years after its inception || ||, [7|, Quantum Chromodynamics (QCD) 
is generally accepted as the theory of strong interactions. It describes the inter- 
action of matter (quarks and antiquarks) with intermediate gauge bosons named 
gluons. Like the other theories of the fundamental interactions, QCD is a gauge 
theory; its symmetry group, SU(3) is associated with three 'color' charges me- 
diating the strong interaction. As the gauge group SU(3) is non-Abelian, the 
(eight) gluons themselves carry color charge and therefore can interact with each 
otherQ Nevertheless, in a somewhat condensed notation, the QCD Hamiltonian 
appears very similar to the one of Quantum Electrodynamics (QED). The former 
is 

H QC B = l{E 2 + B 2 ) + ^{a-B + f3M)^ , (1.7) 
where D denotes the covariant derivative (in the fundamental representation), 

D = iV-#A (1.8) 

1 This explains the name 'gluon' which pictorially expresses the idea that these particles act 
like glue when the bind among themselves ('glueballs') or force quarks to bind into hadrons. 
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and Ai = diag(m u , m^, m s , . . .) the quark-flavor mass matrix. If its mass eigen- 
values were all equal, SU(Np) would be an exact symmetry of QCD. However, 
as already stated, the only exact (local) symmetry is its gauge invariance under 
SU(3) C . 

Note that, in contrast to the Abelian gauge theory QED, the chromo-electric 
and -magnetic fields E and B as well as the gauge potentials A are non-commuting 
(color) matrices so that there are implicit summations over color indices in Q1.7| ). 
The self-interaction of the gluon fields A becomes evident upon expanding the 
chromo-magnetic fields B in terms of the A, which, however, will not be done 
here as it can be found in any textbook on QCD. 

Like its older 'relative' QED, QCD is a renormalizable relativistic quantum 
field theory. Any infinities arising from the point-like (local) nature of the inter- 
action can therefore be consistently absorbed into a redefinition of the physical 
parameters (masses, couplings). As a result, the strong coupling parameter a s 
is not a constant but is running with the typical momentum scale of the phys- 
ical process under consideration. The microscopic reason for this are vacuum 
polarization effects: quarks screen and therefore weaken the color charge (anal- 
ogous to QED), whereas the self-interacting gluons anti-screen the color charge 
which is the dominating effect. Unlike in QED, therefore, the running coupling 
a s (Q) of QCD is weak for high momentum transfer Q (small distances). This is 
the realm of 'asymptotic freedom' where perturbative methods work. For small 
momentum transfer Q (large distances), the coupling is large, perturbation the- 
ory breaks down, and one has to utilize nonperturbative methods. Typical and 
well-established values for a s are QTJ 

a s (M z ) = a s (91.2GeV) = 0.118 ±0.006 , (1.9) 
a s (M T ) = a s (1.78GeV) = 0.33 ±0.03. (1.10) 

Thus, the nonperturbative domain is characterized by a maximum momentum 
scale of approximately 1 GeV. In some loose sense one can therefore speak of two 
relevant phases of QCD, the weak coupling phase or perturbative QCD, and the 
strong coupling phase (nonperturbative QCD). 

Let us return to our first subject, the hadrons. In principle, it is quite clear, 
what a hadron is in QCD: it is an eigenstate of the QCD Hamiltonian, 

H qcd | hadron) = M H \ hadron) , (1-H) 

where Mjj denotes the hadron massQ. The question, of course, is, whether 
the 'QCD Schrodinger equation' ( |1 . 1 1| ) can be solved. If we consider a typi- 
cal hadronic scale like, say, the nucleon radius, — 1 fm, the associated energy 



2 In writing down (1.11) we neglect any possible 'fine splitting' due to electromagnetic or 
weak interactions. 
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scale is Q ~ Hc/Rn — 200 MeV. This number tells us that we are in the low- 
energy regime which implies that the binding of quarks into hadrons is a non- 
perturbative phenomenon. In other words, a perturbative solution of the 'QCD 
Schrodinger equation' will make no sense in general. 

There are two main routes out of this dilemma. Firstly, one can try to per- 
form brute-force ab-initio calculations which involve sophisticated computer sim- 
ulations on the largest machines available. Technically, one makes use of a space 
time discretization leading to a lattice formulation of gauge theories. This pro- 
gram has been developed over more than twenty years now, and has become the 
subject of entire textbooks [|]. 

Secondly, one can rely on a reputable tradition of physics, namely model 
building. There is an abundance of hadron models on the market, the most 
popular one being the constituent quark model and variants thereof. In this 
model, one starts with a non-relativistic phenomenological Hamiltonian of the 
form 

H = H + V coni , (1.12) 

with an ad-hoc confining potential V con ^ which typically is proportional to the 
inter-quark distance r (or r 2 ). The Hamiltonian describes the dynamics of two or 
three constituent quarks with masses M q ~ 300 MeV like in ( |1.6| ) which are treated 
as parameters. The main virtue of the model consists in its rather accurate 
reproduction of the hadron masses ('spectroscopy'). An exhaustive compendium 
of results can be found in the recent reviews || [l(| [ll] . 

However, as any model, also the constituent quark model has its shortcomings. 
First of all, it works where it is not supposed to work. A hadron size of Rh — 
1 fm leads to an average constituent momentum of (p) ~ 200 MeV which is of the 
order of the constituent mass M q . Therefore, hadrons are relativistic bound states 
and a non-relativistic treatment is not appropriate. Furthermore, the relation of 
the model with QCD is rather unclear. In QCD, as in any relativistic quantum 
field theory, one expects a bound state like e.g. the pion to be of the form 

|vr) ~ ^ 2 \qq) + 4>z\qqg) + falqqqq) + ■ ■ ■ ■ (1-13) 

This means that hadrons are states containing an infinite number of quarks and 
gluons, which is consistent with the results of the DIS experiments where, with 
growing resolution Q 2 , an increasing number of partons is observed. This confirms 
that there are non- vanishing amplitudes r/>2 3 "03 ••■ to find two quarks, two quarks 
and a gluon in general to find an arbitrary number of quarks and gluons in a 
hadron. 

Another flaw of the model is that the constituent masses M q do not coincide 
with the masses in the QCD Hamiltonian ( |1 . 7| ) . This is qualitatively explained 
by the idea that the constituent quarks are effective quarks (quasi-particles) with 
effective masses originating from 'medium effects' due to the QCD vacuum. A 
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Table 1.2: The presently observed quark flavors. Q/e is the electric charge in 
units of the electron charge. The (scale dependent!) quark masses are given for 
a scale of 1 GeV. 



flavor 




Q/e 


mass 




down 


d 


-1/3 


10 


MeV 


up 


u 


+2/3 


5 


MeV 


strange 


s 


-1/3 


150 


MeV 


charm 


c 


+2/3 


1.5 


GeV 


bottom 


b 


-1/3 


5.1 


GeV 


top 


t 


+2/3 


180 


GeV 



quantitative derivation of this mass generation within QCD is to a large extent 
still missing. 

A more fundamental drawback of the constituent quark model is the following: 
it does not take into account all symmetries of QCD. This is particularly true for 
one prominent symmetry which is believed to govern most of the physics of light 
hadrons, namely chiral symmetry. This is our next topic. 



1.3 Chiral Symmetry in QCD 



If we have a look at Table |1.2| , which provides a list of all quark flavors, we realize 
that there are large differences in the quark masses as they appear in the QCD 
Hamiltonian. In particular, there is a hierarchy, 

m u , m d < rrh < m c ,m b ,m tj . (1.14) 

light quarks heavy quarks 

As the masses of heavy and light quarks are separated by the same scale as the 
perturbative and nonperturbative regime, one expects different physics associated 
with those two kinds of quarks. This expectation turns out to be true. The 
physics of heavy quarks is governed by a symmetry called 'heavy quark symmetry' 
leading to a very successful 'heavy quark effective theory' |T2[. The physics of 



light quarks, on the other hand, is governed by chiral symmetry which we are 
now going to explain. 

Let us write the QCD Hamiltonian in the following way, 

H QC v = H x + $Mi> , (1.15) 

where we recall the mass matrix for the light flavors, M. = diag(m u , m s ). To 
a good approximation, one can set Ai = 0. In this case, the QCD Hamiltonian H x 



1.3. CHIRAL SYMMETRY IN QCD 



13 



is invariant under the symmetry group SU(3)r® SU(3)l, the chiral flavor group. 
Under the action of this group, the left and right handed quarks independently 
undergo a chiral rotation. Due to Noether's theorem, there are sixteen conserved 
quantities, eight vector charges and, more important for us, eight pseudo-scalars, 
the chiral charges Q§ satisfying 

[Q a 5 ,H x }=0. (1.16) 

This states both that the chiral charges are conserved, and that H x is chirally 
invariant. Under parity, — > —Q§. Now, if \A) is an eigenstate of H x , so is 
Q§\A) with the same eigenvalue. Thus, one expects (nearly) degenerate parity 
doublets in nature, which, however, do not exist empirically. The only explana- 
tion for this phenomenon is that chiral symmetry is spontaneously broken. In 
contradistinction to the Hamiltonian, the QCD ground state (the vacuum) is not 
chirally invariant, 

Q a 5 \0)^0. (1.17) 

For this reason, there must exist a non-vanishing vacuum expectation value, the 
quark condensate, 

(W) = $riPl + Mr) • (1.18) 

This condensate is not invariant (it mixes left and right) and therefore serves 
as an order parameter of the symmetry breaking. Note that in QCD the quark 
condensate is a renormalization scale dependent quantity. A recent estimate can 



be found in |13fl , with a numerical value, 

(0|^|0)(1 GeV) ~ (-229 MeV) 3 . (1.19) 
In terms of the full quark propagator, 

j> + M(p) 

S{P) = f-M\ V ) ' ( 0) 

where we have allowed for a momentum dependent (or 'running') mass M(p), the 
quark condensate is given by 

- ~i J (gjtrSW = -4 Wc / • (1.21) 

We thus see that the involved Dirac trace yields a non-vanishing condensate 
only if the effective quark mass M(p) is non-zero. This links the existence of a 
quark condensate to the mechanism of dynamical mass generation. In this way 
we have found another argument that the bare quarks appearing in the QCD 
Hamiltonian ( |1.15|) indeed acquire constituent masses. Of course we are still 



lacking a microscopic mechanism for that. 
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Table 1.3: Masses of the pseudo-scalar octet mesons (in MeV). 



meson 


7T° 


7T± 


K°, K° 




V 


mass 


135 


140 


500 


494 


549 



Goldstone's theorem []14| now states that for any symmetry generator which 
does not leave the vacuum invariant, there must exist a massless boson with the 
quantum number of this generator. This results in the prediction that in massless 
QCD one should have an octet of massless pseudo-scalar mesons. In reality one 



finds what is listed in Table 1.3 



The non-vanishing masses of these mesons are interpreted as stemming from 
the non-vanishing quark masses in the QCD Hamiltonian which break chiral 
symmetry explicitly. They can be consistently viewed as perturbations within 
an effective field theory dubbed 'chiral perturbation theory', which has become 
one of the standard tools to make quantitative predictions for low-energy strong- 



interaction physics 1 15 



It should be stressed that chiral symmetry has nothing to say about the 
mechanism of confinement which presumably is a totally different story. This is 
also reflected within a QCD based derivation of chiral symmetry breaking, the 
16| . This model explains many facts of low-energy hadronic 

It is therefore possible that 



instanton model 



physics but is known not to yield confinement, 
confinement is not particularly relevant for the understanding of hadron structure 
PL 



1.4 Outline 

This report is divided into two main parts, devoted to the foundations and appli- 
cations of light-cone dynamics, respectively. In the first part, we begin with some 
general remarks on relativistic dynamics (Chapter 2). As a paradigm example 
we discuss the free relativistic particle which is the prototype of a reparametriza- 
tion invariant system. We show that the choice of the time parameter is not 
unique as this corresponds to a gauge fixing, the purpose of which is to get rid 
of the reparametrization redundancies. By considering the stability subgroups of 
the Poincare group, we find that there are essentially three reasonable choices of 
'time' for a relativistic system, corresponding to Dirac's instant, point and front 
form, respectively. 

The latter choice is the basis of light-cone dynamics which we analyze in detail 
in Chapter 3. We find some unique features like a maximal number of kinematical 
Poincare generators, boost invariance of relative coordinates and light-cone wave 
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functions and separation of center-of-mass and relative motion. Many of these 
features are reminiscent of what is known from non-relativistic physics. 

Chapter 4 is devoted to light-cone field quantization. We show how the 
Poincare generators are defined in this case and utilize Schwinger's quantum ac- 
tion principle to derive the canonical commutators. The relation between equal- 
time commutators, the field equations and their solutions for different initial 
and/or boundary conditions is clarified. It is found that light-cone field theo- 
ries, being of first order in the velocity, generally are constrained systems which 
require a special treatment. We use a method of phase space reduction due to 
Faddeev and Jackiw to take this into account. This leads us naturally to the 
construction of Fock space. 

In Chapter 5 we discuss what might be viewed as the most spectacular feature 
of light-cone quantum field theory, the triviality of the vacuum. We show how 
this property comes about and illustrate the difference to standard instant-form 
quantization for some very simple model field theories. 

We conclude the presentation of the foundations with Chapter 6, where we 
introduce the notion of light-cone wave functions. Their special properties make 
them useful tools for investigating hadronic physics. We show how light-cone 
wave functions of hadrons can be obtained by solving the light-cone Schrodinger 
equation. 

The first application we will study in Part II is the phenomenon of spontaneous 
symmetry breaking for scalar field theory (Chapter 7), in a setting which can 
be viewed as a simplified variant of the Higgs mechanism. We analyze how a 
non-vanishing vacuum expectation value can arise in a trivial vacuum. We find 
that the vacuum properties are carried by the field mode carrying vanishing 
longitudinal momentum. This mode satisfies a complicated constraint equation 
which we solve both perturbatively and nonperturbatively. 

As the methods of Chapter 7 do not work for fermionic fields, we have to con- 
sider alternatives which is the subject of Chapter 8. We use the particle spectrum 
of model field theories in 1+1 dimensions to reconstruct their vacuum properties. 
The spectrum is obtained by solving the light-cone Schrodinger equation. This is 
done here with unprecedented accuracy, which is made possible through the use 
of computer algebraic methods yielding analytic solutions even in high orders. 

In Chapter 8 we finally make contact with phenomenology. We calculate the 
light-cone wave function of the pion within the Nambu and Jona-Lasinio model. 
This model is known to provide a good description of spontaneous chiral sym- 
metry breaking, as it is governed by the same symmetry group as low-energy 
QCD. With the pion wave function at hand we predict a number of observables 
like the pion charge and core radius, the electromagnetic form factor and the 
r.m.s. transverse momentum. We conclude with a calculation of the pion dis- 
tribution amplitude which is important for the study of exclusive processes in 
QCD. 



Part I 
Foundations 



Chapter 2 

Relativistic Particle Dynamics 



The nature of elementary particles calls for a synthesis of relativity and quantum 
mechanics. The necessity of a quantum treatment is quite evident in view of the 
microscopic scales involved which are several orders of magnitude smaller than in 
atomic physics. These very scales, however, also require a relativistic formulation. 
As we have seen already, a typical hadronic scale of 1 fm, for instance, corresponds 
to momenta of the order of p ~ hc/1 fm ~ 200 MeV. For particles with masses 
of the order of 1 GeV and below this implies sizable velocities of 0.2 c and larger. 

It turns out that the task of combining the principles of quantum mechan- 
ics and relativity is not a straightforward one. One can neither simply extend 
ordinary quantum mechanics to include relativistic physics nor quantize relativis- 
tic mechanics using the ordinary correspondence rules. Nevertheless, Dirac and 
others have succeeded in formulating what is called "relativistic quantum me- 



chanics", which has become a subject of text books since |17|, [L8|. It should, 
however, be pointed out that this formulation, which is based on the concept 
of single-particle wave-functions and equations, is not really consistent. It does 
not correctly account for relativistic causality (retardation effects etc.) and the 
existence of anti-particles. As a result, one has to struggle with issues like the 
Klein paradox [|17], the definition of position operators [|19| and the like. 

The solution to these problems is provided by quantum field theory, with an 
inherently correct description of anti-particles that entails relativistic causality. 
In contrast to single-particle wave mechanics, quantum field theory is a (relativis- 
tic) many body formulation that necessarily involves (anti-)particle creation and 
annihilation. Physical particle states are typically a superposition of an infinite 
number of "bare" states, as any particle has a finite probability to emit or absorb 
other particles at any moment of time. With the advent of QCD, however, a 
conceptual difficulty concerning this many-particle picture has appeared and al- 
ready been discussed in the introduction. At low energy or momentum transfer, 
hadrons, the bound states of QCD, are reasonably described as bound states of 
two or three constituent quarks and thus as few-body systems. It is rather unclear 
how such a constituent picture can arise in a quantum field theory such as QCD. 



19 



20 



CHAPTER 2. RELATIVISTIC PARTICLE DYNAMICS 



The physically motivated desire to describe hadrons as bound states of a 
small, fixed number of constituents is our rationale to go back and reanalyze the 
relation between Hamiltonian quantum mechanics and relativistic quantum field 
theory. 

Quite generally, bound states are obtained by solving the Schrodinger equa- 
tion, 

^|V(t)} = #|V(t)> , (2.1) 

for stationary states, 

|V>(T)> = e-^iV(T = 0)>. (2.2) 
This leads to the bound-state equation 

W(r = 0)>=£|^(r = 0)>, (2.3) 

where E is the bound state energy. We would like to make this Hamiltonian 
formalism consistent with the requirements of relativity. It is, however, obvious 
from the outset that this procedure is not manifestly covariant as it singles out 
a time r (and an energy E, respectively). Furthermore, it is not even clear what 
the time r really is as it does not have an invariant meaning. 



2.1 The Free Relativistic Particle 

To see what is involved it is sufficient to consider the classical dynamics of a 
free relativistic particle. We want to find the associated canonical formulation 
as a basis for subsequent quantization. We will proceed by analogy with the 
treatment of classical free strings which is described in a number of textbooks 
2"U| , |2"T| . Accordingly, the relativistic point particle may be viewed as an infinitely 
short string. 

We recall that the action for a relativistic particle is essentially given by the 
arc length of its trajectory 

f 2 

S = — ras\2 = — m / ds . (2-4) 

This action[| is a Lorentz scalar as 

ds = y/g^x" (2.5) 
is the (infinitesimal) invariant distance. We can rewrite the action Q2.4|) as 

S = —m / ds y/x~^i~i* = / dsL(s), (2.6) 



1 We work in natural units, h = c = 1. 
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in order to introduce a Lagrangian L(s) and the four velocity x M = dx^/ds. The 
latter obeys 



x 2 = x^x" = 1 , (2.7) 

as the arc length provides a natural parametrization. Thus, x^ is a time-like 
vector, and we assume in addition that it points into the future, x° > 0. In 
this way we guarantee relativistic causality ensuring that a real particle passing 
through a point A will always propagate into the future light cone based at A. 

We proceed with the canonical formalism by calculating the canonical mo- 
menta as 

BL 

p» = -— - = mi' 1 . (2.8) 
dip 

These are not independent, as can be seen by calculating the square using ( |2.7|) , 

2 2-2 2 /r> n\ 

p = m x = m , (2.9) 

which, of course, is the usual mass-shell constraint. This constraint indicates that 
the Lagrangian L(s) defined in ( |2.6| ) is singular, which means that its Hessian 
W^ u with respect to the velocities, 

is degenerate. It has a zero mode given by the velocity itself, 

W^ u x u = . (2.11) 

The Lagrangian being singular implies that the velocities cannot be uniquely 
expressed in terms of the canonical momenta. This, however, is not obvious from 
as we can easily solve for the velocities, 



x^=p fl /m. (2.12) 

But if one now calculates the canonical Hamiltonian, 

H c = -p^ — L = -mx 2 + mx 2 = , (2.13) 

one finds that it is vanishing! It therefore seems that we do not have a generator 
for the time evolution of our dynamical system. In the following, we will analyze 
the reasons for this peculiar finding. 

First of all we note that the Lagrangian is homogeneous of first degree in the 
velocity, 

L(ax") = aL{x^) . (2.14) 
Thus, under a reparametrization of the world-line, 



•s 



x"{s) i-> x^(s'(s)) , (2.15) 
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where the mapping s i— > s' is one-to-one with ds' /ds > (orientation conserving), 
the Lagrangian changes according to 

L{dx^/ds) = L((dx^/ds')(ds'/ds) \ = {ds /ds)L{dx fl /ds') . (2.16) 

This is sufficient to guarantee that the action is invariant under (|2.15|) , that is, 
reparametrization invariant, 

S= dsL(dx^/ds) = ds'- — pL(dx^/ds') = S' , (2.17) 
J S1 y s / ds' ds 

if the endpoints remain unchanged, Si^ = s' 12 . On the other hand, L is homoge- 
neous of the first degree if and only if Euler's formula holds, namely 

BL 

L = —rf = - V ,x^ . (2.18) 

This is exactly the statement fl2.13|) , the vanishing of the Hamiltonian. Further- 
more, if we differentiate ( |2.18| ) with respect to x^, we recover fl2.11| ) expressing 
the singular nature of the Lagrangian. Summarizing, we have found the general 
result |^2J |23[ that if a Lagrangian is homogeneous of degree one in the velocities, 
the action is reparametrization invariant, and the Hamiltonian vanishes. In this 
case, the momenta are homogeneous of degree zero, which renders the Lagrangian 
singular. 

In order to properly incorporate the reparametrization invariance we chose 
an arbitrary parametrization, r ^ £ m (t), so that the infinitesimal arc length 
becomes 

ds 2 = (J^J dr 2 = h(r)dr 2 , (2.19) 

introducing a world-line metric (or einbein) h{r). In what follows we let the 
over-dot denote the derivative with respect to r, e.g. / = df /dr. From ( |2.19|) we 
get that 

i 2 = (t) 2 = HTh (2 ' 20) 

so that, according to ( |2.7| ), h(r) = 1 corresponds to dr = ds. Therefore, the 
world-line metric h provides a 'length' scale for the four- velocity. We write the 
action as 

S = -m J drVi 2 = -m dr ^/h[r) = J dr L(r) . (2.21) 

The canonical momenta are 

dL i M i; M m dx^ ds dx^ 

7r M = -777- = m—= = m . = . — — — = m— — = » M . (2.22) 
dx^ ^ v /7{77) V /7{(Y) ds dr ds 
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We note that they are independent of the world-line metric h, so that, upon 
squaring fl2.22j ) one again finds the mass shell constraint 7r 2 = p 2 = m 2 . However, 
if we want to express the velocities in terms of the momenta, we get 



x" 



Vhp^/m . (2.23) 



Therefore, it is exactly the world-line metric h that expresses the arbitrariness 
involved when solving for the velocities. Unlike for the momenta, the scale for 
the velocities is not fixed, unless one chooses a particular world-line metric as in 



By the same arguments as before, the canonical Hamiltonian still vanishes. 
The Lagrangian one-form thus coincides with the canonical one-form, 

L(r)dr = -p^dx" . (2.24) 

From the latter we read off the non-vanishing Poisson brackets, 

{x» , f} = -sT , (2.25) 

which will be needed in a moment. To proceed we utilize the Dirac-Bergmann 
algorithm for singular systems [[24|, |25|. ^6|. (for monographs, see | 2"2"| , |2"5 , ^ 



by adding the constraint 

= p 2 - m 2 (2.26) 

to the (vanishing) canonical Hamiltonian via a Lagrange multiplier A. This yields 
the primary Hamiltonian 

H p = X(p 2 - m 2 ) . (2.27) 
The constraint (|2.26 ) is the only one arising and has vanishing Poisson bracket 



with itself whereupon it is called first class. Technically, this implies that demand- 
ing consistency, i.e. conservation of the constraint in time r, does not lead to an 
equation determining the Lagrange multiplier. Physically, this means that the 
time development of the system is not completely determined. Such a behavior is 
typical for systems with a gauge or redundancy symmetry, an example of which 
is the reparametrization invariance discussed above. The redundancy consists 
in the fact that a single trajectory (world-line) can be described by an infinite 
number of different parametrizations. The physical objects, the trajectories, are 
therefore equivalence classes obtained by identifying all reparametrizations. An 
infinitesimal reparametrization is given by 

T \r) = T + St(t) , (2.28) 
Sx^( T ) = x^(t + 5t) - x»(t) = x^5t . (2.29) 

In systems with redundancy symmetries, these symmetries are generated by the 
first class constraints. A familiar example is provided by gauge fields (like in 
electromagnetism), where the first class constraints are represented by GauB's 
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law, the generator of (time- independent) gauge transformations. For the case at 
hand, we find using ([2.251) , 

6x» = {x M , 5G} = {xf , 6Se} = -If be = -2m^=5e . (2.30) 

yh 

Identifying 5t = (2m/yh)8e to account for the different dimensionalities, ( [2.29 ) 
and (|2.30|) coincide. Thus, the reparametrization ( |2.29| ) is indeed generated by the 
first-class constraint 9 of (|2.26|) . The Hamiltonian equations of motion become 

x" = {x M , H p } = -2\p» , (2.31) 
p" = {p»,H p } = 0. (2.32) 

It is important to note at this point that the equation for x^ contains the Lagrange 
multiplier A which is undetermined so far. This makes the whole dynamics of 
the system undetermined. If we compare ( |2.31| ) with (|2.23| ), we find that the 
arbitrariness is again encoded in the world-line metric as the multiplier is 

A = -Vh/2m . (2.33) 

To overcome the problem of determining A or h, one has to fix a gauge by choosing 
an auxiliary condition, 

x (x;t) = 0. (2.34) 
Demanding consistency, i.e. conservation of the gauge fixing condition in time, 



X = ^ + {X,H P } = ^ + \{ X ,6} = ^-2\ P ^ = 0, (2.35) 



we obtain the desired equation determining the Lagrange multiplier A, 

We see that it is crucial for the determination of A that the gauge fixing \ 
depends explicitly on the time parameter r and on at least one of the x^ . This is 
necessary for rendering the Poisson bracket {x, 0} nonvanishing. The analogue 
of this condition in the gauge theory context is well known: it corresponds to 
having a nonvanishing Faddeev-Popov matrix. For the relativistic particle, one 
can say that the gauge condition amounts to defining a time parameter r in terms 
of the coordinates x M . This suggests that x wm i n general be of the form 

X (x;t)=t-F(x) , (2.37) 

stating that r is some function F of the coordinates, r = F(x). In this case, 
( |2.36| ) becomes 

A = 1— - . (2.38) 

2(p • d)F(x) v ; 
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The quantity (p ■ d)F(x) is nothing but the projection of the four-momentum p M 
onto the normal = d^F of the hypersurface r = const. The knowledge of the 
Lagrange multiplier (or the einbein h) makes the dynamics uniquely determined, 
as (|2.31| ) leads to 

±" = -2\p» = . (2.39) 

N ■ p 

Squaring this expression or using ( 2.33|) , the world-line metric is found to be 

It thus measures the normal component of the momentum (squared) compared 
to p 2 = m 2 . The most common choice for r is the Galileian time t of, say, the 
observer rest frame, 

x = x°-T = t-T = 0. (2.41) 

This gauge choice corresponds to a particular choice for the world-line metric h, 
given by 

ds 2 = dt 2 - dx t dx i = (1 - v 2 )dt 2 = ^dt 2 = h(t)dt 2 , (2.42) 

with v % = dx l /dt. In this case, the Lagrange multiplier as derived from (|2.3Gf ) 
becomes 



From (|3TD and (ggg ) we get 

x» = -2\pi* =p»/p , (2.44) 
This yields the four-velocity squared 



2 2 1 

^ = -2 = ^ = -^ Kt) , (2-45) 



■„2 _ 

Po Po T 



in accordance with the result fl2.42|) for the world-line metric. The Hamiltonian 
generating the t-evolution is found by solving the mass-shell constraint 9 for 
p° = N ■ p, 

H = p° = vV + m 2 . (2.46) 
It generates the equations of motion 



V = { X \H} = ^, (2.47) 



V 

f = {p\H}=0, (2.48) 

where the first one is consistent with Q2.39D so that the second one coincides with 
the Euler-Lagrange equation following from ( |2.21| ). The Poisson brackets above 
are interpreted as being taken at equal time t. 
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Altogether, the gauge-fixed dynamics is well-determined as can be seen from 
the (unique) solution 

x'Ut) = P-t + x l (0) =vH + x l (0) , (2.49) 
pU 

which corresponds to linear motion with constant velocity. 

We have thus seen that it is the choice of time which makes the dynamics un- 
ambiguous. This problem and its solution are actually typical for reparametriza- 
tion invariant systems, the most important example probably being general rela- 
tivity. There, the problem of choosing a natural time variable is the main obstacle 
for developing a quantum theory of gravity (for a thorough discussion, see p0|). 

In the remainder of this section we will perform a general analysis of the pos- 
sible choices of gauge fixings, that is, of different time parameters r. The choice 
of a particular time variable corresponds to a foliation of Minkowski space into 
hypersurfaces of equal time, r = const, which in general are three-dimensional 
objects, and the time direction 'orthogonal' to them. The time development thus 
continuously evolves the hypersurface E : r = r into Ei : r = T\ > t . Put dif- 
ferently, initial conditions provided on E together with the dynamical equations 
(being differential equations in r) determine the state of the dynamical system 
on Si. 

The gauge fixing ( |2.37|) expresses the time parameter in terms of the coordi- 



nates x^. These are orthogonal coordinates which are essentially tailored to fit 
the common choice of time, r — t. On the other hand, we have explicitly seen 
within this particular example, that the gauge (or time) choice corresponds to 
an associated choice of the world-line metric. It is the latter, on which we now 
want to concentrate. Consider some arbitrary coordinates £ Q = £ a (x), which 
may be curvilinear. We imagine that the £*, i — 1,2, 3, parametrize the three- 
dimensional hypersurface E, so that the remaining one, £°, represents the time 
variable, i.e. r = £°. Then we have 

ds 2 = g, u dx»dx» = g ^ ——dCde = h^iOdCdf ■ (2-50) 

Introducing a vierbein e^(£) the metric h a p(£) is alternatively given by 

M0=a«^a(e)eV0- ( 2 - 51 ) 

This transformation is well known from general relativity, where it corresponds to 
the transformation from a local inertial frame described by the flat metric g^ v to 
a non-inertial frame with coordinate dependent metric h a p(£). For our purposes 
we write this metric in a (3+l)-notation as follows, 



(2.52) 
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Of particular interest is the component hoo, which explicitly reads 

hoo = flW d £ d to = ^ e V V o = n 2 , ( 2 - 53 ) 
where we have defined the unit vector in ^-direction 

n» = W = e"o- (2-54) 



It is related to the normal vector via 



The hij are the metric components associated with the hypersurface. The invari- 
ant distance element thus becomes (setting hoi = hi), 

ds 2 = h 00 dZ°dZ° + 2h oi dZ°dt* + Kjdgdg , 

n +2h i — + h ij — —jdT , (2.56) 

where, in the second step, we have used that £° = r. The quantities w i = d^/dr 
are the velocities expressed in the new coordinates, so that the world-line metric 
can be written as 

h( T ) = n 2 + 2h i w i + h ij w i w j . (2.57) 

For r = £° = t we find n = N = (1, 0), N 2 = 1, h { = 0, w i = v i and h {j = -dy 
leading to the world-line metric h(t) = 1 — v 2 of ( p.42 , |2.45 ). The dynamics of 
the is easily found as follows, 



? = d u C^ = d u C-^— = = , (2.5* 

^ N ■ p N ■ p N ■ p K 



where we have introduced the momenta ir l canonically conjugate to satisfying 

= (2.59) 



2.2 Dirac's Forms of Relativistic Dynamics 

At this point one needs a satisfactory criterion in order to decide which choices 
for the time variable (or the world-line metric h) are sensible. To this end it is 
not sufficient to consider only the r-development and the associated generator 
of time translations (i.e. the Hamiltonian) . Instead , one has to refer to the full 
Poincare group to be able to guarantee full relativistic invariance. The generators 
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of the Poincare group are the four momenta P M and the six operators M^ u which 
combine the angular momenta and boosts according to 

V = -e ijk M jk , (2.60) 
K* = M 0i , (2.61) 

with i, j, k = 1,2,3. These generators are elements of the Poincare algebra which 
is defined by the Poisson bracket relations, 

{P» , p^} = o , 
{M fW , P p } = g up P p - g w P v , (2.62) 
{M^ , M pa } = g^M up - g pp M va - g va M w + g vp M pa . 

It is well known that the momenta P M generate space-time translations and the 
M^ v rotations and Lorentz boosts, cf. (|2.60| , |2.61| ). In the following we will only 
consider proper and orthochronous Lorentz transformations, i.e. we exclude space 
and time reflections. 

Any Poincare invariant dynamical theory describing e.g. the interaction of 
particles should provide a particular realization of the Poincare algebra. For 
this purpose, the Poincare generators are constructed out of the fundamental 
dynamical variables like positions, momenta, spins etc. An elementary realization 
of (|2.62j ) is given as follows. Choose the space-time point x^ and its conjugate 
momentum p M as canonical variables, i.e. adopt ( |2.25| ). 

{x p , f} = -sT ■ (2.63) 

The Poincare generators are then found to be 

P p = p p , = x p f ~ x u p" (2.64) 

as is easily confirmed by checking ( |2.62| ) using ( |2.63| ). An infinitesimal Poincare 
transformation is thus generated by 

SG = -\5u^ u M pu + 5a^P p , (2.65) 

in the following way, 

Sx» = { x » , SG} = Su^Xu + 5a^ , du^ = -5uj vp . (2.66) 
The action of the Poincare group on some function F(x) is thus given by 

5F = {F, 5G} = d^FSa^ - \{x p d v - x u d p )F5u^ . (2.67) 



Though the realization ( 2.64 ) is covariant, it has several shortcomings. It does 
not describe any interaction; for several particles the generators are simply the 
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Table 2.1: All possible choices of hypersurfaces £: r = const with transitive 
action of the stability group d denotes the dimension of that is, the 
number of kinematical Poincare generators, xj_ = (x 1 ,^ 2 ). 



name 


£ 










r 


d 


instant 


x° 


= 








t 


6 


light front 


x° 


+ X 3 


= 






t + a; 3 /c = 


7 


hyperboloid 


rp 2 

Xq 


-x 2 


= a 2 


>0, 


x° > 


(t 2 - xVc 2 - aVc 2 ) 1 / 2 


6 


hyperboloid 


X 




= a 2 


> 


x° > 


(t 2 - xi/c 2 - aVc 2 ) 1 ^ 


4 


hyperboloid 


rp 2 

■' 


,->-» 2 

X 2 


= a 2 


>0, 


x° > 


{?-xl)l(*-a 2 l<?yi 2 


4 



sum of the single particle generators. This point, however, is of minor importance 
to us, and will only be touched upon at the end of the next chapter. The solution 
of the problem, as already mentioned in the introduction to this chapter, is the 
framework of local quantum field theory. More importantly, the representation 
(|2.64|) does not take into account the mass-shell constraint, p 2 = m 2 , which we 
already know to guarantee relativistic causality as it generates the dynamics. 

To remedy the situation we proceed as before by choosing a 'time' variable r, 
i.e. a foliation of space time into essentially space-like hypersurfaces E with time- 
like or light-like normals. £ should be chosen in such a way that it intersects all 
possible world-lines once and only once. Apart from this necessary consistency 
with causality this foliation appears quite arbitrary. However, given a particular 
foliation one can ask the question which of the Poincare generators will leave the 
hypersurface £ invariant. The set of all such generators defines a subgroup of 
the Poincare group called the stability group G^, of E. The associated generators 
are called kinematical, the others dynamical. The latter map £ onto another 
hypersurface £' and thus involve the development in r. One thus expects that 
the dynamical generators will depend on the Hamiltonian (and, therefore, the 
interaction) which, by definition, is a dynamical quantity. 

It is clear, however, that the stability group corresponding to a particular 
foliation will be empty if the associated hypersurface looks very irregular and thus 
does not have a high degree of symmetry. One therefore demands in addition that 
the stability group acts transitively on E: any two points on E can be connected 
by a transformation from G-£ ■ With this additional requirement there are exactly 
five inequivalent classes of hypersurfaces |HJ which are listed in Table [2.1| . 



The first three choices have been found by Dirac in his seminal 1949 paper [32 



on 'forms of relativistic dynamics'. He called the associated forms the 'instant', 
'front' and 'point' forms, respectively. These are the most important choices as 
the other two forms have a rather small stability group and thus are not very 
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useful. We have only listed them for the sake of completeness. 

It is important to note that for all forms one has lim^oo r = t, which means 
that in the non-relativistic case there is only one possible foliation leading to the 
absolute Galileian time t. This is consistent with the fact that there is no limiting 
velocity in this case implying that particle trajectories can have arbitrary slope 
(tangent vector). Therefore, the hypersurface S nr : t = const is the only one 
intersecting all possible world-lines. 

To decide which of the Poincare generators are kinematical we use the general 
formula (|2.67 ) describing their action. Imagine that S is given in the form X : 



t = F(x) as in Table |2~T . If P M or are kinematical for some /i or then, 
for these particular superscripts, we must have 

d^F = , {aPST - x u d")F = , (2.68) 

saying that the gradient and the 'angular derivative' of F have to vanish. Let us 
now discuss the different forms in more detail along these lines. 

The Instant Form 

The choice of Galileian time r = t is of course the most common one also in 
the relativistic case, and we have discussed it briefly in the preceding subsection. 
To complete this discussion, we construct the associated representation of the 
Poincare generators on E : t = 0. The idea is to explicitly saturate the constraint 
p 2 = m 2 by eliminating the variable conjugate to r = t, i.e. p° = N ■ p = 
a/p 2 + m 2 , and setting a; = in ( |2.64|) . Alternatively, one can follow Dirac [32] 
and add the constraint to any of the Poincare generators, 

P» = p v + - m 2 ) , (2.69) 
M v» = rtftf + + X ^(p 2 - m 2 ) . (2.70) 

Now we demand that the right-hand-sides should be independent of p°. Differ- 
entiating the right-hand-side of (|2.69|) with respect to p° one finds 

1 + 2X°p° = , (2.71) 

which is exactly the consistency condition (|2.35| ) for the case at hand. Thus A is 



determined as A = —l/2p° which coincides with ( |2.43j ) upon identifying A = A. 
The multipliers X 1 turn out to be zero. The same procedure applied to the M^ v 
of Q2.70 ) yields as the nonvanishing Lagrange multipliers the A * with 



A* = . (2.72) 

Altogether we obtain the following (3+l)-representation of the Poincare genera- 
tors, 

P l = P\ M ij = xV-xV, 

P° = Up, M i0 = x l u p , u p = (pV + m 2 ) 1 / 2 . 1 ' 
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This result is as expected: Compared to fl2.64| ), p° has been replaced by u p , and 
x° has been set to zero. It should, however, be pointed out that for non-Cartesian 
coordinates the construction of the Poincare generators is less straightforward, 
as will be seen in a moment. 

As a general rule we note that only the generators containing p° in the repre- 
sentation ( |2.69| , |2.70| ) lead to non-vanishing Lagrange multipliers. These are just 
the dynamical generators which are, for the instant form, given by the ordinary 
Hamiltonian P° and the boost generators M l °. In agreement with ( [2.68| ), one 
has 

d l F = (x i d j - ^d^F = , (2.74) 

so that E is both translationally and rotationally invariant confirming that the 
dimension of its stability group is six (cf. Table On the other hand, 

d°F = 1 ^ , (2.75) 
{x°d l - x i d j )F = -x l + , (2.76) 

from which we read off that, apart from the Hamiltonian, also the boosts are 
dynamical, i.e., E is not boost invariant. The latter fact is, of course, well known 
because the boosts mix space and time. Under a boost along the n-direction with 
velocity v (n = v/v), t transforms as 

t — > t' = t cosh uj + (n • x) sinh uj , (2.77) 

where u is the rapidity, defined through tanhtu = v. From ( |2.77|) it is evident 
that the hypersurface E : t = is not boost invariant. 

In obtaining the representation ( 2.73 ), we make the Poincare algebra compat- 
ible with the instant-form constraint x° = 0. An elementary calculation, using 
{x l , p 7 } = S^, indeed shows that the generators (|2.73| ) really obey the bracket 
relations ( [2.62 ). We have already seen that the Hamiltonian P° = uj p generates 
the correct dynamics. 



The Point Form 

As might have been expected, the situation is somewhat more complicated for 
the point form. The hypersurfaces of equal r = (x 2 — a 2 ) 1 / 2 are hyperboloids 
implying that curvilinear coordinates are involved. For the following discussion 
we chose a = for simplicity. Note, however, that in this case the hypersurface 
E : r = is the light-cone which is not transitive and has a smaller stability 
group of dimension four. We therefore assume r > in what follows. 
It is natural to introduce hyperbolic coordinates defined via 

x° = rcoshcj , 

x 1 = r sinh uj sin 6 cos <fi , 

x 2 = r sinh uj sin 9 sin <fi , 

x 3 = r sinh uj cos 9. (2.78) 
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The coordinates can thus be written in the compact form 

x M = rn M = r(coshc<j, e r sinha;) , (2.79) 

where e r is the radial unit vector in R 3 . The new coordinates are a straightforward 
generalization of the (three-dimensional) polar coordinates. Using the general 
formula (|2.50| ), they lead to a world-line metric 

h( r) = i- S inh* „ (f ) ' - r> S i„h* a S i„* (g ) ' . (2,0) 

The normal vector on E is iV^ 1 = n fl = x^jr which can be interpreted as a 
four- dimensional radial vector. The Lagrange multiplier is therefore given by 

A = "2iV^ = "27^' (2 ' 81) 
so that the equation of motion for the x^ becomes 

T 



■pf*. (2.82) 
p ■ x 

The Hamiltonian, that is, the generator of r-development, is 

H = N ■ p = x ■ p/r . (2.83) 

Obviously, we cannot set r equal to zero. Note that this Hamiltonian does not 
belong to the Poincare generators as it is essentially the dilatation operator gen- 
erating scale transformations [3"3fl. Its explicit r-dependence indicates that the 
dynamical system is not scale-invariant which is obvious as we consider a massive 
particle. 

Let us analyze the Poincare generators. To decide which of them are dynam- 
ical, we calculate with r = F(x) = Vx?, 

d^F = — — ^ , (2.84) 

T 

(x^d u -x u d")F = -{x»x v -x v x»)F = 0. (2.85) 
r 

Not unexpectedly, all momenta are dynamical, but Lorentz transformations are 
kinematical. This is obvious from the fact that r is a Lorentz scalar rendering S 
boost and rotation invariant. The stability group G^, thus has again dimension 
six. 

Let us finally mention that the point form, though already introduced by 
Dirac, has only rarely been used in physical applications [0, |33|, |35], |36]]. The 
reason is, of course, the curvilinear nature of S which particularly hampers the 
quantization. 

We are left with the front form which is at the basis of the present work. We 
reserve the next chapter for its investigation. 



Chapter 3 

The Front Form 



3.1 Generalities 



fa , a , a 2 , a 3 ) 



(a + , a , a 2 , a" 



For an arbitrary four-vector a we perform the following transformation to light- 
cone coordinates, 

(3.1) 
(3.2) 



where we have defined 



a T = a + a 3 , a 
We also introduce the transverse vector part of a as 

a_|_ = (a 1 , a 2 ) . 

The metric tensor becomes 



a -a 3 



( 1/2 \ 
0-10 

0-1 

\ 1/2 ) 



The entries 1/2 imply a slightly unusual scalar product 



a ■ b = g^ofW 



\a + b- 



1 „- 



1,2 . 



(3.3) 



(3.4) 



(3.5) 



According to Table [2.1| , the front form is defined by choosing the hypersurface 
X : x + = 0, which is a plane tangent to the light-cone. It can equivalently be 
viewed as the wave front of a plane light wave traveling towards the positive 
z-direction. Therefore, £ is also called a light-front. It corresponds to a gauge 
fixing 



x(x, r) = t — N ■ x 



N 



1,0,0,-1 



N 2 = 



(3.6) 



where N has been written in ordinary coordinates. We see that iV + = iV°-|-iV 3 = 
which implies that the normal N to the hypersurface lies within the hypersurface 
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Figure 3.1: The hypersurface X : x + = defining the front form. It is a 
null-plane tangential to the light-cone, x 2 = 0. 





. As N is a light-like or null vector, £ is often referred to as a null- 



plane |37|, |39|]. We have depicted the front-form hypersurface £ together with the 



light-cone in Fig. |3J. 



The unit vector in ^-direction is another null- vector, 

dx^ 



Ifl, 0,0,1) 



dx + 2 

so that n ■ N = 1 as it should. Given the scalar product 
invariant distance element 

/ dx~ dx^ dx 1 

ds 2 = g iiV dx iL dx v = dx + dx~ — dx l dx l = 

from which the metric h can be read off as 

h(x + ) = x 



dx^ 



dx + dx^ 



X^ X^ 



(3.7) 

5|), we infer the 
(3.8) 

(3.9) 



dx + dx~ l 
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Within the Hamiltonian formulation we find that (p ■ d)x = N ■ p = p + , so that, 
from Q2.38| ), the Lagrange multiplier A becomes 

A = -l/2p + , (3.10) 

implying the equation of motion x^ = p^/p + . Note that p + must be non-zero. 
The equations of motion are generated by the Hamiltonian 

H = n.p = p-/2 = ^±P-, (3.11) 



by means of the Poisson brackets 



&- = {x-,H} = ^, &* = {x l ,H} = ^. (3.12) 



and 

The solution is 



p+ = { p + , H} = , p 1 = {p i , H) = . (3.13) 



x~(x + ) = —x + + x-(0) = v~x + +x-(0) (3.14) 
p+ 

x\x + ) = —x + + x i {0)=v i x + + x i {0), (3.15) 
p+ 

and again corresponds to linear motion with constant (front-form) velocity. Note 
that the Hamiltonian (|3.11|) is not the normal projection N -p of the momentum, 
because N • p lies within E and thus corresponds to a kinematical direction. As 
for the instant form, the light-cone representation of the Poincare generators can 
be obtained by solving the constraint p 2 = m 2 for p~ , inserting the result into 
the elementary representation ( p.64j ) of the generators and setting x + = 0. The 
kinematical generators are 

P i =p\ p+ = p + , M +i = -x l p + , M+- = -x~p + , M 12 = x 1 p 2 - x 2 p x . 

(3.16) 

The dynamical ones (or Hamiltonians) are given by 

p _ = p\ + m 2 ^ m- { = x-p { -x l p~ . (3.17) 
p+ 

We thus have found seven kinematical generators, so that the front form leads 
to the largest stability group. Furthermore, it is important to note that the 
Hamiltonian P~ does not contain a square root as already pointed out by Dirac. 
However, there is a singularity at p + = which therefore is a somewhat peculiar 
point. In the gauge theory language, it corresponds to a vanishing Faddeev- Popov 
'matrix'. We will quite frequently be concerned with this issue in what follows. 
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3.2 Lorentz Boosts 

Consider the following boost in ^-direction with rapidity u written in instant-form 
coordinates, 

t — > t cosh oj + z sinh u , (3.18) 
z — ► t sinh u; + z cosh u; . (3.19) 

As stated before, such a boost mixes space and time coordinates z and t, respec- 
tively. However, if we add and subtract these equations, we obtain the action of 
the boost for the front form, 

x + -> e u x + , (3.20) 
x~ -> e^aT . (3.21) 

We thus find the important result that a boost in ^-direction does not mix light- 
cone space and time but rather rescales the coordinates! Note that x + and x~ are 
rescaled inversely with respect to each other. The scaling factor can be written 
as i 

^ , 0.22) 
1 + v 

if the rapidity u is defined in the usual manner in terms of the boost velocity v, 
tanhu; = v. One should note in particular, that one has the fixed point hyper- 
surface £ : x + = which is mapped onto itself, so that the relevant generator, 
M + ~ = 2M 30 = —2K 3 , is kinematical. However, we see explicitly that this is not 
true for x + ^ 0, where we get a rescaling. Stated differently, the transformation 
to light-cone coordinates diagonalizes the boosts in ^-direction. Therefore, the 
behavior under such boosts becomes especially simple. A pedagogical discussion 
of some elementary applications can be found in Q . 

In our discussion of the boosts we are still left with the other two kinematical 
generators, the M + \ These act on the transverse coordinates as 

5x i = {x i , M +J } 5(3 3 = x + 5fr = , (3.23) 

where the last identity holds on E. Translated to instant-form coordinates, M +l 
is a combination of boosts and angular momenta, 

M +i = K l + e ij L j . (3.24) 



The finite transformation corresponding to ( |3.23|) is found via exponentiation. In 



general, if the infinitesimal transformation is written as 

5 B A = {A,B} = t adB(A) , (3.25) 
the finite transformation is given by 

A' = exp(i adB)(A) . (3.26) 
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From ( |3.23| ) we thus obtain 

x i -»• x i + x + f3 i , (3.27) 

which coincides with the infinitesimal transformation because the exponential 
series terminates after the term linear in the boost parameter f3 l . Again, on E 
the x l remain unchanged. 

We are actually more interested in the transformation properties of the mo- 
menta, as these, being Poincare generators, are more fundamental quantities than 
the coordinates, in particular in the quantum theory ||31|| . As P M transforms as a 
four-vector we just have to replace x M by P M in the boost transformations ( p. 20 



3.21| , |3.27|) and obtain, 



P+ _^ e -p+ ; (3.28) 
P- -> e^P- , (3.29) 
P i -»■ P i + (3 { P + . (3.30) 

We remark that P + = is a fixed point under longitudinal boosts. In quantum 
field theory, it corresponds to the vacuum. Note further, that in the third identity 
describing the action of M + \ longitudinal and transverse momenta (which are 
both kinematical) get mixed. 

We can now ask the question how to boost from (P + , P*) to momenta (Q + , Q % ). 
This can be done by fixing the boost parameters u and (5 % as 

+ O 1 — P i 

- = -log^, P = %i~- (3.31) 

Obviously, this is only possible for P + ^ 0. We emphasize that in the construction 
above there is no dynamics involved. For the quantum theory, this means that 
we can build states of arbitrary light-cone momenta with very little effort. All we 
have to do is applying some kinematical boost operators. The simple behavior of 
light-cone momenta under boosts will be important for the discussion of bound 
states. The same is true for the subject of the next section. 



3.3 Galilei Subgroups and Consequences 

Let us analyze the Poincare algebra in light-cone coordinates in more detail. We 
are particularly interested in the Poisson bracket relations of the seven generators 

{M 12 ,M +i } = e lj M +j , 

{M 12 , P } = e lj P j , 

{M +i ,p-} = -2P\ (3 ' 32) 

{M +i ,P j } = -5 ij P + 
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All other brackets of the seven generators vanish. Consider now the two-dimensio- 
nal Galilei group. Its generators (for a free particle of mass m) are: two momenta 
k l , one angular momentum L = e % ^x l k\ two Galilei boosts G l = mx l , the Hamil- 
tonian H = k 1 k l /2m and the mass m, which is the Casimir generator. Upon using 
{x* } k j } = 5 ij and identifying P i «-> k\ M 12 <-> L, M +i «-> -2G*, P+ <-> 2m and 
P~ ^ H, one easily finds that ( [3.32 ) forms a subalgebra of the Poincare algebra 



which is isomorphic to the Lie algebra of the two-dimensional Galilei group. (A 
second isomorphic subalgebra is obtained via identifying M~ % <-> 2G l and ex- 
changing P + with P .) The first two identities in (|3.32p , for instance, state that 
M +l and P % transform as ordinary two-dimensional vectors. 

One thus expects that light-cone kinematics will partly show a non-relativistic 
behavior which is associated with the transverse dimensions and governed by the 
two-dimensional Galilei group. 

A very important example where these features are profitably at work is pro- 
vided by a system of many particles which, for the time being, will be assumed as 
non-interacting. Let the i th particle have mass and light-cone four momentum 

Pi = (pt,P±i,Pl) ■ ( 3 - 33 ) 

As the particles are free, the total four momentum is given by the sum of the 
individual momenta, 

p = J2^- ( 3 - 34 ) 

i 

Note that both the individual particles and the total system are on-shell, 

K = ^4^. (3-35) 

Pi 



P~ = , (3.36) 



P\^M 2 
P^ 

where M denotes the total mass. We introduce relative momentum coordinates 
X; and k i ,- via 



Pi = x,P + , (3.37) 
Pxi = XiP ± + k ±i . (3.38) 

Thus, Xi and \t±i denote the longitudinal momentum fraction and the relative 
transverse momentum of the i th particle, respectively. Comparing with (|3.34|) we 
note that these variables have to obey the constraints 

5>i = l, ^k ±l = 0. (3.39) 

i i 

A particularly important property of the relative momenta is their boost invari- 
ance. To show this we calculate, using fl3.2£j ), 

x '. = e w pt/e w P + = Xi . (3.40) 
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From this and ( |3.30| ) we find in addition 

k'_u = V'u ~ ^P'± = P±i + m - x l (P± + f3P + ) = k ±l , (3.41) 

which indeed proves the frame independence of Xi and k_y. 

Let us calculate the total light-cone energy of the system in terms of the 
relative coordinates. Making use of the constraints (|3.39|) , we obtain 

P±i + m i _ ( X « P -L + k ^i) 2 + m i 



i 

.2 



pt ~ XiP+ 




= K + E ^1 - ? om + ^ ■ (3-42) 

This is another important result: the light-cone Hamiltonian P separates into 
a center-of-mass term, 

Pcu = Pl/P + , (3-43) 
and a term containing only the relative coordinates, 

, II • (3.44) 

The second identity, which states that P~ is essentially the free invariant mass 
squared, follows upon multiplying ( |3.42| ) by P + , 

/ 2 2 

p+p- = p+p- _ p2 = m 2 = U + mi ■ (3.45) 

% 

To simplify things even more, one often goes to the 'transverse rest frame' where 
P_i_ and therefore the center-of-mass Hamiltonian Pq M from ( |3.43| ) vanish. 

Summarizing we note that the special behavior under boosts together with 
the transverse Galilei invariance leads to frame independent relative coordinates 
and a separation of the center-of-mass motion reminiscent of ordinary non-relati- 
vistic physics. This is at variance with the instant form, where the appearance of 
the notorious square root in the energy, P° = (P 2 + Mq) 1//2 , prohibits a similar 
separation of variables. 

For a system of interacting particles, we add an interaction W to the free 
invariant mass squared defined in (|3.45|) , 

P 2 = M 2 = Mq + W . (3.46) 

The total light-cone energy P~ is thus no longer given by the sum of the individual 
energies, but rather by 

P- = ^=E^H-^. (3.47) 



40 



CHAPTER 3. THE FRONT FORM 



Though we will not enter an intensive discussion of quantization at this point, 
we note that, when translated to the quantum theory, ( |3.46|) gives rise to the 



following light-cone Schrodinger equation, (P 2 — M 2 )\ip) = 0, or 

(M 2 -M 2 M) = W\ij) , (3.48) 

where M 2 now denotes the mass-squared eigenvalue. In later chapters we will 
derive and solve such equations within quantum field theory. 

Of course, in the interacting case, all dynamical Poincare generators differ 
from their free counterpart by some 'potential' like W. This has already been 
pointed out by Dirac |Q, who also stated that finding potentials which are 
consistent with the commutation relations of the Poincare algebra is the "real 
difficulty in the construction of a theory of a relativistic dynamical system" with 
a fixed number of particles. 

It has turned out, however, that Poincare invariance alone is not sufficient to 
guarantee a reasonable Hamiltonian formulation. There are no-go theorems both 



for the instant [41] and the front form [42], which state that the inclusion of any 



potential into the Poincare generators, even if consistent with the commutation 
relations, spoils relativistic covariance. The latter is a stronger requirement as it 
enforces particular transformation laws for the particle trajectories (world-lines). 
Thus, covariance imposes rather severe restrictions on the dynamical system []3I] . 

The physical reason for these problems is that potentials imply an instanta- 
neous interaction-at-a-distance which is in conflict with the existence of a limiting 
velocity and retardation effects. Relativistic causality is thus violated. This is 
equivalently obvious from the fact that a fixed number of particles is in conflict 
with the necessity of particle creation and annihilation and the appearance of 
anti-particles. 

Nevertheless, with considerable effort, it is possible to construct dynamical 
quantum systems with a fixed number of constituents which are consistent with 
the requirements of Poincare invariance and relativistic covariance |JT], WA . 



At this point one might finally ask the question whether the different forms 
of relativistic dynamics are physically equivalent. From the point of view that 
different time choices correspond to different gauge fixings it is clear that equiva- 
lence must hold. After all, we are just dealing with different coordinate systems. 
People have tried to make this equivalence more explicit by working with co- 
ordinates which smoothly interpolate between the instant and the front form 
[J5, £6, 47, |3]. In the context of relativistic quantum mechanics, it has been 



shown that the Poincare generators for different forms are unitarily equivalent 
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We are, however, more interested in what might be called a 'top-down ap- 
proach'. Our aim is to describe few-body systems not within quantum mechanics 
but quantum field theory to which we now turn. 



Chapter 4 

Light-Cone Quantization of 
Fields 

4.1 Construction of the Poincare Generators 

We want to derive the representation of the Poincare generators within field the- 
ory and their dependence on the hypersurface £ chosen to define the time evo- 
lution. To this end we follow ]33| and describe the hypersurface mathematically 
through the equation 

S : F(x) = r . (4.1) 
The surface element on £ is implicitly defined via 

da^u(x) = [ d i xN^5{F(x) - t)u{x) , (4.2) 
s J 

where, as before, = d^F(x) is the normal on E and u some arbitrary function. 
We will write this expression symbolically as 

da^ = d A x N»8(F(x) - r) . (4.3) 

The central object of this section will be the energy-momentum tensor, 

dC 

T» v = j^r-rdy - g^C , (4.4) 

with C being the Lagrangian which depends on fields that are collectively denoted 
by 0. With the help of the energy-momentum tensor we can define a generator 

A[f] = [ daMx)T^(x) , (4.5) 

where A and / can be tensorial quantities carrying dummy indices a, (3, . . . which 
we have suppressed. A[f] generates the infinitesimal transformations 

S f B(x) = f^x)d^B{x), (4.6) 
6 f x» = /"(I) , (4.7) 
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where / is now understood as being infinitesimal. In the same way as for a finite 
number of degrees of freedom, the generator A is called kinematical, if it leaves 
E invariant, that is, 

S f F = f^F = . (4.8) 

Otherwise, A is dynamical. With the energy-momentum tensor T^ v at hand we 
can easily show that kinematical generators are interaction independent. We 
decompose T Ml/ , 

into a free part Tq U = T^ u (g = 0), g denoting the coupling, and an interacting 
part (we exclude the case of derivative coupling). If A is kinematical, we have 



from fl43|, [D| fl 



A int [f] = - J d 4 x5(F - r^J^F = - J d A x5{F - r)C int 5 f F = , (4.10) 

which indeed shows that A does not depend on the interaction. Of course, we are 
particularly interested in the Poincare generators, P a and M a ^ . The momenta 
P a and Lorentz transformations M^ v correspond to the choices /" = and 
fa = x0l 9a ~ xl3 9^i respectively so that, from (|4.5| ), they are given in terms of 
T" u as 

C da^ a , (4.11) 



M 



J da^T^ - x p T^) . (4.12) 



From ( |4.8|) it is easily seen that the Poincare generators defined in (|4.11| , [4.12| ) 



act on F(x) = r exactly as described in Q2.671) . The remarks of Chapter 2 on 
the kinematical or dynamical nature of the generators in the different forms are 
therefore equally valid in quantum field theory. Let us discuss the three main 
forms of relativistic dynamics explicitly. 

The Instant Form 

We recall the hypersurface of equal time, 

£ : F(x) = iV • x = x° = t , (4.13) 
which leads to a surface element 

da» = d A x N^5(x° - r) , N» = (1, 0) . (4.14) 
Using ( 4.11| , |4.12p , the Poincare generators are obtained as 

P" = [ d 3 xT°» , (4.15) 

(4.16) 



= J d 3 x {x»T 0v - x v T°") . 



4.2. SCHWINGER'S (QUANTUM) ACTION PRINCIPLE 
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The Front Form 

Hypersurface and surface element are given by 

£ : F(x) = N ■ x = x + = r , da 11 = d A x iV M 5(x + - r) , (4.17) 

where N is the light-like four-vector N = (1, 0, 0, —1), N 2 = of ( |3.6| ). In terms 
of T Ml/ , the Poincare generators are 

P> = ij dx~d 2 x ± T +>M , (4.18) 

M"" = \ I dx~d 2 x ± {x»T +v - x u T +IM ) . (4.19) 
Jt 

The somewhat peculiar factor 1 /2 is the Jacobian which arises upon transforming 
to light-cone coordinates. 



The Point Form 

As we have already seen, the point form is slightly more involved. The hypersur- 
face of equal time, 



£ : F(x) = Vx 2 = t > , x l > , (4.20) 

has the radial vector = x^/t as its normal. In terms of the hyperbolic 
coordinates (|2.78|), the Poincare generators become 



P» = t 2 dudddcp smh^uj sin 9 XaT^, (4.21) 
Jt 

M" u = t 2 f duid6d<Psmh 2 u sinO Xa^T™ - x u T a ») . (4.22) 
Jt 

Note again that the limit r — > does not make sense. 



4.2 Schwinger's (Quantum) Action Principle 

Our next task is to actually quantize the fields on the hypersurfaces S : r = F(x) 
of equal time r. There is more than one possibility to do so, and we will explain 
a few of these. We begin with a method that is essentially due to Schwinger 



4£, 50, 51]. We define a four-momentum density 



<9(<9 M 0) 



so that the energy-momentum tensor T^ v can be written as 



T^ u = Wd v (j) - g^C 



(4.24) 
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In some sense, this can be viewed as a covariant generalization of the usual 
Legendre transformation between Hamiltonian and Lagrangian. Using the normal 
iV M of the hypersurface E, we define the canonical momentum (density) as the 
projection of IF, 

n = N ■ II . (4.25) 

Schwinger's action principle states that upon variation the action S = f d d x C 
changes at most by a surface term which (if E is not varied, i.e. dx^ = 0) is given 
by 

SG(t) = J da lx W54>= J d d x5(r - F)n5<p (4.26) 

taken between the initial and final hypersurfaces Ei and E2. The quantity SG is 
interpreted as the generator of field transformations, so that we have 

50={4>,6G}, (4.27) 



in case that E is entirely space-like (with time-like normal) |4"9| , |50| . We note in 
passing that the generator 5G is a field theoretic generalization of the canonical 
one-form used in analytical mechanics. In the last expression we have been a 
bit sloppy with the space-time arguments, so let us improve on this. As in the 
preceding chapter, we parametrize space-time with foliation coordinates with 
£° = t and the £ 4 living on E. Then we calculate 

mo = m),sG(r)} 



I d d -^'j{om),<0}m') , , (4.28 



where J(£) denotes the Jacobian associated with the change of variables x M — ► £ a . 
From the last expressions we infer that the canonical equal-time bracket must be 

{0(0 , *(o w = m) , <ow=r = J ^(o ^(t - o • (4.29) 

As usual, commutators are abstracted from the basic bracket by invoking Dirac's 
correspondence principle, that is, by replacing the bracket by i times the com- 
mutator. For arbitrary classical observables, A, B, this means explicitly, 

\A,B] = i(A^B} , (4.30) 

With minor modifications, the approach leading to (|4.29|) can also be used for 
the front form, resulting in what is called light-cone or light-front quantization. 
The canonical light-cone momentum is 

7T = N ■ n = iV • dcj) = <9 + = 2-^—6 , (4.31) 

ox~ 
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which is peculiar to the extent that it does not involve a (light-cone) time deriva- 
tive. Therefore, 7r is a dependent quantity which does not provide additional 
information, being known on E when the field is known there, so that ir is merely 
an abbreviation for d + <j). Again, the reason is that the normal iV M of the null-plane 
E lies within E. As a result, <fi and d + (f) have to be treated on the same footing 
within Schwinger's approach which leads to an additional factor 1/2 compared 
toll, 



\5<p= {0, 5G} , (4.32) 

with a front-form generator 

SG(x + ) = \ [ dx-d 2 x ± d + ^5(p . (4.33) 

The appearance of the peculiar factor 1/2 in ( |4.32j ) has been discussed at length 
by Schwinger in |5l| (see also ||52|| ). This factor cancels the light-cone Jacobian 
J = 1/2 in ( |4.33|) , so that we are left with the equal-time commutator (for d = 
4), 

[(f)(x),ir(y)] x + =y+=T = iS(x - y )5 2 (x ± - y ± ) (4.34) 

As the independent quantities are the fields themselves, we invert the derivative 
d + and obtain the more fundamental commutator 

[00), <j>(y)]x+=y+=T = -fsgn(aT - y-)5 2 (x± - y ± ) . (4.35) 



In deriving ( |4.35|) we have chosen the anti-symmetric Green function sgn(x ) 
satisfying 

^-sgn(ar) = 25{x-) , (4.36) 
so that Q4.34|) is re-obtained upon differentiating (|4.35|) with respect to y . We 



will see later that the field commutator ( |4.35[ ) can be directly obtained from 



Schwinger's method. Before that, however, we will study the relation between 
the choice of initializing hypersurfaces, the problem of field quantization and the 
solutions of the dynamical equations. 

4.3 Quantization as an Initial/Boundary- Value 
Problem 

As a prototype field theory we consider a massive scalar field in 1+1 dimensions. 
Its dynamics is encoded in the action 

S[</>] = j d 2 xC = J d 2 x (±<9 M 0<9> - |m 2 2 - V[(j)\) , (4.37) 
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where V is some interaction term like e.g. A0 4 and C = Cq + V. By varying the 
free action in the standard way we obtain 



5S= / da^W5ct) + 

'dM JM 



dC _ B dC 



(4.38) 



If we do not vary on the boundary of our integration region M, S^qm = 0, the 
surface term in 5S vanishes and we end up with the (massive) Klein-Gordon 
equation in 1+1 dimensions, 

(□+m 2 )0 = O. (4.39) 

We are going to solve this equation by specifying initial and/or boundary con- 
ditions for the scalar field on different hypersurfaces E. In particular we will 
clarify the relation between the associated initial value problems and the deter- 
mination of 'equal-time' commutators. 

It may look rather trivial to consider just the free theory, but this is not really 
true. Let us analyze what quantization of a field theory means in the light of 
the different forms of relativistic dynamics. One specifies canonical commutators 
like [(f>(x), 0(O)] x . e s ; where the equation for the hypersurface E: r = const defines 
an evolution parameter r and the notion of 'equal time' as already pointed out 
several times. E should be chosen such that the interaction does not change the 
commutator on it, i.e. the latter is a kinematical quantity. 

Now, if is a free field, the commutator, 

[0(x),0(O)]=*A(x), (4.40) 

is exactly known: it is the Pauli- Jordan or Schwinger function A which is a special 
solution of the Klein-Gordon equation ( [4.39 ). It can be obtained directly from 



the action in a covariant manner as a Peierls bracket |53 , |54|. Alternatively, one 
can find it from its Fourier representation which yields |55 

A(x) = -— [ d 2 P 5(p 2 -m 2 )sgn(p°)e- tp x 
2ir J 

= — ^sgn(x°) 6{x 2 ) Jo(mv^) 

= — i[sgn(x + ) + sgn(ar)] J (mVx+ar) , (4-41) 

We note that it is antisymmetric, A(a;) = — A(— x), and vanishes outside the 
light-cone, i.e. for x 2 < 0. 

If <f) is an interacting field one still has some information from causality. For 
space-like hypersurfaces (instant and point form) the equal-time commutator 
vanishes like in the free field case, 



[<f>(x),my<o = o 



(4.42) 
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This expresses the fact that fields which are separated by a space-like distance 
cannot communicate with each other. For the front form with hypersurface £ : 
x + = 0, the situation is slightly different. In 1+1 dimensions, £ is the light-cone 
and therefore entirely light-like. In higher dimensions, £ still contains light-like 
directions namely where x~ = xj_ = 0. This can also be seen from the (free) 
commutator ( |4.35|) which does not vanish at these points corresponding to a 
light-like separation. It may therefore well be that the interaction changes the 
light-cone behavior of the commutator. This is actually seen in deep-inelastic 
scattering, where the assumption of free-field commutation relations leads to the 
parton model and scaling [j5"6|1 , which, however, is violated in the interacting 
theory, i.e. QCD (though only weakly). 

Concluding, we can say that, if the fields are separated by a space-like dis- 
tance, the free field 'equal-time' commutators are boundary conditions also for 
the interacting theory and thus of sizable interest. 



4.3.1 The Cauchy Problem 



Before discussing the case of the front form, we would like to recall that the 
conventional quantization on a space-like surface (based on the instant form) 
corresponds to a Cauchy problem: if one specifies the field and its time deriva- 
tive on £ : x° = 0, the solution of the Klein-Gordon equation is uniquely 
determined as |57|, [37 1 



<f>(x) 



2/0=0 



d 

dy 1 (f)(y)-^A(x - y) 



(4.43) 



This solution expresses the well-known fact that the creation and annihilation 
operators in the Fock expansion of the quantum field are given in terms of the 
field and the velocity on £, 



a(p l ) 



dx e 



1 — ip x 



u p <p(x° = 0, x 1 ) + i(f)(x° = 0, x 1 ) 



(4.44) 



with Up = (p^ + m 2 ) 1 / 2 . If we compare (|4.43j ) with the two independent canonical 
commutators 



<f>(x),<j>(0) 

0(*),0(O) 



J x°=0 



J x"=0 



iA(x)\ x o =0 = , 
iA(x)\ x o =0 = -%5{x l ) , 



(4.45) 



we see that they can be viewed as Cauchy data for the Pauli- Jordan function A. 
Using the data Q4.45D in formula ( |4.43| ) we obtain a consistency condition, 







dy\<j>{y),m\^-o^-V) 



(4.46) 
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which expresses the commutator for arbitrary time differences x° in terms of its 
initial values given by ( f4.45|) . 



4.3.2 The Characteristic Initial- Value Problem 

In the following we will perform an analogous discussion for the hypersurfaces £ : 
x^ = const, which, in d + 1 + 1, constitute the light-cone. In Dirac's classification, 
the light-cone corresponds to a degenerate point form with parameter a = 0. 
One thus does not have transitivity as points on different 'legs' of the cone are 
not related by a kinematical operation. The Hamiltonians are P + and P~ , so 
that a Hamiltonian formulation appears a little bit strange and will involve two 
time parameters. We will not pursue this issue here. A Lagrangian formulation, 
however, is no problem; it just leads to the Klein-Gordon equation. Let us discuss 
how to solve it using initial data on S. 

The light-fronts x ± = are characteristics of the Klein-Gordon equation ||58 . 
Therefore, one is dealing with a characteristic initial- value problem [^, |60fl : if 
one specifies the field <fi on x + = Xq and x~ = Xq , say 4>(x + ,Xq ) = f(x + ) and 
4>(xq , x~) = g(x~), the solution of the Klein-Gordon equation is [[37], [38| 







x) = dy + f(y + )——A(x-y) + dy g(y )——A(x-y) 



y =z 



\ d 



(4.47) 

This amounts to quantization on two characteristics, x = 0, i.e., in d = 1 + 1, 
really on the light cone, x 2 = 0. The following two independent commutators, 

[<f>{x), 0(0)1^=0 = iA{x)\ x ± =0 = -|sgn(x T ) , (4.48) 

are then characteristic data for the Pauli- Jordan function. Again this can be cast 
in the form of a consistency condition, 



o 



+ 

y-=0 



+ / dy-[<j>(y),<t>{0))J^A(x-y) , (4.49) 



o 



y+=0 



which expresses the commutator in terms of its characteristic initial values. 

It turns out that, in case the field is massless, the above quantization pro- 
cedure is the only consistent one (in d=l+l), if one wants to use light-like hy- 
persurfaces [pT|]. 

4.3.3 The Hamiltonian Formulation 



In this subsection we want to analyze whether a (Hamiltonian) quantization on 
a single light-front is possible in the massive case. Note that in the instant form 
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with its corresponding Cauchy problem everything is quite clear. The data for 
and translate into data for and the momentum n, which is independent of 0. 
We already know that within the front form there is no independent momentum 
as the normal of a null-plane lies within the null-plane. If one does not want to 
use a second light-front in order to avoid an awkward Hamiltonian formulation 
with two times, it seems that one only needs half of the data as compared to the 
Cauchy or the characteristic initial- value problem ]52|. This seems strange and, 
indeed, is not the case as we will see. 



771 

x + ,x') = — — / dy~G(x~,y~)(j)(x + ,y~) , (4.51) 



Let us write the Klein-Gordon equation ( |4.39| ) in terms of light-cone coordi- 
nates 

(d + d~ +m 2 )0 = O , (4.50) 
and solve for = d<p/dx + (which amounts to inverting d + = 2d/dx~). This gives 

i~ 
4 

where the Green function G is an inverse of the derivative d/dx~ , 

l-G{x-,y-) = 5{x--y-). (4.52) 

The Hamiltonian equation of motion is given by the Poisson bracket with H = 
\m 2 j dx~(p 2 , 

777? f 

0(x+,a; _ ) = — / dy' {(j)(x + ,x~) , 4>(x + ,y~)} 4>(x + ,y~) (4.53) 

with a presently unknown bracket of the fields 0. Clearly, Euler-Lagrange and 
Hamiltonian equation of motion, (|4.51|) and (|4.53|) have to be equivalent. This is 
realized if one identifies 

(0(x + ,x-) , 0(x+,y-)} = -lG(x-,y~) , (4.54) 



Upon quantization, this should, of course, coincide with ( |4.35| ) specialized to 



d = 1 + 1. To make the whole procedure of solving the equations of motion 
unambiguous, we need an initial condition for on S and a further condition 
which makes the inversion of d + unique. Thus, as announced, we need again two 
conditions. 

After quantization the basic bracket ( |4.54|) becomes the free field commutator 
according to the correspondence principle. The free field commutator on the light 
front x + = thus essentially is the Green function G. In momentum space, G 
is given by (a regularization of) the expression l/p + . In this sense, quantization 
on x + = amounts to inversion of the light-cone momentum p + (viewed as a 
distribution). It is well known that there is no unique inverse (or regularization) 
of l/p + [JB2| . The possible regularizations differ by terms proportional to S(p + ) 
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which correspond to homogeneous solutions of ( |4.52| ). To solve this problem, 
one has to specify additional conditions with respect to the variable x~ which is 
conjugate to p + . 

It turns out that data on a second characteristic, say (p(x + , x~ = Xq) = f(x + ), 
do not work: they do not lead to an antisymmetric Green function in ( |4.51| ) 
which is needed for the commutator. Therefore one does not obtain equivalence 
between Euler-Lagrange and Hamiltonian equations of motion, which is a severe 
inconsistency. As a result, there is no Hamiltonian formulation with data on two 
characteristics and a single light-cone Hamiltonian. 

There exists, however, an alternative. We follow the philosophy of an approach 
called 'discretized light-cone quantization' [|63], |64], |65| and enclose the system in 
a finite volume, — L < x~ < L. Upon demanding periodic boundary conditions 
for the field, 

4>(x + ,x~ = -L) = 4>(x + ,x~ = L) , (4.55) 

which are to hold for all light-cone times x + , the longitudinal momenta become 
discrete, 

kl = 2nn/L , n E Z . (4.56) 

As the spatial boundaries are not separated by a space-like but a light-like dis- 
tance, one has to check whether the boundary conditions ( ^4.55| ) are consistent 
with the equation of motion and thus with relativistic causality. This will be done 
in the next subsection. Assuming this consistency, the Euler-Lagrange equation 
( |4.51| ) becomes 



TV? f 

-— / c/y~G p (x~,?/~)0(x + ,|/" 



(4.57) 



where G p is uniquely determined as the periodic sign function, 



G p (x ,y~ 



x 



-sgn(x 



y 



y 



2L 



(4.58) 



so that there is no arbitrariness left. From ( f4.54| ) we obtain a unique and consis- 
tent commutator, 



[x),<f)(0)] x+=0 



--G p (x~ 



x 



2Sgn(x 



2L 



(4.59) 



Apart from the second term, which is a finite volume correction, this coincides 
with ( |4.35| ) . We are left with checking the consistency of the procedure above. 



4.3.4 The Initial-Boundary- Value Problem 

The attempt to solve the Klein-Gordon equation with periodic boundary condi- 
tions in x~ poses an initial-boundary-value problem which is explicitly defined 
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through the conditions 



\x = 0, x ) 
(x + , x~ = L) 



9{x ) , 
d)(x + , x' 



-L) 



(4.60) 
(4.61) 



where ( 4.6C ) represents the initial data and ([4.61 ) the boundary conditions. The 
unique solution of the Klein-Gordon equation obeying these conditions is obtained 
via Fourier transformation, 



(x,x 



(4.62) 



The on-shell energy is k n = m 2 /k+ = m 2 L/27rn, and the g n are the Fourier modes 
of the function g(x~) specified in (|4.60|) , i.e. 



9n 



1 

2L 



dx~ g{x-)e in ™~ /L 



(4.63) 



Investigating the Klein-Gordon equation in momentum space one finds that the 
zero mode go is vanishing (see also Section 4.4). It is therefore already omitted 
in ( 4.62Q . From this solution it is obvious that data on a second characteristic 
are not necessary. In the quantum theory, ( |4.63|) corresponds to the peculiar fact 
that the light-cone Fock operators are expressible in terms of the field alone; the 
velocity is not needed |66|, |52j . 

Another consistency check is to demand periodic boundary conditions for the 
solution Q4.47D of the characteristic initial-value problem. Together with the data 
f(x + ) this would over-determine the problem, which can only be avoided, if / is 
expressible in terms of g. This is indeed the case, 



-r pc+ioo 

Alt t J C — lOO 



dpe 



H L [g,p] 



exp(—m 2 L/2p) 



(4.64) 



where Hi\g,p\ is a complicated but known functional of the data g [55 



In addition, let us check the consistency of the solution (|4.53|) for an important 
example, the periodic Pauli- Jordan function, 



Ar(x + ,X ) = 



E 



_ — ik n x + 12— in-KX /L 

47m 



(4.65) 



In analogy with ( 4.49|) this leads to the light-cone commutators 

-1 

x+=0 



<f>(x),<P(0) 
0(x),0(O) 



"2 L2 Sgn(x_) - k 



= gix 



J_v- (-1 

47T ^ 



n 



exp 



m L , 

i- x 

Aixn 



(4.66) 
f(x+), (4.67) 
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where ( f4.66|) coincides with ( f4.59| ), as expected. Due to the periodic boundary 



conditions the second commutator given by f(x + ) is not independent but satisfies 
exactly (|4.64 ) with g given by ( }4.66|) . As in the case of the Cauchy problem 



and the characteristic initial-value problem one can write down a consistency 
condition, 



[0(x),0(O)] = / dy-[<f>(y),cf>{0)} ^A L (x - y) 



(4.68) 

/+=0 



which expresses the commutator through its value on one characteristic (y + = 0) 
and the periodic Pauli- Jordan function. 

4.3.5 Fermions 

For fermions, light-cone quantization is slightly more involved because not all of 
the fermion field components are independent. This is easily seen by considering 
the free massive Dirac equation (again, for simplicity, in 1+1 dimensions), 

(i$ - m)ifj = . (4.69) 

If we write the two-spinor ip in terms of a lower and upper component as 

*=(£)■ (*■*» 

the Dirac equation (|4.69|) splits up into 

d'ipi = —imip2 , (4-71) 
d+^ 2 = -imipr . (4.72) 



The important point to note is that (|4.72 ) represents a constraint as d + is a 



spatial rather than a temporal derivative. It should be solved for the dependent 
lower component ip2 m terms of the independent upper component xpi, 

^ 2 = -zm(«9 + )-Vi • (4.73) 

Obviously, at this point, the ubiquitous inverse of d + arises again, and we have to 
make sure that it becomes well defined. Like for the scalar field, the mentioned 
ambiguity affects the quantization procedure. The covariant anti-commutator of 
two fermion fields with arbitrary space-time separation x is given in terms of an 
invariant function S, 

{iP(x),$(0)} = iS(x) , (4.74) 
where S can be inferred from the Pauli- Jordan function as 



S(x) = [i(p + m)A(x) 



(4.75) 
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From this representation it is obvious that S obeys the Dirac equation (as it 
must), as A is a solution of the Klein-Gordon equation. We thus expect that 
S will also contain dependent degrees of freedom. To make these explicit, we 
introduce a two-by-two matrix notation and write ( |4.74|) in the form 

[iPi(x) , Vj(0)} = iSijix) , i,j = 1,2 . (4.76) 

According to ( |4.75| ), the components iSij are 

iSn(x) = -d + A(x) , (4.77) 
iSi 2 (x) = iS 2 %(x) = imA(x) , (4.78) 
iSn(x) = -d~A{x) . (4.79) 

Comparing the second and the third identity we learn that S22 is a dependent 
quantity satisfying the constraint (|4.72|) , 

d + S 22 = -imS 12 . (4.80) 



Of course, taking into account ([4.76 ), this is consistent with the fact that ijj 2 is 
the dependent field component. There is, however, one subtlety that deserves 
a separate analysis. If one performs the derivatives in ( |4.75| ) and writes the 
Pauli- Jordan function (|4.41| ) as 



A(x) = D(x) J (mVx 2 ) , (4.81) 
the iS^ are explicitly found to be 

iS n (x) = 6(x~) + mD(x)^x+/x~J 1 (mVx 2 ) , (4.82) 

^12 (#) = imD(x)Jo(mVx^) , (4.83) 



^22(2;) = 5(x + ) +mD(x)^x-/x+J 1 (mVx 2 ) . (4.84) 

There are several remarks in order. First of all, the first identity provides the 
canonical anti-commutator of the independent field components, 

Mx) , Vl(0)) =iS n (x + = 0,x-) =6(x-) . (4.85) 

J x+=0 

In addition, we note that Sn and S22 are related by an exchange of x + and x~ . 
Going back to instant-form coordinates, this amounts to the spatial inversion 
x 1 — > —x 1 , that is, a parity transformation. Obviously, as parity exchanges light- 
cone space and time, it is a dynamical and therefore complicated transformation. 
Furthermore, the leading light-cone singularities of the Sa are mass independent, 
as is well known from, say, the space-time analysis of deep-inelastic scattering 
\ $% . In particular, we have 

lim iS 22 (x) = iS% 2 = iS 22 (x~ = 0) = 5(x + ) (4.86) 

m— >0 



54 



CHAPTER 4. LIGHT-CONE QUANTIZATION OF FIELDS 



On the other hand, if we formally solve ( |4 . 8 0| ) for S22, we obtain 

S22(x) = -im(d + )~ 1 Si2 , 



(4.87) 



which, at first glance, seems to be in conflict with ( f4.86| ) upon performing the limit 
m — > 0. One might be tempted to add the homogeneous solution S® 2 to the right- 
hand-side of ( [4.871 ). This would imply that S22 is not entirely determined in terms 
of Su, and one would have to provide S22(x + ,x~ = 0) as a boundary condition 
for all light- cone times x + . Analogous arguments apply to the relation between 
ip2 and ipi. In the language of the preceding subsections this would correspond 
to a characteristic initial-value problem. Fortunately, these difficulties can be 
circumvented, if one makes use of parity, as we are now going to show. 

We assume knowledge of the field tpi at light-cone time x + = 0. For arbitrary 
times x + it is then given by 



ipi(x) 



b{k + )e~ i% - x + d\k + ) 



Sk-x 



where 



rv Jb — t^th Ju \ cyf\i JL j rv — lit j rv ; 

and the Fock operators obey the anti-commutation relations 

{b{k + ) , tf{q + )} = {d{k + ) , d\q + )} = 5{k + - q 
Plugging (|4.88| ) into ( [4.73| ), the dependent field becomes 

dk 



2,7T 



b(k + )e- lk - x + d\k + )i 



ik-x 



from which one infers the anti-commutator 



iS 2 2 (x) = — 



rrv 



dk' 



(k 



+^2 



r-ik-x + e ik-xl 



(4- 



(4.89) 



(4.90) 



(4.91) 



(4.92) 



Note the m? in front of the integral which is consistent with ([4.87|) as S12 is of 
order m. The trick is now to invoke parity by substituting k + = m 2 /k~ with k~ 
as the new integration variable. This absorbs the m? in front of the integral and 
results in 

dk~ 

exp(—im 2 x~/k~ —ik~x + ) . (4.93) 



2tt 



iS 2 2{x) = 

Evaluating this for x~ = 0, we finally obtain 

iS 22 {x + , x' = 0) = {^2(2:) , ^2(0) } x - =0 



5(x' 



(4.94) 



We reemphasize that this result was obtained by only assuming initial data (in 
terms of the usual Fock expansion) for the independent field component ipi, and 
parity. Data on the second characteristic are once more unnecessary. 
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Altogether, we have thus shown that, for massive bosons and fermions in 
space-time dimension d = 1+1, quantization on a single light front x + = is 
possible. However, there are special requirements to be met. In the bosonic case, 
periodic boundary conditions in the spatial direction lead to a consistent and sim- 
ple light-cone Hamiltonian formulation. Data on a second characteristic are not 
necessary. The result justifies the approach of discretized light-cone quantization 
f)3l , [64], |65|. For fermions, we have used parity as an additional input. 

For d = 3+1 analogous arguments apply because the solution of the Klein- 
Gordon equation is still of the form (|4.47|) with just additional integrations over 



the transverse dimensions [37, 38 



4.4 First-Order Systems 

In this section we will encounter yet another method of Hamiltonian reduction, 
i.e. of identifying and quantizing the physical degrees of freedom. It will serve 
us to construct the light-cone Fock space. We will restrict our attention to the 
bosonic case. According to ( f4.37|) , the Lagrangian of a scalar field in 1+1 di- 



mensions, when written in terms of light-cone coordinates, takes the following 
form, 

L[(j>] = \ J dx- {\d + (j)d-(t) - \m 2 ct) 2 - V[(j))) . (4.95) 

The important thing to note is the fact that the Lagrangian is linear in the 
light-cone velocity d~<p = 2d(p/dx + . This means that the appropriate method 



for quantization is the one of Faddeev and Jackiw for first order systems []68 



called "(Constrained) Quantization Without Tears" by Jackiw |69|]. It avoids 



many of the technicalities of the Dirac-Bergmann formalism and is in general 
more economic. It reduces phase-space right from the beginning as there are 
no 'primary constraints' introduced. The method is essentially equivalent to 
Schwinger's action principle, especially in the form presented by Schwinger in 



The approach becomes particularly transparent for the case at hand if we 
again discretize the system in momentum space 0, by enclosing it in 
a finite spatial volume, — L < x < L. We proceed by Fourier expanding the 
classical free field as in ([4.62 ), 



<P(x + ,x-) = a + J2g n e~ i{ktx ' /2 + knX+/2) ■ (4.96) 

Note that we have allowed for a zero-momentum mode a in the expansion 
( |4.96| ). In what follows we will show that it actually vanishes. If we define 
a n = (47r|ri|)" 1//2 , separate positive and negative modes in (|4.96j) , and use the 



fact that a„ n = a* (as <fi is real), the summation can be restricted to run over 
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only positive modes. At light-cone time x + = the field thus becomes 

0(x+ = 0, x~) = a + Y (a n e- in ^~ /L + a* n e mnx ~ /L ) , (4.97) 

Plugging this into the free part of the Lagrangian we get (discarding a total time 
derivative) 

nx 2 L 

L [a n , a ] = —i}, a n a* n - m 2 La\ - - — a* n a n = -i/^ a n a* n - H , (4.98) 

n>0 n>0 n>0 

with if denoting the Hamiltonian. The Euler-Lagrange equations are 

- ia n + a n = , (4.99) 

Ann 

2m 2 La = , (4.100) 

where we have defined a n = da n /dx + . The first equation, Q4.99j ), is just the free 
Klein-Gordon equation which can be seen upon multiplying by The second 
identity, ( |4.10G| ), is non-dynamical and thus a constraint which states that a free 
field does not have a zero mode, a = 0. 

The method of Faddeev and Jackiw is based on equivalence of the Euler- 
Lagrange and Hamiltonian equations of motion. The latter are 

yri J~j 

a n = {a n , H} = 22 ^-^ {a n , a* k } a k , (4.101) 

fc>0 

Obviously, this coincides with Q4.99Q , if the canonical bracket is 

K,<} = -i4n- (4.102) 
The constraint ( |4.100| ) is obtained by differentiating the Hamiltonian, 

OH 

— = 2m 2 La Q = . (4.103) 
oa 

Let us briefly show that the approach presented above is equivalent to Schwinger's 



5I|. From ( |4.98[ ) we read off a generator 

5G = -i^ajal (4.104) 

n>0 

effecting the transformation 

Sa* n = {< , 5G} = -iJ2 i< » o*} Sa* k > (4-105) 

fc>0 
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which in turn implies the canonical bracket ( [4.102| ). 

Quantization is performed as usual by employing the correspondence princi- 
ple, so that, from ( 4.102| ), the elementary commutator is given by 



The Fock space expansion for the (free) scalar field <fi thus becomes 
<f){x + = 0, or) = — ( a n t 

* — _ \ Ann, V 



-irnrx I L i f inizx /L 



(4.106) 



(4.107) 



n>0 



Note that the Fock 'measure' 1 / yAwn does not involve any scale like the mass m 
or the volume L. This is at variance with the analogous expansion in the instant 
form which reads 



{x,t = 0) = -= Y] 



(ct„ 



+ oLe 



t r> ik n x 



) 



(4.108) 



where — L < x < L, k n = nn/L, and [a n ,ay = 5 mn . Obviously, does depend 
on m and L. We will discuss some consequences of this difference in the next 
chapter. 

We can use the result (|4.107|) to calculate the free field commutator at equal 
light-cone time x + , 



•«<>)], 



E 



-imrx /L 



Ann 



-sgn(aT 



X 

2L 



(4.109) 



which coincides with (|4.35|) , ( [4.59 ) and ( 4.66|) . The commutator ( 4.109| ) has orig- 
inally been obtained in using the Dirac-Bergmann algorithm for constrained 
systems. The Faddeev-Jackiw method, however, is much more economic and 
transparent. In particular, it makes clear that the basic canonical variables of 
a light-cone field theory are the Fock operators or their classical counterparts. 



The a n with, say, — N < n < N in ( 4.96 ) can be viewed as defining a (2N + 1)- 
dimensional phase space. A phase space, however, should have even dimension. 
This is accomplished by choosing a polarization in terms of positions and mo- 
menta, here a n and a^, with n > 0, and by the vanishing of the zero mode, a = 0. 
We will later see that this vanishing is a peculiarity of the free theory. 



Chapter 5 

The Light-Cone Vacuum 



5.1 Basics 

One of the basic axioms of quantum field theory states that the spectrum of the 
four-momentum operator is contained within the closure of the forward light-cone 
(in momentum space) 0, [71]. The four- momentum P M of any physical, that is, 
observable particle thus obeys 

P 2 > , P° > , (5.1) 

which is, of course, consistent with the mass-shell constraint, p 2 = m 2 . In other 
words, there are no tachyons having imaginary mass (m 2 < 0). The tip of the 
cone, the point P 2 = P° = 0, corresponds to the vacuum. From the spectrum 
condition ( |5.1|) we infer that 

P 2 - P 2 > Pi > or P° > |P 3 | . (5.2) 

This implies for the longitudinal light-cone momentum, 

P+ = p° + p 3 > |P 3 | + p 3 > . (5.3) 

We thus have the important kinematical constraint that physical states must have 
non-negative longitudinal momentum, 

(phys|P + |phys) > . (5.4) 

The spectrum of P + is thus bounded from below. Due to Lorentz invariance, the 
vacuum |0) must have vanishing four- momentum, and in particular 

P + |0) = 0. (5.5) 

Therefore, the vacuum is an eigenstate of P + with the lowest possible eigenvalue, 
namely zero. We will be interested in the phenomenon of spontaneous symmetry 
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breaking, i.e. in the question whether — roughly speaking — the vacuum is degen- 
erate. We are thus going to analyze whether there is another state, \p + = 0), 
having the same eigenvalue, p + = 0, as the vacuum. If so, it must be possible to 
create this state from the vacuum with some operator U, 

\p+ = 0) = U\0) , (5.6) 

where U must not produce any longitudinal momentum. Note that within or- 
dinary quantization such a construction is straightforward and quite common, 
for example, in BCS theory. A state with vanishing three-momentum can be 
obtained via 

|p = 0) = J d 3 kf(k)a\k)a j (~k)\0) , (5.7) 

where / is an arbitrary wave function. Evidently, the contributions from modes 
with positive and negative momenta cancel each other. It is obvious as well, 
that within light-cone quantization things must be different as there cannot be 
an analogous cancelation for the longitudinal momenta which are always non- 
negative. Instead, one could imagine something like 

dk + J d 2 k ± f(k ± )5(k + )a j (k + , k ± )a\k + , — kj_)|0> . (5.8) 

The problem thus boils down to the question whether there are Fock operators 
carrying light-cone momentum k + = 0. As we have seen in the preceding section 
there are no such operators, and a construction like ( |5.8| ) is impossible. 

The only remaining possibility is that, if U contains a creation operator 
a^(k + > 0) carrying longitudinal momentum k + ^ 0, there must be annihila- 
tors that annihilate exactly the same amount k + of momentum. Thus, after 
Wick ordering, U must have the general form 

C/=(0|[/|0) + J dk + f 2 {k + )a){k + )a{k + ) 

k+>0 

+ J dp + J dk + f 3 (k + ,p + )a\p + + k + )a(p + )a(k + ) 

p+>0 k+>0 

+ J dp + J dk + h{k + ,p + )a){p + )a){k + )a{k + + p + ) 

p+>0 k+>0 

+ ... . (5.9) 

It follows that the light-cone vacuum |0) is an eigenstate of U, 

U\0) = (0\U\0)\0) . (5.10) 

As we only deal with rays in Hilbert space, the action of U on the vacuum does 
not create a state distinct from the vacuum. This result is actually rather general. 
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Any quantity that is obtained by integrating some functional of the fields over 
longitudinal space, i.e., 

F[<j>\ = [ dx-F[(j>] , (5.11) 



is of the form ( |5.9| ), because the integration can be viewed as a projection onto the 
longitudinal momentum sector k + = 0. The most important examples for such 
quantities are the Poincare generators, as is obvious from the representations 
(|4.18| , |4.19|) . This implies in particular that the trivial light-cone vacuum is an 



eigenstate of the fully interacting light-cone Hamiltonian P , 

P"|0) = 2#|0) = (0|P-|0)|0) . (5.12) 

This can be seen alternatively by considering 

P + P~\0) = p-p + \0} = , (5.13) 

which says that P~\0) is a state with k + = 0, so that P~ must have a Fock 
representation like U in (|5.9|) . 

The actual value of (0|P"|0) is not important at this point as it only defines 
the zero of light-cone energy. Note that within the instant form the Fock or trivial 
vacuum is not an eigenstate of the full Hamiltonian as the latter usually contains 
terms with only creation operators where positive and negative three-momenta 
compensate to zero as in Q5.7| ). The instant-form vacuum thus is unstable under 
time evolution. Such a vacuum, a typical example of which is provided by ( |5.7| ), 
is called 'nontrivial'. 

In the next section we perform an explicit comparison of nontrivial instant- 
form vacua and their trivial light-cone counterparts. 



5.2 Trivial versus Nontrivial Vacua 

Let us discuss some simple field theoretic models which are exactly solvable so 
that one can explicitly determine the ground state and the particle spectrum 
built on top of it. 



5.2.1 Constant Source 

Consider first a massive scalar field in d = 1 + 1 coupled to a constant external 
source p. We solve this model field theory both in the instant and front form. 
Using coordinates interpolating between the two forms the same model has been 
studied by Hornbostel . We work in a finite spatial volume to circumvent any 



possible infrared problems. The Lagrangian density is 

£ = \d^<P - \m 2 ^ + p0 = C + p<p . (5.14) 
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Note that the external source explicitly breaks the reflection symmetry — > — 
of the free Lagrangian C. Upon completing the square, the Lagrangian ( 5.14Q 
becomes 

£ = \d^<p - ±m 2 (0 - p/m 2 ) 2 + p 2 /2m 2 . (5.15) 

The minimum of the (classical) potential is thus located at = p/m 2 = 0o, which, 
in the quantum theory, corresponds to a non-vanishing vacuum expectation value, 

(n\<p\n) = O • (5.16) 

The state denotes a 'non-trivial' vacuum in the following sense. If we expand 



the field in Fock space like in (|4. 108 



(x, t = 0) = -L= -1= + A ie~ lknX ) , ( 5 - 17 ) 



where — L < x < L, k n = irn/L, u 2 = k 2 + m 2 and 

[^,4,]=^, (5.18) 

we must have 

A n |fi) ^ (5.19) 

for at least one n. The vacuum state \Q) must therefore contain particles cor- 
responding to the A-modes of 0. If we denote the 'trivial' vacuum, i.e. the one 
annihilated by the A n , as usual by |0), the two vacua must be related as in (|5.7f) 
via 

|Q) = £/L4 n ,4]|0>- (5.20) 

Let us determine the explicit form of U. To this end we note that the system 
under consideration is a field theoretic version of the shifted harmonic oscillator. 
If we define a translated field 

^ = 0-0 o , (5.21) 
this does not have an expectation value in as 

= (n\(f>\n) -0o = . (5.22) 

Thus, from the Fock expansion of <p 

<p{x, t = 0) = ^=J2 ~W= + a ^~ lknX ) . ( 5 - 23 ) 



2L ^ 

we infer that the a n annihilate the non-trivial vacuum \Q) and thus correspond to 
quasi-particle operators built on the ground state Put differently, ip can be 
viewed as a free field with vacuum \Q), as is also obvious from the Klein-Gordon 
equation, 

(□ + m 2 )0 = (□ + m 2 )ip + p = p , (5.24) 
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which implies that tp obeys a free equation of motion. The relation between the 
A- and a-modes is easily found. As and if are related by a translation, we must 
have 

<p = W<pU , (5.25) 
where U is the translation operator 

U = exp(-ipU/m 2 ) = exp(-i(j) U) . (5.26) 

Accordingly, the operator U from ( |5.20| ) is determined in terms of the momentum 
II conjugate to 0, which is the integral over the momentum density 

tt(z, t = 0) = (A n e tKx - Aie~^ x ) . (5.27) 



The momentum II is thus essentially the zero mode of the density 7r, 

L 

IL(t = 0) = J dx7r(x, t = 0) = iVml (4 - A ) . 



(5.28) 



Putting everything together we obtain the explicit representation of U in terms 
of Fock operators, 

U[A n , 4] = expfv^L 0o (4 - A))] • (5-29) 

Note that this quantity is entirely determined by the zero momentum modes A\ 
and Aq. Plugging ( |5.29|) into ( |5.20|) , one finds that \ fl) is a coherent state of zero 
momentum modes, 

|0) = exp ( — mL0g/2) exp (VmL<j> A^ |0) , (5.30) 

where we have used the Baker-Campbell-Hausdorff formula to factorize U. Be- 
ing a coherent state, \Q) must be an eigenstate of the annihilators A n . This is 
confirmed by a simple calculation yielding 



A n ^ \ty = , A \Q) = V^Lc^o \n) . (5.31) 

The non-vanishing vacuum expectation value thus arises from the action of the 
zero mode operators, 

L 

(n\<t>\n) = ^Jdx (Qim = -±={si\a + Al\n) = <p . (5.32) 

—L 

The quasi-particle operators are obtained by the unitary transformation ( |5.25 ), 

a n = UA n U^ = A n , (n^0) 

a = UA U^ = A Q - V^L(po , (5.33) 
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As might have been expected, only the zero modes get shifted through the action 
of U, namely by the c-number —yJmL 0o- This transformation is a particular 
example of a canonical (or Bogolubov) transformation. It conserves the commu- 
tation relations ( |5.18|) 



i.e. 



5„ 



(5.34) 



and diagonalizes the Hamiltonian, as is shown by the following calculation, 



H 



- dx 
2 J 

-L 



7T 2 + 



dx A 



+ m' 



dx i 



Vw„44 - —VmZ(A + A 
^— ' m 



22 u n an a n - (2L) 



2m 2 



Obviously, the Hamiltonian becomes diagonal in the a-modes. 
find a vacuum energy 

E vac = (n\H\n) = -(2L) : '^ 



(5.35) 



In addition, we 



2m? 



(5.36) 



which is consistent with the shift in energy density in (|5.15|) . We thus have solved 
the problem entirely: the Hamiltonian is diagonal and its ground state as well as 
the excited states are determined. The entire effect of adding a constant source 
is a rearrangement of the vacuum state accompanied by a lowering of the ground 
state energy. The spectrum of excited states which is determined by the mass m 
(i.e. the curvature of the oscillator potential) remains unchanged. 

Let us see now whether we can do the same within light-cone quantization. 
Again, we enclose the system in a finite spatial box, — L < x < L and demand 
periodic boundary conditions in x~. We perform the same decomposition of the 
fields, 

= if + 0o , (5.37) 

where 0o represents the vacuum expectation value, and look for the translation 
operator U relating and (p. The momentum (density) conjugate to is 

7T = d + (p . (5.38) 

We know from the discussion of the initial- value problem that tc does not provide 
any new information: knowledge of implies knowledge of it as d + is a spatial 
derivative. Its zero mode is therefore determined through the spatial boundary 
conditions according to 

L L 

d 

dx~ —-<f>(x-) = 2[(f>(x- =L)- 4>(x~ = —L)\ 
ox~ 

-L -L 



n 



dx 



IT (X 







(5.39) 
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and we note that it is vanishing. This has far-reaching consequences. First of all 
we note that U is the identity, 

U = exp(-z0 o n) = 1 , (5.40) 

so that the vacuum does not change under action of U, 

= |0> , (5.41) 

nor do the Fock operators. How can this observation be reconciled with the 
decomposition (|5.37|) , i.e. the existence of a non- vanishing vacuum expectation 
value 0o ? To this end we consider the Klein-Gordon equation for 0, 

(□ + m 2 )0 = (d + d- + m 2 )<p = p , (5.42) 

where the external source p provides an inhomogeneous term. Upon integrating 
this over x~ the derivative term vanishes as in (|5.39[) , and we are left with 

L 

1 



— / dx (f)(x ) = p/m =0o- (5.43) 



which is nothing but a straightforward generalization of (|4.100| ). The zero mode 
of thus coincides with the vacuum expectation value 0o- It is a c- number, 
entirely determined by the source (or 'interaction') p. This is at variance with 
ordinary quantization where we have 



L 

1 J dx ^ x ) = ^L=(A + 4) ^ 0o . (5.44) 



Here, the zero mode is operator valued and yields the vacuum expectation value 
only after acting on the non-trivial vacuum |n) as is evident from 



The decomposition ( 5.37|) , when written in Fock space, becomes 

<& x ~) = E ( a ne lkix ~ /2 + aie- iktx -' 2 ) + O . (5.45) 

n>0 

Obviously, the Fock operators in and if are the same, and the zero mode of 
ip, which is a free field, vanishes in accordance with ([4.100 ). As a result we 
conclude that, in order to determine the vacuum expectation value of 0, we have 
to calculate its zero mode which is constrained by the Euler-Lagrange equation 
of motion. This will be discussed exhaustively in Chapter 7. 
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5.2.2 Mass Shift: Bosons 



Let us consider two free scalar fields, and <pi with different masses, m and 
mi [0- We stay in d = 1 + 1 but go over to a continuum formulation. The 
Heisenberg fields obey the Klein-Gordon equations 



(□ + mf) <fn = 0, 

and have the instant-form Fock expansions 

dk 1 



0,1, 



00, t) 



(H(k) e~ iki - x + a\{k) e iki - x 



(5.46) 



(5.47) 



2tt v ^0) 

with Ui(k) = (k 2 + mf) 1//2 , ki ■ x = &i(k)t — kx, and the commutation relations 

[<H{k),4(p)] = 2n6(k-p) . (5.48) 

It is important to note that the instant-form fields even at t — (i.e. the 
Schrodinger fields) depend on the masses rrii both explicitly via u>i and implicitly 
via the Fock operators Oj and a\. Accordingly, also the (Fock) vacua defined 
through 

a i (k)\n i ) = Q, (5.49) 

are mass dependent. This can be seen explicitly in the functional Schrodinger 
picture |73| in which the states are functionals *&[(p] = (<p\^) of the classical field 
ip, and the momentum operator tt(x) acts as the functional derivative —i5/5ip(x). 
In this very intuitive formulation of quantum field theory, the free field vacua 
( J5.49 ) are represented by Gaussian functionals [ 75 ] , 



(<p\fti) = det 1 ^ 4 (u;j/7r) exp 
= det 1 ^ 4 (u;j/7r) exp 



dxdy <p(x) iOi{x - y) <p(y) 



dk 
2^ 



<p(-k) uJi(k) <f(k) 



(5.50) 



where the mass dependence resides in the 'covariances' cjj. If we rewrite the 
Klein-Gordon equation for <pi as 



(□ + ml)(j)i = -(ml - mg)0i = 5m 2 0i 



(5.51) 



we can equivalently view (f>i as an interacting field with the interaction given by 
the Hamiltonian 



H mt = i / dxbm 2x2 



(5.52) 



The field <pi in the interaction picture, and the Heisenberg field <f)Q thus satisfy 
the same free field equation of motion. This implies that the associated fields in 



66 



CHAPTER 5. THE LIGHT-CONE VACUUM 



the Schrodinger picture are the same. If we choose t = as the time where the 
different dynamical pictures coincide, we have the common initial conditions, 



(po(x, t = 0) = 0i (x, t = 0) 



(5.53) 



and analogously for the momenta 7Tj. This relates the Fock operators according 
to 

ai(k) = a k a (k) + (3 k a\(-k) , (5.54) 
with the coefficients a k and (3 k given by 



a k = \ uji(k) / uj (k) + \/u! (k)/ui(k)j , 
Ac = \ (\A>i(fc)M)(fc) - \/ uo(k) I ui(k) 



(5.55) 
(5.56) 



As in the preceding subsection, we can interpret the a\{k) as the annihilators 
of quasi- or dressed particles which by some (unspecified) mechanism acquire an 
effective mass m\ ^ m . Using the fact that a\ — /3| = 1, one can define an angle 
6 k such that 

cosh^ = a k , sinh^fc = (3 k . (5.57) 
Therefore, the transformation ( 5.54Q can be rewritten as 



Q i(fc) — a>o(k) cosh^ + a (—k) sinh^ 



(5.5J 



which is the standard Bogolubov transformation for bosons. As this transforma- 
tion is unitary, it can equivalently be written with the help of a unitary operator 
U, 

ai (k) = Ua (k)W , 



given by 



U = U{d k ) = exp 



dk 

47T 



a (k)a (-k) - a\(k)al(-k) 



(5.59) 



(5.60) 



The vacuum |f^i) is thus found to be 



dk 



= U\Q, ) = exp / — 9 k a (k)a {-k) - aJ(fc)aJ(-A;) |fi ) 



(5.61) 



In quantum optics, such a state is called a 'squeezed' state J74|. It corresponds to 
a harmonic oscillator where the parabola representing the quasi-particle potential 
is 'squeezed' compared to the parabola of the particles (a -modes). It is exactly 
this 'squeezing' which leads to a different curvature at the minimum thus implying 
a shift in the mass, m —>■ m. From the point of view of the particles, the vacuum 
is non-trivial, as ao|^i) ^ 0. 
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Again, within light-cone quantization, the situation is drastically different. 
We have seen that the initial field at x + = does not involve any scale, in 
particular no mass scale. As a consequence, both the light-cone Fock operators 
and the light-cone vacuum are mass independent [|66"1 , |75| , |76| . Heuristically, one 
can argue that the front form is equivalent to a formulation in the 'infinite- 
momentum frame' [77j which is obtained from the instant form through the limit 

In this frame, all masses can be 
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of a Lorentz boost with velocity v — > c 
neglected compared to the longitudinal momentum, and thus should only play a 
minor role |76| . 



Technically, for the case at hand, the initial condition analogous to ( p.53| ), 



4>i(x ,x + = 0) = <j> Q {x~ 



0) 



(5.62) 



implies that the a - and ai-modes coincide so that U becomes the identity. Once 
again we thus find that the light-cone vacuum is not changed by the interaction, 
but remains trivial. 



5.2.3 Mass Shift: Fermions 



Consider finally the same situation for fermions, however, now in d — 3 + 1 ]79] . 
Our starting point are the two Dirac equations 

(i# - mi) V>i = , * = 0, 1 . (5.63) 
We expand the Heisenberg fields as 
f d 3 k 1 r 

Mx) = E J (27r)3/2 ^(k, A)ui(k, A)e-*- + d\(k, X)v t (k, A)e*-» 

(5.64) 

with LUi an d fct as m ( p-47p . The elementary spinors u and v are those of Itzykson 
and Zuber |80| , multiplied by a factor (2m i ) 1 / 2 to allow for a straightforward 
massless limit. The Fock operators obey the canonical anti-commutation rela- 
tions, 

{&i(k,A),&I(p,A)} = {dj(k,A),4(p,A)} =5 xx ,5 3 (k-p) . (5.65) 

As before, the Fock operators with subscripts i = and i = 1 (corresponding to 
bare and dressed particles, respectively) are connected by a Bogolubov transfor- 
mation, 

&i(k, A) = Ub (k, A)f/ f , di(k, A) = Ud (k, A)f/ f , (5.66) 

which conserves the canonical anti-commutation relations. Like in the bosonic 
case, U is found from the common initial condition, 



V>oOM = 0) = if> 1 (x,t = 0) 



(5.67) 
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A somewhat lengthy calculation involving a number of spinor matrix elements 
finally yields |79[ 

MM) = a fc & (M)+&4(-k,A) , (5.68) 
di(k,A) = a k do(k,\)-p k bl(-k,\), (5.69) 



with coefficients 



«fc = sJW +7fc) 
ft = 



1(1 



7fc 



7fc) 
m mi + A; 2 



(5.70) 
(5-71) 
(5.72) 



Obviously, we have a 2 , + /3 
to write 

&i(k,A) = 
di(k,A) = 



1, so that this time there is an angle 9 k allowing 



&o(k, A) cos^fc + d\{— k, A) sin^ 
c?o(k, A) cos 6^ — &o(— k, A) sin 6^ 



The unitary operator U is found to be 



(5.73) 
(5.74) 



(5.75) 



defining the nontrivial vacuum \fli) = U\Q ), with &o|^i)> d \Qi) 0- 

Upon initializing the fields on the null-plane x + = 0, however, things again 
become very different. The fermion fields ipi are decomposed into ipi = 
where only the are independent fields (see Section 4.3.5 and App. A). These 
obey a Klein-Gordon equation, 



U = ex P \-J2 d 3 k9 k [&o(-k, A)do(k, A) + &J(k, A)4(-k, A) 



= (id + )-\-di + m 2 ) , i = 0, 1 



(5.76) 



which is obtained after eliminating the The (independent) Heisenberg fields 
have the Fock expansion 

d 2 k 



°° dk + 



167T 3 



bi(k, X)u it+ (k, X)e- ik ^ + d\(k, X)v i& (k, X)e^ x 

(5.77) 



where k = (fc + ,kj_) and 

\c ■ . nr* — — lc~ nr* Jr I — lc~^~ T 1 . "V 1 — — \c I Is* . nr* 



(5.78) 



with the fcr being on shell, fcr = (kj_ + m 2 )/k + . Again, it is crucial to note that 
both the Fock operators and the spinors and Vi t+ are mass independent (see 
App. A). As a result, the common initial condition, 



ipo,+ (x,x + = 0) = ipi, + (x,x + = 0) 



(5.79) 
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leads to a trivial Bogolubov transformation with U = 1 [ 79|, ^TJ , like in the bosonic 
case . The light-cone vacuum, therefore, is stable under mass generation and the 
same for both fields i/>i,+, namely the trivial Fock vacuum. This is corroborated 
by taking the infinite- momentum limit of the coefficient 7^, 



lim 7 fc 

|k|-+oc 



(5.80) 



which also makes the Bogolubov transformation ( |5.68| , |5.69| ) trivial. In Chapters 8 
and 9 we will discuss how one can obtain a fermion condensate in the trivial light- 
cone vacuum (see also [fTD]). 



Chapter 6 

Light-Cone Wave Functions 



In this chapter we collect some basic facts about the eigenvalue problem of the 
light-cone Hamiltonian, or, in other words, about the light-cone Schrodinger equa- 
tion and its solutions, the light-cone wave functions. We follow the discussion 
of Brodsky and Lepage [3JJ and also use their conventions. Accordingly, this 



chapter contains only few original results, but we think it worthwhile to include 
to make our presentation more self-contained. 



6.1 Definition 

Let us first stick to a discrete notation and, for the time being, stay in 1 + 1 
dimensions. We thus have a Fock basis of states 

|n) = 4|0) , 
\m,n) = 4,4)0) , 

\m, ...,nj>r) = 4 1 ---4j ) ■ (6- 1 ) 
This leads to a completeness relation defining the unit operator in Fock space, 

1 = [Q)(Q| + ^ \n)(n\ + \ ^ \m, n) (m, n\ + . . . 

n>0 m,n>0 
oo 

= + E \ni,...,n N )(n u ...,n N \. (6.2) 



N=l m,...,njv>0 



An arbitrary state can thus be expanded as 

Wl = ^H^ 71 ) + I (m,n\ip)\m,n) + ... . (6.3) 



n>0 m,n>0 
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The sums are such that the longitudinal momenta in each Fock sector add up to 
the total longitudinal momentum of Note that the vacuum state does not 
contribute as it is orthogonal to any particle state, (O^) = 0. The normalization 
of this state is obtained as 



<v#> 



£|«>I 2 + ! E 



n>0 
oo 



m,n>0 



Ea?7 E IK 

N=l rti,...,n JV >0 



(m, n 



+ . . . 



(6.4) 



Let us assume that the state \ip) corresponds to a bound-state obeying the light- 
cone Schrodinger equation (|3.48| ), 



(M 2 -M*)\ij) = W\ 



(6.5) 



We want to project this equation onto the different Fock sectors. For this we need 
the eigenvalues of the free invariant mass squared when applied to an iV-particle 

state 

\N) = \n u ...n N ) . (6.6) 



We find that Mq is diagonal as expected, 



N 



i=l 



(N\M 2 \N') = (N\P + E^l^') = EtT^I^ = M n(N\N') 



N 



m 

.1,; 



(6.7) 



The light-cone Schrodinger equation thus becomes a system of coupled eigenvalue 
equations, 



(M 2 -M?) 
(M 2 -M|) (kl\ip) 



(l\W\m) (l\W\mn) ... 
(kl\W\m) (kl\W\mn) ... 



{km 



(6.8) 

Clearly, this represents an infinite number of equations which in general will prove 
impossible to solve unless the matrix is very sparse and / or the matrix elements are 
small. The former condition is usually fulfilled as the interaction W in the light- 
cone Hamiltonian typically changes particle number at most by twoQ. Assuming 
the matrix elements to be small amounts to dealing with a perturbative situation. 
This will be true for non-relativistic bound states of heavy constituents, but not 
for light hadrons which we are mainly interested in. We shall, however, encounter 
situations where the magnitude of the amplitudes decreases enormously with the 



^^Note that terms with only creation operators are forbidden by fc + -conservation. Still, in 
1+1 dimensions, things can become messy as interactions with polynomials of arbitrary powers 
in 6 are allowed 
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particle number N, so that it is a good approximation to restrict to the lowest 
Fock sectors. In instant form field theory this has long been known as the Tamm- 



Dancoff method [§4|, |§ . 

Let us turn to the more realistic case of 3+1 dimensions in a continuum 
formulation. In the discussion of the Klein-Gordon equation we have seen that 
the solutions are of the form 

" ' ' (6.9) 

where ft denotes the on-shell four- momentum with p 2 = m 2 . From this ex- 
pression we read off that the invariant normalization of a momentum eigenstate 
\P + , Pj_) = \P) is given by 

(P\K) = lQn 3 P + 5 3 (P - K) . (6.10) 

We already know that the bare Fock vacuum is an eigenstate of the interacting 
Hamiltonian. It thus serves as an appropriate ground state on top of which we 
can build a reasonable Fock expansion. If we specialize immediately to the case 
of QCD, we are left with the Fock basis states 

|0>, 

\qq : = b\k 1: ai)d)(k 2 , a 2 ) |0) , (6.11) 

\qqg : k h ai) = b j (k 1 ,a 1 )d j (k2,a 2 )a t (k 3 ,a 3 )\0) , 

: (6.12) 

In these expressions, tf, and a) create quarks q, antiquarks q and gluons g 
with momenta ki from the trivial vacuum |0). The <X{ denote all other relevant 
quantum numbers, like helicity, flavor and color. 

In a more condensed notation we can thus describe, say, a pion with momen- 
tum P = (P + ,P±), as 

n,Ai i » 1 

where we have suppressed all discrete quantum numbers apart from the helicities 
Aj. The integration measure takes care of the constraints Q3.39 ) which the relative 
momenta in each Fock state (labeled by n) have to obey, 

\\dxi = Y[dXi5(l -^Xj) , (6.14) 

i i j 

l[d 2 k ±l = 16vr 3 J]rf 2 A; ± ^ 2 (^k ±J ) . (6.15) 

i i j 
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As a mnemonic rule, we note that any measure factor d 2 k±i is always accompanied 
by l/16vr 3 . 

The most important quantities in ( 6.13Q are the light-cone wave functions 



(n : XiP + , XiP± + k ±i , \i\n(P)) 



(6.16) 



which are the amplitudes to find n constituents with relative momenta p\ = XiP + , 
p ±i = £jPj_ + k±i and helicities A, in the pion. Due to the separation properties 
of the light-cone Hamiltonian the wave functions do not depend on the total 
momentum P of the pion. Applying ( |6.10| ) to the pion state ( |6.13| ), we obtain 
the normalization condition 



e / 

n,Xi i 



d 2 k±i 
16tt 3 



\ip n / n (xi,k ±i , Ai)p 



(6.17) 



The light-cone bound-state equation for the pion is a straightforward generaliza- 
tion of (IE! 



(M 2 

(M 2 - M, 



M 2 -) 

11' 



119 ) 



(qq\ir) 



(qq\W\qq) (qq\W\qqg) ... 
(qqg\W\qq) (qqg\W\qqg) ... 



(qq\n) 
(6. 





If a constituent picture for the pion were true, its wave function would be en- 
tirely given by the projection (qq\ir) onto the valence state. All the higher Fock 
contributions would vanish and the unitarity sum ( |6.17|) would simply reduce to 



E 

AA' 



dx 



d 2 k A 
167T 



A I 2 



\i>qq/n(x,k±,\,\')\ 



1 . 



(6.19) 



We will later discuss a model where this is indeed a good approximation to reality. 



6.2 Properties 

Let us rewrite the light-cone bound-state equation ( |6.8|) by collecting all light- 
cone wave functions ip n = (n\ip) into a vector 



From this expression it is obvious that all light-cone wave functions tend to vanish 
whenever the denominator 

e = M 2 — Mq = M 2 — (J> t ) 2 = M 2 ~Y1 ^ + m * (6-21) 
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becomes very large. We know from Section 3.3 that this quantity measures how 
far off energy shell the total system, i.e. the bound-state is, 

P--J2p; = e/P + . (6.22) 



For this reason, e is sometimes called the 'off-shellness' [55, 0. We thus learn 



from ( |6.20| ) that there is only a small overlap of the bound-state with Fock states 



that are far off shell. This implies the limiting behavior 

ip(xi, kj_j, Xi) — > for Xi — * , Ar[j — > oo . (6.23) 

These boundary conditions are related to the self-adjointness of the light-cone 
Hamiltonian and to the finiteness of its matrix elements. Analogous criteria have 
been used recently to relate wave functions for different Fock states n (via "ladder 
relations" , [jSgfl ) and to analyze the divergence structure of light-cone perturbation 
theory 

Suppressing spin, flavor and color degrees of freedom, a light-cone wave func- 
tion will be a scalar function </>(xj, kj_j) of the parameter e. This is used for 
building models, the most common one being to assume a Gaussian behavior 



<P(x l ,k ±l ) = Nexp(-\e\/(3 2 ) , (6.24) 

where (3 measures the size of the wave function in momentum space. Note, 
however, that a Gaussian ansatz is in conflict with perturbation theory which 
is the appropriate tool to study the high-k^ behavior and indicates a power 
decay of the renormalized wave function (up to possible logarithms, ||82]| ). For 
the unrenormalized wave functions the boundary conditions ( |6.23| ) are violated 
unless one uses a cutoff as a regulator. 

As the off-shellness e is the most important quantity characterizing a light- 
cone wave function let us have a closer look by specializing to the simplest possible 
system, namely two bound particles of equal mass. One can think of this, for 
instance, as the valence wave function of the pion. The off-shellness becomes 



M 



2 kj_ + m 2 



x(l — x) x(l — x) 




M\x-\Y + +k 2 ± . (6.25) 



The second term on the right-hand-side is positive because, for a bound state, 
the binding energy, 

E = M — 2m , (6.26) 

is negative so that 2m > M. As a result, the off-shellness is always negative. 
Only for free particles it is zero, because all momentum components (including 
the energy) sum up to the total momentum. In this case, each individual term 
in (|6.25|) vanishes, 

M = 2m, x = | , k ± = 0. (6.27) 
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It follows that the light-cone wave function of a free two-particle system (with 
equal-mass constituents) is of the form 

4>(x, kj.) ~ S(x - 1/2) 5 2 (k ± ) . (6.28) 

One expects that for weak binding, in particular in the non-relativistic case, the 
wave functions will not look too differently, and thus will be highly peaked around 
x = 1/2 (in the equal mass case) and kj_ = 0. Let us check this explicitly. We go 
to the particle rest frame with P = or P + = P~ = M and = 0, implying 
p ±i = k±i. In this frame, the non-relativistic limit is defined by the following 
inequalities for the constituent masses and momenta (in ordinary instant-form 
coordinates) , 



p 2 

p° -mi ~ -p- < IpJ < mi . (6.29) 

IVTLi 

The prototype systems in this class are of hydrogen type where we have for 
binding energy and r.m.s. momentum, 

(p) = ma , (6.31) 

with a = 1/137 the fine structure constant and m the reduced mass. In this case, 
the hierarchy (|6.29|) becomes 

a 2 

— < a < 1 , (6.32) 

which is fulfilled to a very good extent in view of the smallness of a. 
Consider now the longitudinal momentum of the i th constituent, 

Pt =P°i +P 3 i -mi+ ^ + p\ = %iP + = XiM . (6.33) 

We thus find that we should replace pf in instant-form non-relativistic wave 
functions by 

Pi = XiM - m { , (6.34) 

where we neglect terms of order p?/mj . Let us analyze the consequences for the 
off-shellness. The latter is in ordinary coordinates 

e = M 2 — Mq = (M + ]>>°)(M - . (6.35) 

i i 

We thus need 

Erf = Eli: = <^ 
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with E = M — J2 m ii so that the off-shellness ( |6.35| ) becomes 

k 2 ±l + {Mxi - mi) 2 



e ~ 2M(E - ^ P 2 i/ 2m i) - 2M 



2rrij 



(6.37) 



where we have performed the replacement ( |6.34| ) in the second identity. The 
light-cone wave functions will be peaked where the off-shellness is small, that is, 
for 

x i = m i /M 1 and kj_j = , (6.38) 

as expected from the non-interacting case. 

As an explicit example we consider the ground state wave function of positro- 
nium, which in ordinary momentum space is given by the Fourier transform of 
ijj(r) = Nexp(-mar). Performing (|6.34j) once more and normalizing to one, we 
obtain 

87T7 



k\ + (xM - m e ) 2 + Y 



2 ' 



(6.39) 



where m e is the electron mass and 7 = ma the r.m.s. momentum. The result 
( |6.39| ) is valid for small momenta, i.e. when k"j_, (xM — m e ) 2 m 2 . It is obvious 
from (|6.39|) that the positronium wave function is sharply peaked around x = 
m e /M ~ 1/2 and k\ = 0. 



Part II 
Applications 



Chapter 7 



Zero Modes and the Light-Cone 
Vacuum 



Our first application will not (yet) deal with light-cone wave functions but rather 
with the issue of Chapter 5, the triviality of the light-cone vacuum. We will, 
however, be concerned with a non-trivial example of an interacting field theory, 
namely </> 4 -theory in 1+1 dimensions. As before we will address the question how 
a non-vanishing vacuum expectation value of the field can arise in the light- 
cone vacuum. As a further motivation we note that the dynamical system under 
consideration may be viewed as a low-dimensional caricature of the Higgs sector 
of the standard model which is also given by a 4 -type scalar field theory. 

Historically, it was Weinberg in 1966, who realized that within "old-fashioned" 
Hamiltonian perturbation theory the infinite-momentum limit of many diagrams, 



in particular vacuum diagrams, is vanishing f94fl . This limit was later on shown to 
be (at least perturbatively) equivalent to light-cone quantization [78 , [T7j| . Within 
the context of current algebra |)5], ^6], [J?] the success of this limit was traced 



to the fact that light-like charges always annihilate the vacuum, irrespective of 



whether they are conserved or not [pq]. Thus, Coleman's theorem ||98||, "the 



symmetry of the vacuum is the symmetry of the world" , does not apply PS| . As 
we have seen in Chapter 5, the Fock vacuum, i.e. the ground state of the free 
Hamiltonian, is stable under interaction. In the same manner as the light-like 
charges, the fully interacting Hamiltonian annihilates the vacuum, the formal 
reason being that both quantities conserve longitudinal momentum and thus are 
of the form ( |5.9| ). 

With the advent of QCD as the theory of strong interactions, however, people 
began to feel somewhat uneasy. There was (and is) growing evidence, that many 
of the phenomenological aspects of hadron physics, like confinement and chiral 
symmetry breaking are related to the non-trivial features of the QCD vacuum 
(within standard instant-form quantization). Let us only mention features like 
quark- and gluon condensates, instantons, monopole condensation etc., which all 
indicate that the vacuum is densely populated by non-trivial quantum fluctua- 
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tions, which furthermore are not accessible to perturbation theory. 

The concern that arose at these points can be put into the question: can 
the existence of these large vacuum fluctuations be reconciled with the triviality 
of the light-front vacuum? For QCD, the answer to this question is not (yet) 
known. For model field theories, some answers have already been found. Not 
surprisingly, it is the zero modes that are to some extent the carriers of non- 
trivial vacuum properties. One generally distinguishes between constrained and 
unconstrained zero modes. The latter typically arise in (low-dimensional) gauge 
theories ||100| , |101| , |102| , |103| , |104j , |105| , |106|| . The former we have already met in 
the discussion of the constant source (Subsection 5.2.1). For 4 -theory in 1+1 
dimensions, it has been shown that the zero modes are responsible for spontaneous 
symmetry breaking of the reflection symmetry — > — ||107| , |108| , |109| , |110| , |1 1 1 
Let us review the basic features of what is going on 1 1 12 



7.1 The Constrained Zero Mode 

We recall that the Lagrangian of a scalar field theory in 1+1 dimensions, written 
in terms of light-cone coordinates, cf. ( f4.95| ), 



L[<f),d-(f>) = \ J dx-\d + cj)d-(t)-U[(t)\ , (7.1) 



is linear in the velocity d~<p. According to the discussion of Chapter 4, the 
elementary bracket between the fields is essentially given by the inverse of the 
spatial derivative, d + = 2d/dx~. In order to uniquely define this inverse, one 
has to specify the domain of d~ and boundary conditions. As before, we enclose 
our spatial variable x~ in a box, — L < x < L, and impose periodic boundary 
conditions on our fields. We know that in this case the operator d + has zero 
modes, namely all spatially constant functions. To deal with this aspect, we split 
our field <fi into a zero mode u and its complement <p, 

4>(x + ,x~) = oj(x + ) + (p(x + , x~) , (7.2) 

L 

u(x + ) = I dx~4>(x + ,x~) , (7.3) 
such that ip does not have a zero mode, 



-L 



yW,-) = 0. (7.4) 

—L 

The Lagrangian ( pM|) can then be rewritten as 

L 

L[<p,u] = J dx-±(p(-2d_)<p-H[<p,w] . (7.5) 



—L 
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For brevity, we use the derivative <9_ = d/dx~ and (p = d(p/dx + . Note that uj 
does not appear in the kinetic term. The Hamiltonian H = P~ coincides with 
the potential U from ( |7.1|) after the replacement ( [7.2| ) has been performed. We 
thus find that the zero mode u is constrained by its 'equation of motion', 



e = 



SH_ 

5uj 



. 



(7.6) 



The basic bracket can be read off from the kinetic term in the Lagrangian (a la 
Schwinger or Faddeev/Jackiw, cf. Chapter 4), 



1 i 



~sgn(af 



y-) 



x -y- 



2L 



7.7) 

which, of course, coincides with ( [4.59 ) upon quantization via the correspondence 
principle. 

As is well known, not all classical observables can be quantized unambiguously 
due to possible operator ordering problems | |113| , |114| , |115| , |116|| . Such problems 
do not arise for the field ip and the bracket (|7. 7|) , where the field-independent 
right-hand-side leads to a c- number commutator. The constraint (|7.6|) , however, 
implies a functional dependence of the zero mode u on ip and thus a non-vanishing 
commutator of uj with (p. This can be explicitly verified by calculating the asso- 
ciated Dirac bracket within the Dirac-Bergmann algorithm ||117| , |118j , |119|| . For 
the quantum theory, it results in an ordering ambiguity with respect to to and 
if [ IP7 , |110| |. Therefore, a definite ordering has to be prescribed. We chose 



Weyl (or symmetric) ordering ||113| , |120|| which is explicitly hermitian. Using this 



prescription, one finds the quantum Hamiltonian for light-cone 0^ +1 -theory, 



H = dx 



L 

d,X~ -m 2 CU 2 + — (u 4 ' + COLf 3 + ipujip 2 + ip 2 uj(p + Lf 3 LU + 

-L 

+ uj 2 (p 2 + ip 2 u 2 + uxpujip + ipuipuo + u{p 2 u + {pu) 2 ipj 



■ (7-8) 



Note that we have chosen the sign of the mass term(s) in such a way that there 
is no spontaneous symmetry breaking at tree level as the classical potential does 
not have any turning points. With this Hamiltonian, equation ( |7.6j ) for the 
constrained zero mode u reads explicitly 



1j 

= -z— = m uj + —uj + — — ax \tp + tp uj + (pujip + ujip = . 
ouj 31 31 2L J L J 



(7.9) 
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This is nothing but the zero mode of the Euler-Lagrange equation of motion for 
the total field decomposed into u and <p |63], |121| , |110|| . The remainder of this 
chapter is concerned with different approaches to solve this equation for uj. 



7.2 Perturbative Solution 

To obtain a perturbative solution for uj we expand it in a power series in A, 

oo 

uj = ^\ n u n (7.10) 

n=0 



Inserting this into ( |7.9|) determines the coefficients u n recursively. For the first 
three we find 



u = 



1 1 



Qm 2 2L 



dx <p 3 (x) 



(7.11) 
(7.12) 



1 1 

36m 4 (2L~y 



dx dy~ 



V 2 (x)uj 3 {y)+uj{x)uj 3 {y)uj{x) + oj 3 (y) if 2 (x) (7.13) 



All higher orders may be obtained similarly. Unfortunately, however, we have 
not been able to find a closed formula for u n which would allow for summing up 
the whole series ( |7.10|) . 

If we expand the quantum field <p in terms of Fock operators, cf. ( |4.107|) , 



n=l 



'Ann - 



a n e iknX 12 + ale iknX 12 



(7.14) 



with the discretized longitudinal momentum fc+ = 2im/L, one notes that the 
vacuum expectation value of uj is zero to all orders in A, since u n contains an odd 
number of Fock operators, i.e. (0|a> n |0) = 0. Thus 



|0> = 5^A»<0|w B |0> = 0. 



(7.15) 



n=0 



This seems to imply that a non-vanishing vacuum expectation value for u can 
only arise non-perturbatively and must be non-analytic in the coupling A. We 
will discuss this issue in the next sections and continue for the time being within 
the framework of perturbation theory. In particular, we will study the effect 



7. 2. PERTURBATIVE SOL UTION 



S3 



of the zero mode u on the mass renormalization. To this end, we split up the 
Hamiltonian H from (|7.8|) into two pieces, 



H = H + H UJ , (7.16) 

where H is independent of lu, 



L 

^ 2 2 i ^ , „ J 



H = I dx [ -mV + ) , (7-17) 



and Hu is the a;-dependent interaction given by the last two lines in (|7.8| ). H w 
can be simplified using the constraint equation (|7.9|), 

L 

= J dx" (—u 4 + (pujip 2 + (p 2 uj(p + (puj 2 ip — uoip 2 uo) . (7.18) 

— L 

Because u is of order A, the zero mode dependent part H w is of order A 2 . Explic- 
itly, we find, neglecting terms of order A 3 , 



H w = — J dx" {ipUJip 2 + Lp 2 UHp) 



L 

A 2 1 



144m 2 2L 



dx"dy" [ V (x) V 3 (y)ip 2 (x) + V 2 (x)ip 3 (yMx)] ,(7.19) 



where we have used the first order term ( [7.12 ) . This induces a mass shift of order 
A 2 which is given by ||119|| (see also | 108|| ) 

5m 2 n = ^((n\H u \n) - (0|#J0>) . (7.20) 

Here, |n) = c4|0) denotes a one-particle state of longitudinal momentum fc+ = 
27m/ L. We have subtracted the constant vacuum energy, (0|i?^ |0), which di- 
verges linearly with the volume. This is sufficient to render the mass shift finite. 
After inserting the Fock-expansion ( |7.14| ) into ( [7.191 ) one obtains for ( |7.20| ) 

A 2 L lr^ 1 .A 1 i2t™ 



Sm 2 = —- — _____ h4^ 



6m 2 (47r) 3 n L — — (n — m)m — — (n + m)m 

v ' m=l v ' m=l v ' 

A 2 



6m 2 16/T 



V [37 + *(n) + 2tf(l + n)l < , (7.21) 
7r z r2 L J 
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where \l/ denotes the Digamma-function | 122 |, and 7 = — = 0.577216... 
is Euler's constant. The question now is, whether the result ( |7.21|) is a finite 
size effect which vanishes in the infinite- volume limit. The latter is obtained by 
replacing 

2ttN/L 

lim — / dk + f(k + L/2n) (7.22) 
l^oo 2n J 

n=1 2n/L 

in such a way that /c+ approaches a finite limit k + , i.e. n, L — > 00 with n/L finite. 
For the case at hand, this amounts to replacing the Digamma-function by its 



A' 



asymptotics ||122| 

ty(z) ~ log(z) 



1 

2z 



+ 0(l/z 2 



00 in |arg(z)| < n) 



Denoting the infinite-volume limit of the mass shift (|7.20| ) by 5 



obtain 



5m 



X 2 



171' 



lim — 



\og{k + L/2-n) + 7 + 0{l/k + L) 



m 



0. 



(7.23) 
we finally 

(7.24) 



The vanishing of this expression implies that the zero mode induced second order 
mass shift 5m 2 n is indeed a finite size effect. We do not have a general proof 
that this is also true for higher orders in A. However, it seems to be plausible 
that a single mode like to constitutes a "set of measure zero" within the infinite 
number of modes as long as one applies perturbation theory. From the theory of 
condensation, however, it is well known that single modes, especially those with 
vanishing momentum, can significantly alter the perturbative results. To analyze 
this possibility one clearly has to use non-perturbative methods. This will be 
done in the next sections. 



7.3 Mean-Field Ansatz 

As </4 +1 -theory is super-renormalizable, the number of divergent diagrams is fi- 
nite, namely one: the tadpole resulting from a self-contraction of two fields at the 
same space-time point. Therefore, renormalization can, at least within pertur- 
bation theory, be done via normal-ordering, which is nothing but making use of 
Wick's theorem: expanding the appearing powers of in a sum of normal-ordered 
terms with more and more self-contractions, one separates the finite term (with no 
contraction) from the divergent ones (with at least one contraction). The latter 
terms then are just the negative of the required counter-terms. For conventional 
</> 4 -theory one finds 

+ = ^m 2 (:0 2 :+T) + A( ;0 4 :+6T:02;+3r2 j 
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where 

T=(0 2 > = / - 1 (7.26) 
J 4tt Vfc 2 + m 2 

formally denotes the logarithmically divergent tadpole contribution (in d = 1+1) 
which coincides with the vacuum expectation value of <fi 2 or a self-contraction of 
the field. 

It is obvious that only mass and vacuum energy get renormalized. The renor- 
malized Hamiltonian is obtained by adding the counter-terms 

SU = --T : 2 : --m 2 T - —T 2 . (7.27) 
4^28 V ' 

If the vacuum expectation value of <fi 2 is taken in a Fock vacuum corresponding 
to the bare mass m, the latter coincides with the renormalized mass to order A. 

This rather trivial renormalization procedure cannot be straightforwardly ex- 
tended to light-cone field theory, simply because we do not know the Hamiltonian! 
The zero mode u is a complicated functional of the Fock operators a n , a* , which 
has to be found from the constraint (|7.9|) before we can normal-order. There is, 
however, a way around this obstacle, such that an exact knowledge of to is not 
needed for renormalization. 

To this end we use the following general ansatz for uj, 



n>0 m,n>0 m,n>0 

Ulmn 0'i +m+n dl(l m a n + ^Imn a l a m a n a l+m+n + 

l,m,n>0 l,m,n>0 

+ 8kM,m+nWklmn a k a 'l a >m a n + • • ■ j (7.28) 



k,l,m,n>0 



where u , u> n , ... are c-numbers to be determined. Obviously, u is the vacuum 
expectation value of the operator uj, 

coo = (0|u;|0} , (7.29) 

whereas the coefficients ui n etc. contribute to matrix elements of uj in higher Fock 
sectors. The ansatz ( 7.28Q is hermitian and takes care of the fact that uj cannot 



transfer any momentum. Thus, it is a discrete analogue of ( |5.9|) and can be 
understood as a Wick expansion written in the opposite of the usual order: the 
first term ujq is the sum of all contractions, the second the sum of all contractions 
but one and so on. Accordingly, each individual term in the expansion is a normal- 
ordered operator. Plugging the ansatz ( [7.28j ) into H w one notes that the latter 
may be viewed as an effective Hamiltonian where additional operators ip n , n > 2 
have been induced by the zero mode. This is reminiscent of a renormalization 
group analysis in which the same operators, being all marginal, are generated 



86 



CHAPTER 7. ZERO MODES AND THE LIGHT-CONE VACUUM 



by the renormalization group flow ||123| , |83f . In what follows, we will, however, 
truncate the ansatz (|7.28j ) after the second term (~ a* a n ) which only introduces 
additional if 2 contributions and thus corresponds to a mean- field treatment. 

Inserting the ansatz ( [7.281) into ( [7.8|) and ( [7.9|) we obtain for the constraint 
and the Hamiltonian 



9 = 9o + ^2 °nala n + ... , 

n>0 

H = H + ^2 H n(4i a n + ■■■ 



(7.30) 
(7.31) 



n>0 



where, in accordance with the truncation of our ansatz ( [7.28| ), we have omitted 
terms containing more than two Fock operators. The coefficients H , H n and 9 , 
9 n are functions of u and u> n . Thus, ( |7.31| ) is an effective one-body or mean-field 
Hamiltonian describing the influence of the zero mode u. Explicitly, one finds for 
the coefficients of the constraint, 



0o 



A. 



A 



-T ) u + — \ LO, 



m 2 + —T ) uj. 



3! 

A 
3! 



E 

n>0 



4irn 



n 11 6CC?Q 6C(J n 

u n + 3^ + + — + — 



(7.32) 
(7.33) 



and of the Hamiltonian (scaled by 2L) 



Ho 
2L 



m 2 + ^T)^H 



A 
4! 



n>0 



n>0 



4im 



+ — T + —T 2 , 
2 8 

^ = -(m 2 + -T 
2L 2 V 2 



+ A ( (o; Q + u n y - ^ + — ( Wo + co n ) 2 + J" — ) ■ (7-35) 

ZLI V irn Vim ^— » irrn / 



(7.34) 



cj„ + 2u; ^n + 



27m 



+ 



4! 



27m 



m>0 



irm 



T now denotes the light-cone tadpole (in discretized form) 

^ 1 ^ 4vm 



n>0 



(7.36) 



which is mass independent in contrast to ( |7.26|) (unless one regularizes with a mass 
dependent cutoff, cf. Chapter 9). Note that the one-particle matrix elements of 
9 and H are given as the sum of two coefficients, 

(n\9\n) = 9 + 9 n = ( m 2 + —T] (u + u n ) 
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+ 



A 
3! 



(u + u n ) +6- 



47m 



fe>0 



ink 



(7.37) 



and 



(n\H/2L\n) = (H + H n )/2L 



+ "77 



A, 



1 



m 2 + -T) (u + uj n y + -[m 2 + -T 



A, 



2nn 



{u + u n f + (12wjj + 24cj cu n + Uuj 2 n )- — + 

47m 



m 2 m A ~ 
— T + — T . 
2 8 



+ 



+ 



(7.38) 



From the divergence structure above it is clear that all coefficients can be made 
finite by adding the counter-term 



6H/2L = --T y —a{a n - -Tu 2 - -m 2 T - -T 2 
1 4 ^27rn n 4 2 

n=l 



8" 



(7.39) 



which can be obtained from ( |7.27| ) by integrating over x~ and decomposing the 
field. The renormalization is thus standard, i.e. performed by normal-ordering 
and formally achieved by setting T = in the expressions above. 
It is convenient to rescale the coefficients 



LU 



U), 



UJQ 



'An 



and define a dimensionless coupling g as 

A 



9 



247rm 2 



> 



such that ( |7.32| - |7.35| ) become 



H /2L 
HJ2L 



+ 



m 

47T 

2 

m 
Air 

9 
4 



1 2 , 9 



n=l 



4 \ u ^ n ^ n 



OO o 

a;: 



n=l 



(7.40) 
(7.41) 

(7.42) 
(7.43) 



u n + 2u u n + — J + 
12 



fc=i 



n f--' k n 



(7.44) 
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and 



n=l 

\ n 
This system of equations has a trivial solution 



. 



(7.45) 
(7.46) 

(7.47) 



corresponding to the symmetric phase with vanishing vacuum expectation value 
of the field. Note further, that uo n = implies ujq = 0, which means that it 
is the operator part of u that drives a possible phase transition. In the latter 
case, where u n ^ 0, there are non-trivial solutions of (7.45, 7.46| ) which, however, 
cannot be obtained exactly. If we assume that there is a critical coupling where 
the field starts to develop a vacuum expectation value, then, very close to this 
coupling, ojq should be small so that we can expand u n as a power series in u 

ram, 

u n = a n (g) luq + (3 n (g) u% + O(u 5 ) . (7.48) 



Plugging this into ( |7.46| ), the expansion coefficients are determined recursively 

as, 

6# 



a,, 



3n 



n + 6g 
ng 

n + 6g 
a i 6g 



(7.49) 



a n (3 + 3a n + a\ 
18<? 



6 \n + 6g 



+ 



n + 6g \n + 6g 



(7.50) 



Inserting u n into (|7.45|) , one obtains ujq as a function of the coupling g via 



a (g) + a 2 (g)u 



n>0 \ n>0 



UJn 



. 



(7.51) 



The coefficients ao and a 2 can be evaluated analytically with the result 

a (g) = 1- g[V(l + 6g) +j] , (7.52) 



g 



1 +18/ [^ (1) (1 + 6g) + 3^ (2) (1 + 6g) +2g 2 ¥ 3 \l + 6g)] (7.53) 



where the are Polygamma functions [122]. As (3 n and therefore a 2 are posi- 
tive, equation (|7.51|) develops two real solutions for uq (with opposite sign) if 



ao(g) < . 



(7.54) 



7.3. MEAN-FIELD ANSATZ 



89 



The critical coupling g c is determined if equality holds. Using the explicit form 
of do from ( |7.52|) , one finds 

a (g c ) = 1 - g c [*(l + Qg c ) + j] = , (7.55) 

with a numerical value for the critical coupling of 

g c = 0.53070059 . . . , (7.56) 

which in terms of the original parameters is (cf. ||124|| ) 

A c = 2Aixg c m 2 ~ 40.0 m 2 . (7.57) 

Above this coupling, the vacuum expectation value u is non-vanishing and ac- 
quires one of two possible signs so that the reflection symmetry is spontaneously 
broken. 

It is very illuminating to visualize the discussion above in terms of an effective 
potential V e s [|107j |125|| . To this end we interpret (|7.51|) , which determines the 
vacuum expectation value u , as the equation localizing the minima of V c g. The 
latter is thus obtained via integration, 

V cS [Q] = \a»(g) Q 2 + \a 2 {g) ft 4 . (7.58) 

The integration constant has been chosen such that V^ff[f2 = 0] = 0. At the 
critical coupling, g c , the coefficient a of Q 2 changes its sign from positive to 
negative. Therefore, below g c , the effective potential has no turning points and 
its minimum is located at u = (symmetric phase). On the other hand, above 
g c , V e e develops two non-trivial minima at ±Uq (broken phase). As the order pa- 
rameter ujq changes continuously, the phase transition is second order. Obviously, 
Kff has precisely the form of a Landau- Ginzburg free energy which corroborates 
our mean-field interpretation. The final confirmation is provided by the calcula- 
tion of the critical exponents. In a thermodynamic/magnetic analogy we view g 
as a temperature, Q as a magnetization, and dV^/dVL as an external magnetic 
field. The coefficient is then an inverse susceptibility which, near the critical 
coupling, should behave like ao ~ (g — gc) 1 ■ (Evidently, the exponent 7 has noth- 
ing to do with Euler's constant.) From ( |7.52|) and ( |7.55|) we obtain via Taylor 
expansion around g c , 

ao(g) -(9- 9c) a' (g c ) , (7.59) 

yielding the expected mean field result, 7 = 1. The order parameter goes like 
- gc) 13 . Using ( |7.51| ) and ( |7.59| ), we find 

wb = ± [ - a (g)/a 2 (g)] 1/2 ~±(jg- g c ) 1/2 [ - a (g c )/a 2 (g c )] 1/2 , (7.60) 

i.e. a mean field exponent (3 = 1/2. The dependence of the rescaled vacuum 
expectation value u on g is plotted in Fig. |7Tl) . From both (|7.60|) and the figure 
it is evident that uo is non-analytic in the coupling g. 
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The last exponent, 5, is denned through the behavior of the magnetic field at 
g c , dV e s/dQ ~ Q 5 . We calculate 



dV cS /dQ\ g=gc = a 2 (g c )n 3 , (7.61) 
so that 5 = 3. As a cross-check, we note that our critical exponents obey the 

(7.62) 



correct scaling law [|123|| , 



1 + 



2 





At this point it is instructive to make contact with perturbation theory. For 
this purpose, we expand a (g) as given in ( |7.52| ) in powers of g using the series 
representation ||122|| 



*(l + z) = - 7 + j^(-l)»C(n)« 



n-l 



(7.63) 



n=2 



from which we read off that + z) is the generating function of the sequence 
C(n), C being the Riemann zeta function. The convergence region of the right- 
hand-side of ( |7.63| ) is \z\ < 1. For a (g) this yields, 



ao(g) = 1-6 ((2) g 2 + 36 C(3) g 3 + 0(g 4 



(7.64) 
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which converges for g < 1/6. Keeping just the first two terms, one can find a 
critical coupling perturbatively ||107| , |125| , |126j| via 



a (^ P ) = l-6C(2)^ =0 



(7.65) 



leading to 

9cp = l/vr = 0.3183098... > 1/6 (7.66) 

which is beyond the radius of convergence of (|7.64| ). If we include the third order 
in g, the equation ao(g) = does no longer have a positive root like g cp . We 
conclude that in the present framework, perturbation theory is not applicable for 
studying the critical behavior, and the value (|7.66|) for g cp has to be abandoned. 
In other words, upon using perturbation theory one stays in the symmetric phase 
with u Q being zero, in agreement with the discussion of Section 7.2. 

Finally, from the one-particle energies ( |7.44| ), we calculate the mass-shift in- 
duced by the non-perturbative zero mode u as given by the ansatz (|7.28|) , 



5ml 



2-RTl 



m 



m 



n 



(u 2 n + 2u u n ) + 



no 

+ T 



12 



U + U n ) 4 — O>o H (^0 + ^n) 2 

n 



k=l 

Using the expression (|7.48| ) for u n this becomes (near the critical coupling 



(7.67) 



5ml 



3gm 2 u%t 1-2 



1 + 6g/n 



(*(1 + Gg) +7 



+ 



+ 6<7- 



n + 7g 
n + 6gY 



(7.68) 



In contrast to the perturbative result (|7.21| ), this expression is not vanishing in 
the continuum limit, 

5m 2 = 3gm 2 ul 1 - 2a Q {g) . (7.69) 

Due to ( |7.54|) and (|7.55| ), a (g) is a small (negative) quantity. The mass shift 
induced by to is thus positive. We want to emphasize that the zero mode ui has 
a non-trivial influence on the spectrum. 



7.4 The Tamm-Dancoff Approximation 
7.4.1 General Remarks 

The Tamm-Dancoff approximation was originally designed within the conven- 
tional formulation of quantum field theory [Q. The principle of the method was 
to enormously reduce the particle number to a finite (and small) one. Due to 
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the complicated many-body structure of the equal-time vacuum, however, the 
approximation failed to lead to quantitative results and was abandoned there- 
after. Although it has been noted rather early that light-cone field theory might 
be better suited for a Tamm-Dancoff approximation 



127 1 



it was not until re- 
cently that people began to start systematic investigations ||128|| . The hope was 
(and is), of course, that many of the problems of the original Tamm-Dancoff 
approximation can be avoided due to the simplicity of the light-cone vacuum. 
Specifically, people first tried to calculate bound states for 1+1 dimensional field 
theories where the method turned out to work extremely well (see Chapter 8 and 
127 , |129| , |130| , |131| , |132|| ). Considerably more effort is needed in 3+1 dimensions 
133| , |134|| , the essential problem to solve being renormalization within a Hamil- 

[i36i, nm , nn |i39 



tonian formulation [113 



as the structure of counter-terms 
as well as their required number is unclear. Furthermore, dynamical symmetries 
like rotational invariance or parity are violated in general []14U[| . Nevertheless, the 
first phenomenological applications have appeared recently ||141| , |142| , |143j . 



In the following we will use the model at hand to shed some light on the 
issue of renormalization, in particular on the construction of the counter-terms 
needed. These could, in principle, be different for different particle number sec- 
tors ("sector-dependent renormalization" ||128|| ). For a super-renormalizable the- 
ory like which conventionally is renormalizable by normal-ordering, this 
would be a rather undesired feature since it would complicate the renormaliza- 
tion procedure enormously. Therefore, in the following, we will try to keep the 
renormalization as simple as possible. To this end, we attempt to incorporate 
the normal-ordering prescription into the Tamm-Dancoff approximation. First we 
need some definitions. The projection operators on the lowest particle number 
states are 



Po 
Pi 



|0)(0| 
^2\n)(n\ 



n>0 



Pi = \ Yl k m >( 



m, n\ 



(7.70) 
(7.71) 

(7.72) 



n,m>0 



P 



N 



\ni,...,n N )(m,...,n N \ 



ni,...,n N >0 



(7.73) 



with the Fock basis states defined in (|6.1| ). Most important to us will be the 
Tamm-Dancoff projector 



(7.74) 



a=l 



which projects onto the direct sum of all sectors with particle number less than or 
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equal to N. Before we perform any detailed calculation, let us make a few remarks 
about the relativistic invariance of the Tamm-Dancoff approximation [|144j| . 

Let V denote the Poincare group. In d space-time dimensions the number of 
Poincare generators is 

dimP ^+^ D. (7.75) 

If the Poincare generators are Go, G\i ■ ■ ■ Gd-i> the Tamm-Dancoff approximation 
is relativistically invariant if 

[PnG^PnG.Pn] = P N [G t ,G fc ]P N . (7.76) 

It is easy to see, that this expression can only hold if at most one of the Poincare 
generators fails to commute with Pfj. So we must have, e.g. 

[G ,Pn] + 0, (7.77) 
[G fc ,P N ] = 0, k = l,...D-l. (7.78) 

The last identity implies that the Gk conserves particle number and therefore 
must be kinematical, i.e. interaction independent. Thus, Go must be the only 
dynamical generator. Clearly, this can only happen if the dimension of the sta- 
bility group (of kinematical generators) is D — 1. There is only one single case 
where this requirement is met, namely the front form in d = 1+1, which is the 
basis of the present discussion. Here, according to ( |7.75| ), there are only three 
generators; the dynamical generator is Go = P~ , the light-cone energy; the kine- 
matical generators are the momentum P + and the (longitudinal) boost M + ~ . 
We think that this unique feature is one of the reasons, why the Tamm-Dancoff 
approximation (or more generally: Fock space truncation methods) work so well 
for light-cone field theories in d — 1+1. 

Encouraged by this observations we continue and analyze the impact of the 
Tamm-Dancoff approximation on the quantum nature of the scalar field. This 
will be important for the issue of normal-ordering. We imagine that any oper- 
ator O can be built from the elementary Fock operators a n , a' n , no matter how 
complicated its form may be. Thus, we define the iV-particle Tamm-Dancoff 
approximation by the replacement 

a n -> P n ci„Pn , (7.79) 
4 -> P N 4P N . (7.80) 

In this way, however, one is mutilating the quantum structure of the theory. This 
can be seen by calculating the iV-particle Tamm-Dancoff approximation of the 
elementary commutator, 

[a m ,al} N = [PnOotPn, PnOhPn] • (7.81) 
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A somewhat tedious calculation leads to the final result ||112 

[o-mi a V\U = ^mnPlM-l ~ Pn^Ii^Pn • (7.82) 

Taking matrix elements of this expression one readily sees that the correct re- 
sult (S mn ) for the commutator within iV-particle Tamm-Dancoff approximation 
is reproduced only up to the (iV-l)-particle sector. The additional term on the 
right-hand-side of ( |7.82[) is acting in the iV-particle subspace only. In other words, 



matrix elements of expressions involving the elementary commutator, which are 
calculated within iV-particle Tamm-Dancoff approximation, should not be trusted 
beyond the (iV-l)-particle sector. This will be relevant for the problem of normal- 
ordering to be discussed shortly. 

To be a little bit more explicit, we list the lowest order expressions for 



[flm,4]i = $mn 1 0)(0 1 - \n)(m\ , (7.83) 



[ a m ; ^]2 — S mn 



\0)(0\ + ^\l)(l\\-^2\l,n)(m,l\. (7.84) 

l>0 l>0 

What are now the implications of all that for the renormalization, in particular 
mass renormalization? As can be seen from ( |7.25| ) and ( |7.27| ), the latter is encoded 
in the normal-ordering prescription of the expression 

L L 

J dx~ip 2 (x) = J dx~ : <p(x) 2 : +T 

-L -L 

n>0 ra>0 

Normal-ordering thus amounts to splitting off the divergent tadpole contribution 
T, which, in the Fock space language, is given by an elementary commutator (or 
contraction). Thus the remarks above, leading to ( [7.821) , apply. Let us calculate 
the iV-particle Tamm-Dancoff approximation of ( [7.85| ), 

L 

1 / , O , s iVPTDA V-^\ 1 „ + „ „ „ „ 

— / dx-if\x) ~ = 2Z^P N -l4PN-ianPN-l + PN-lT- 
—L n>0 7771 

- Y^^—PNalanPN . (7.86) 

n>0 

The same remarks as for (|7.82|) are in order. Matrix elements of the last expres- 
sion should not be expected to be consistent beyond the (N-l) -particle sector. 
Furthermore, one should note that, as if 2 is a one-body operator, there is only 
one commutator (or contraction) involved in the above normal-ordering. For a 
fc-body operator we therefore conjecture that its renormalization will be correct 
only up to the (N — A;)-particle sector. For example, in order to get the renor- 
malization of the two-body operator yj 4 correct up to the one-particle sector, a 
three-particle Tamm-Dancoff approximation will be needed. 



7.4. THE TAMM-DANCOFF APPROXIMATION 



95 



7.4.2 One-Particle Tamm-Dancoff Approximation 

The one-particle light-cone Tamm-Dancoff approximation is defined by the re- 
placements 



a n ~ P 1 a n P 1 = |0}(n| (7.87) 
PiaJ,Pi = |n)(0| . (7.88) 



al ~ 



From our general results we expect that within one-particle light-front Tamm- 
Dancoff approximation we will get a consistent renormalization of the one-body 
operator (p 2 in the vacuum (zero-particle) sector only. 

To solve for the zero mode u we make the following Tamm-Dancoff ansatz 

ujtd = co|0)(0| + c » H • ( 7 - 89 ) 

ra>0 

We similarly expand the constraint 6 and the Hamiltonian H 

9 ~ 9 \0)(0\+J2W(n\ , (7.90) 

n>0 

H ~ H \0)(0\ + J2Hn\n}(n\, (7.91) 



n>0 



where the bars simply indicate the distinction of the coefficients above from those 
of the mean field ansatz ( |7.30|) and ( |7.31|) . The coefficients can be found as 



functions of cq, c n upon inserting the ansatz ( |7.89| ) into the constraint (|7.9| ) and 
the Hamiltonian (|7.8| ) yielding 



and 



So= (^ + ^)c + A U + E ^ , (7 , 2) 

7 \ n>0 / 

°" = I - 2 + iL) + a ( c » + • (7J3) 



21 - K™ !+ i^ + 5( c:+ ^s + E^ 

v 7 \ n>0 n>0 

H n If , A \ 1 / A \ 1 



H, 



2L 2 \ 16vrn J n 2 \ 12 J Ann 

A ( 4 , CoC r , 



+ ^ K + T^ + fH- ( 7 - 95 ) 
4! V n 2nn Ann 1 
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At a first look, these expressions for the coefficients appear to be a disaster: 
the infinities in form of the tadpole T do not appear systematically, the mass 
gets renormalized in the vacuum sector only, but differently for the constraint 9 
and the Hamiltonian H . Both expressions differ from the standard expression 
m 2 + XT/2. There is also a divergent contribution from the (p 4 term to Hq which 
differs from the usual AT 2 /8 of (|7.25 ). As stated above, we do not believe the 



coefficients 9 n and H n to be correct within one-particle light-cone Tamm-Dancoff 
approximation. They will be discussed in the next subsection, when we go to 
higher order. 

The way to remedy the situation (for the coefficients 9q an d So) is the follow- 
ing. We insist, firstly, on the standard mass renormalization, m 2 + XT/2, however, 
according to our general discussion, in the vacuum sector only. Secondly, we do 
not assume that the coefficients Co and c n are independent, and use this freedom 
to redefine them in the following way 

c = ujq , (7.96) 
c n = u + u n . (7.97) 

Inserting this into (|7.92|) and ( [T.94|) one finds 




k = ^ 2 + 2 T J^ + 3T^ + E^j , (7-98) 



+ -7T T +77 T • ( 7 - 99 ) 
2 4! 



Remarkably, the simple redefinition (|7.97 ) has led to the desired results. The 



mass renormalization is standard and the same for 9q and H . The divergences 



thus can be made to vanish by adding the counter-term ( |7.39| ). Both equations 



(7.98) and ( |7.99| ) coincide with the lowest order results from the mean-field ansatz 



( |7.32| ) and ( |7.34| ) (up to the constant (^-contribution to H Q given by AT 2 /4! 
Note that there are no two-body (T 2 ) contributions to the constraint. This is 
obviously true to all orders, so the renormalization of 9 is slightly simpler than 
that of H, namely just mass renormalization. 

The coincidence with the mean field results is not accidental. If one calculates 



the lowest order matrix elements of the mean field ansatz ( |7.28| ), one finds 



uo = (0\u MF \0) = (0\u TD \0) = c , (7.100) 
lu + u n = (n\u MF \n) = (n\u TD \n) = c n . (7.101) 

Analogous relations hold for the matrix elements of the constraint and the Hamil- 
tonian, 

(O|0|O) = 9 = 9 , (7.102) 
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(0|if |0) = H = H , (7.103) 
(n\6\n) = 9 Q + 9 n = 9 n , (7.104) 
(n\H\n) = H + H n = H n . (7.105) 



Thus, after the redefinition (|7.97 ), the zero- and one-particle matrix elements 



of oo calculated within the mean field ansatz and Tamm-Dancoff approximation 
coincide. As the renormalization within mean field ansatz was conventional and 
straightforward, it is not too surprising that the behavior of the redefined Tamm- 
Dancoff approximation under renormalization gets improved. This will be an- 
other guideline in the following. 

7.4.3 Two-Particle Tamm-Dancoff Approximation 

If we now go one step further and include also two-particle states via 

a n ~ P 2 a n P 2 = |0)(n| + ^|m)(m,n| (7.106) 

m>0 

4 ~ P 2 4P 2 = |n)(0| + ^2 \n,m){m\ , (7.107) 



m>0 



we should further improve our renormalization program. We expect a consistent 
renormalization of ^-contributions in the vacuum- and one-particle sector, and 
of (/^-contributions in the vacuum sector. The extended ansatz for u becomes 



UJ TD — c |0>(0| + ^c»(n| + 



n>0 



+ | 2j Cmn l m ' n )( m + n \ + I ^2 C *mn\ m + n )( m i n \ + 
m,n>0 m,n>0 

\k,l)(m,n\. (7.108) 

k,l,m,n>0 

It is now very plausible (though we cannot prove it a priori) that a consistent 
renormalization requires a redefinition also of the coefficients c mn , c* mn) and Cki mn - 
To proceed as before, we would need the matrix elements 9q, 9 n , Ho, and H n with 
all two-particle contributions in order to just get a consistent renormalization up 
to the one-particle sector. This is very tedious and inefficient, as we are only 
interested in the divergent contributions from higher order terms to lower order 
matrix elements. Fortunately, there is an alternative: we simply demand that the 
matrix elements of ujtd and ujmf coincide also in the two-particle sector. This 
gives us the desired redefinitions, namely 

Cmn 2cJ mn , (7.109) 

C mn — ^mn 5 (7.110) 

Cklmn 

(7.111) 
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Thus, essentially, only the two-particle matrix elements cumn get redefined. Ex- 
pression (|7.108|) becomes 

Utd = w„|0)(0| + ^2{u + u ri )\n){n\ + 



+ E 



n>0 



,n)(m + n\+ V} u* mn \m + n)(m, n\ + 



m,n>0 



m,n>0 



+ I (^o + Um + uj n )\m,n)(m,n\ + ^ 4+z,m+n ^ mn |/c, Z) (m, n| 



m,n>0 



k,l,m,n>0 



The diagonal two-particle term in the last line above contributes to the equations 
determining the coefficients ujq and u n and crucially alters the renormalization 
behavior. In ||110fl , it was noted that our mean field ansatz amounts to including 
two-particle matrix elements, and it is just these terms that we have now explicitly 
displayed. Omitting them, as was done in ||110|| , leads to an uncanceled tadpole 
contribution and thus to a spurious logarithmic divergence. 

If we neglect all terms containing u mn , oj* mn , and Ukimn, which are not of inter- 
est within two-particle Tamm-Dancoff approximation, we find for the constraint 



A. 



= {m' + -T)LU 

n = (rn 2 + ^T) (uj + u n ) + — 

and for the Hamiltonian 
Ho/2L 



\ n>0 / 



(UJ + UJ n ) + 6— 

Aim 



E 

fc>0 



Airk 



(7.112) 
(7,113) 



A, 
2 



m 2 + ^T) Wo 2 + |f 



A 



m 2 m A ~ 
+ — T + —T 
2 8 



H n /2L 



1 



A 



+ 



- ( m + -T ) (u + u. 

A 
4! 



n>0 



- [m 2 + -T 



n>0 



2 
n 

Ann 



+ 



(7.114) 



27m 



+ 



(u + LO n f + I 12(W + Un)* + 20J. 



1 



47r/c ^— ' 47rA; 

fc>0 fc>0 



Ann 
+ 



A 



(7.115) 



Several remarks are in order. As 9 does not contain two-body components, the 
renormalization is correct up to the one-particle sector. (|7.112 ) and ( |7.113| ) thus 
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coincide with ( |7.32| ) and ( [7.37] ). The coefficient 6q is not even changed by includ- 
ing the two-particle contributions as can be seen by comparing with (|7.98|) . In the 
vacuum coefficient Hq, the mass renormalization (due to the tp 2 contributions) 
and the vacuum energy m 2 T/2 + \T 2 /8 (with the T 2 contribution stemming from 
the <y9 4 -term) are correct, as expected. So Ho is consistently renormalized. In the 
one-particle coefficient H n , which should be compared with ( |7.38|) , only the mass 
renormalization in the a»-sector is correct, as this is due to one-body contributions 
like uj 2 ip 2 . As anticipated, mass renormalization and vacuum energy stemming 
from the <£> 4 -term differ from the correct values by numerical factors. To get these 
correctly, one would have to perform a three-particle Tamm-Dancoff approxima- 
tion. Presumably, this would only change the coefficients of the divergent terms, 
whereas the coefficients of the finite terms would remain the same. This would 
then be analogous to the change in Hq by going from one-particle Tamm-Dancoff 
approximation ( |7.99| ) to two-particle Tamm-Dancoff approximation ( |7.114| ). 

Summarizing, we can say that, in order to obtain a consistent renormalization 
within a iV-particle light-cone Tamm-Dancoff approximation, one has to include 
contributions from (N + l)-particle matrix elements by appropriately redefining 
the coefficients in the Tamm-Dancoff ansatz. In this way, the Fock ansatz method 
(|7.28|) and the Tamm-Dancoff approximation become completely equivalent. 



7.5 Discussion 

In this chapter we have analyzed the vacuum structure of light-cone 4 +1 -theory 
by comparing different methods of solving for the constrained zero mode of the 
field operator. Within perturbation theory, the zero mode induces a second order 
mass correction which is vanishing in the infinite volume limit. We believe that 
to all orders in perturbation theory the zero mode only induces finite size effects, 
although we do not have a general proof. 

We have presented two non-perturbative methods to obtain a solution for 
the zero mode. An ansatz in terms of an increasing number of Fock operators, 
which we have truncated after the one-body term, seems to be the most economic 
procedure. With considerably more efforts, exactly the same results can be ob- 
tained within a light-cone Tamm-Dancoff approximation, if the renormalization 
procedure is properly designed. 

With either method we find a non-vanishing vacuum expectation value <p c 
of the field if the coupling A exceeds a critical value of A c ~ 40 m 2 , implying 
spontaneous breakdown of the reflection symmetry <p — > —(f). As the vacuum 
expectation value changes continuously, the associated phase transition is of sec- 



ond order, which has been rigorously established for the model at hand ||145 
The order parameter C , for instance, shows a square-root behavior as a func- 
tion of the coupling, so that the associated critical exponent is ft — 1/2. The 
critical behavior is thus of mean field type, which is only qualitatively correct, 
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as the 4>f + i model is in the universality class of the two-dimensional Ising model; 
thus, (3 should be 1/8. Because the calculations tend to become very messy, it is 
difficult to say, whether this shortcoming can be removed if one extends our ap- 
proximations to higher orders. It is, however, plausible that the critical behavior 
will change if more terms in the ansatz ( |7.28| ) are retained as these correspond 
to marginal operators ip n with n > 2. This issue, of course, deserves further 
investigations. 

Another problem we have to face is the absence of any volume dependence of 
the phase transition. We have been working in a finite spatial volume of length 
2L, the length scale L, however, which does not appear in the Fock measure, 
drops out of the equation ( [7.55 ) determining the critical coupling. On the other 
hand, there cannot be a phase transition in a finite volume due to topological 
fluctuations (kinks and anti-kinks) which have non-vanishing statistical weight for 
L < oo. We have not incorporated these fluctuations by our choice of (periodic) 
boundary conditions, so it is perhaps not too astonishing that we do not obtain 
a volume dependence. It should, in addition, be pointed out that the light-cone 
spatial volume, —L < x < L, is of a very peculiar nature. It corresponds to 
a light-like Minkowski region, the invariant length of which vanishes. This may 



give rise to features that are not present if the volume is space-like JF7 . 

A better understanding of this problem is certainly desirable. As a first step 
towards doing better, a pure continuum formulation of the problem has been 



developed recently that leads to a slightly different critical coupling ||146|| , namely 



g c = 0.69812669 . 



(7.116) 



Note that in ||146|| a different coupling normalization, g = 6g, is used so that 
g c = 4.18876012. . .. The value (|7.116|) can be straightforwardly obtained from 
our equation ( [7.55 ) if the continuum limit is performed as in Section 7.2, i.e. via 
replacing the Digamma function by its asymptotics (|7.23|) . Instead of (|7.55| ), we 
find, 

1 - g c \og(Qg c ) = , (7.117) 



which is solved by Q7.116 ). This agreement of two independent calculations pro- 
vides additional confidence in our method. We also point out that the critical 
coupling ( [7.116Q is not too different from the one obtained within sophisticated 



renormalization group treatments [Q, which, using our normalization, is given 

by g c ,RG 



0.84 ±0.11 P6 |. 



Chapter 8 

't Hooft and Schwinger Model 



In the last chapter we have seen that zero modes carry at least some of the vacuum 
structure for a scalar field theory. Unfortunately, the methods used there cannot 
be straightforwardly generalized to the case of fermionic fields. There, a non- 
trivial vacuum shows up through the appearance of a fermion condensate, which 
in QCD, for example, signals spontaneous chiral symmetry breaking, as we have 
discussed in the introduction,. The point now is, that a fermion condensate is 
bilinear in the fields, and we do not expect fermionic zero modes to play any 
particular role in the formation of the condensate. This is also supported by the 
fact that a single fermionic zero mode cannot condense as it does not have vacuum 
quantum numbers. We will therefore follow a different route in this chapter and 
try to reconstruct vacuum properties from the particle spectrum. This is where 
light-cone wave functions enter the stage. 



8.1 Prelude: Fermion Condensates 

Assume that we have a symmetry which is explicitly broken so that the associated 
current is not conserved, 

d^f{x) = A{x) . (8.1) 



With the help of the corresponding Ward identity [fL47 ] one can derive the fol- 
lowing formula for an arbitrary operator B, 



(0\SB/5a\0) = -ij d 4 x(0\TA(x)B(0)\0) = -J^ 



(0\A(0)\n){n\B(m _ (g2) 



m 2 n 



Here, 6B/6a denotes the change of B under the symmetry transformation, and 
in the last step a complete set of states |n), each of mass m n , has been inserted. 
For the case of chiral symmetry, choosing 

A = B = 2m^z7 5 V> , (8.3) 



101 



102 



CHAPTER 8. 'T HOOFT AND SCHWINGER MODEL 



one finds for a single quark flavor of mass m 



<0|fr|0) = _ mE l»|n>| 2 



In ||148|, p.49|| , this expression was used as a definition for the quark condensate in 



the 't Hooft model (see the next sections). From the above derivation, however, 
it is clear that flS.4| ) holds quite generally. Using, for instance, the PCAC relation 
for the axial vector current, j£ , 

9n35 {x) = f v ml ?r(x) , (8.5) 



where n(x) is an interpolating pion field and f w the pion decay constant, Q3.2| ) 
involves the pion two-point function and is easily evaluated with the result 

flml = -4m(0|^|0) . (8.6) 



This is the famous Gell-Mann-Oakes-Renner relation ||150|| (to lowest order in the 
quark mass m [fL51[| ), which relates the QCD parameters, m, the current quark 
mass, and the fermion condensate to the observables f n and m n . We have written 
everything in terms of bare quantities since the right-hand side does not change 
under renormalization ||152|| . Thus, m n in (513) is the physical pion mass. 



The Gell-Mann-Oakes-Renner relation can of course also be derived from ( |S.4|) 
if one replaces 

${x)i^il)(x) = — Umlirix) , (8.7) 
2m 

and assumes that, for small m^, the sum is saturated by the pion. 

In any case, we want to stress the fact that in (|8.4j) and flSlf ) above a vacuum 
quantity, the condensate, is expressed in terms of the particle spectrum. So, 
once the spectrum is known, after, say, diagonalizing the light-cone Hamiltonian 
by one of the various methods on the market, we can translate back properties 
of the spectrum into properties of the vacuum. In view of that, we suggest to 
deemphasize the role of the vacuum, which is natural to the extent that most 
of its properties are not directly observable. This is particularly true within the 
light-cone framework, where the vacuum state seems to decouple completely from 
the particle states. Similar ideas have been put forward long ago, in the context 
of chiral symmetry in the (light-cone) parton model, by Susskind et al. "In 
this framework the spontaneous symmetry breakdown must be attributed to the 



properties of the hadron's wave function and not to the vacuum" ||154|| . A related 



point of view has also been taken more recently in 

Before we pursue the program just outlined we would like to remark that the 
'master equation' ( |8.2| ) cannot be derived by strictly sticking to the light-cone 
framework. To this end note that the first term in ( |3.2|) can be written with the 
help of the charge 

Q( x °) = [ d 3 xj°(x°,x) , (8.8) 
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the generator of the symmetry, as 

(Q\6B/6a\0) = -i(0\[B,Q]\0) , (8.9) 

where the commutator is evaluated at equal time x° = 0. This expression can- 
not be directly translated into the light-cone language by replacing the ordinary 
charge Q(x°) by the light-cone or "light-like" charge, 

(W./drAxmr^, (8-10) 

and evaluating the commutator at equal light-cone time. The reason for this 
is a peculiar property of light-like charges which was briefly mentioned in the 
introduction to Chapter 7: they annihilate the vacuum, irrespective of whether 



they generate a symmetry or not |66|, |155||, which is in accordance with the 



triviality of the light-cone vacuum. Thus, the right-hand side of ( |8.9[) , evaluated 
on a null-plane, is always zero. Hence, that part of the operator 5B/5a having a 
non-vanishing vacuum expectation value cannot be obtained by an infinitesimal 
transformation generated by the light-like charge Q. For example, in the light- 
cone sigma model, the relation 

[n,Q 5 ] = -ia (8.11) 

does only hold for those modes of the field operators, 7r and a, having non- 
vanishing light-cone three-momenta, {p + ,p±) ^ 0. These non-zero modes do not 
have a vacuum expectation value. Thus, the vacuum expectation value of ( 8.11] ) 
consistently vanishes on both sides, as it should. 



The moral is that we have to assume the validity of the identity ( |3.2|) for any 
possible choice of quantization hypersurface, in particular for null planes. Ac- 
cordingly, we can use a complete set of eigenstates of the light- cone Hamiltonian 
in the last term of (|8.2[ ), as was done in | ]148| , |149j | . This amounts to solving the 
light-cone Schrodinger equation for the corresponding light-cone wave functions. 

In the following, we will test the ideas outlined above for two particular mod- 
els due to 't Hooft and Schwinger, respectively. Our presentation is a slightly 
modified version of 



8.2 State of the Art 



The number of papers on the Schwinger ||157|| and 't Hooft [ 127 , |158|| models is 
legion - for a twofold reason. On the one hand, the two models are particularly 
simple and become, for special choices of parameters, even exactly solvable. On 
the other hand, despite their simplicity, both models contain interesting physics 
analogous or similar to properties of gauge theories in higher dimensions. The 
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models therefore serve as interesting theoretical laboratories for studying phe- 
nomena like confinement and chiral symmetry breaking, to mention only the most 



prominent ones. For recent reviews on the two models we refer the reader to [159] 



for the Schwinger model and to ||160|| for the 't Hooft model and generalizations 
thereof. 

While confinement (or charge screening) is realized in the same way in both 
models, via a linearly rising Coulomb potential, the second feature, chiral sym- 
metry breaking, is not. In the 't Hooft model, 1+1 dimensional QCD in the limit 
of large Nc, chiral symmetry is 'almost' spontaneously broken | |161| , |162| | and a 



massless bound state arises in the chiral limit of vanishing quark mass [127, 158 



In the massive Schwinger model, QED in d = 1+1 [|163| , [164] , chiral symmetry is 



anomalously broken, and the contribution of the anomaly to the mass gap sur- 
vives the chiral limit yielding the free, massive boson of the massless model. This 
boson becomes interacting in the massive model and again can be viewed as a 
bound state of fermionic constituents. For the sake of brevity we will call the low- 
est bound state of both models the 'pion' (although, in the real world of d — 3+1, 
the Schwinger model boson rather corresponds to the r/-meson). In this chapter 
we will try to determine its mass and (light-cone) wave function as accurately 
as possible. We mention in passing that the Schwinger model chiral anomaly is 
closely related to the appearance of a vacuum 6*-angle parameter ||165| , |166|| which 



also affects the particle spectrum ||164||. Throughout this chapter, however, we 



implicitly assume that 9 is set to zero. 

The calculation of bound state masses and wave functions for the two models 
at hand has a long history, beginning with the exact solution of the massless 
Schwinger model [ 157| |. For non- vanishing 'electron' mass, m, the model is no 



longer exactly solvable and one has to resort to approximations. One of them 
is to assume the mass m being small and expand around the massless solution 



163], |164|| . For the 'pion' mass squared one expects an expansion of the form 

M 2 (m) = 1 + M x m + M 2 m 2 + M 3 m 3 + 0(m 4 ) , (8.12) 

where all masses are measured in units of the basic scale /z = e/y/n, the mass 
of the boson in the massless Schwinger model, which is thus represented by the 
'1' in (|3.12| ). The first order coefficient Mi was obtained analytically in ||167|| via 
bosonization, 

Mi = 2e 7 = 3.56215 . (8.13) 
with 7 = 0.577216 being Euler's constant. Shortly afterwards Bergknoff, using 



light-cone Hamiltonian techniques (see below), found a value [[129 



Mi = 2n/V3 = 3.62760 , (8.14) 



which differs from (|8.13|) by 1.8 %. One topic of this chapter will be the analysis of 
this discrepancy and an attempt to do better by refining 't Hooft 's and Bergknoff 's 
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methods. This is particularly important as the coefficient Mi is directly related 
to the vacuum structure of the Schwinger model by ||168| 



Mi 



-4tt(0|^|0) cos 9 



(8.15) 



where (0| ?/>?/> |0) is the condensate of the massless Schwinger model [|169| , |1 70jl (in 
units of /io), 

1 .. 



(o\W\o) 



2tt 



-0.28347 



(8.16) 



and 9 the vacuum 0-angle. The result ( |3.13D thus corresponds to 9 = 0. In |p9|l , 
the relation (|8.15|) (for 9 = 0), which is a 1+1 dimensional analogue of the Gell- 
Mann-Oakes-Renner formula fl3.6|) , has been used to determine the condensate 
from the 'pion' mass squared via 



(0|^|0) 



l__d_ 

Att dm 



M 2 (m) 



m=0 



—Mi . 
4tt 



(8.17) 



Any inaccuracy in the determination of Mi thus immediately affects the value of 
the condensate [[4^, p^Hf). 

Recently, the calculations of M 2 (m) have been extended to second order. 



Using functional integral techniques and mass perturbation theory, Adam ||175 
found the analytical expression 



M 2 =4 7 r 2 (0|^|0) 2 (Acos2^ + J B) 
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where A = —0.6599 and B = 1.7277 are numerical constants, given in terms of 
integral expressions containing the 'pion' propagator for m = 0. Inserting the 
values for A, B and the condensate, and setting 9 = 0, (|3.18 ) becomes 



M 2 = 3.3874 



(8.19) 



This result has been confirmed independently by Fields et al. ||176| , who derived 
the same integral expressions by summing up all relevant Feynman diagrams in 
the bosonized theory using near-light-cone coordinates. 

To first order in m, the mass-squared of the 't Hooft model 'pion' has already 



been calculated by 't Hooft ||158|| by solving the associated light-cone bound-state 
equation. He derived this equation by projecting the covariant Bethe-Salpeter 
equation onto three-dimensional hypersurfaces of equal light-cone time, x + . Soon 
afterwards, the equation was rederived using light-cone Hamiltonian techniques 



177|| . The result for the 'pion' mass-squared is 



M 2 (m) 



71 



2—m + 0(m 2 ) 
v3 



(8.20) 



1 At this point it should be mentioned that within light-cone quantization there have been 
many attempts to calculate the condensate of the massless model alternatively by solving for 
its vacuum structure |7l[ ^ [ml], [l05| ^72[ [fi)e| [l73], |l~74 |. 
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As the light-cone bound state equations of the two models at hand differ only by 
an additive contribution due to the anomaly (see |178|| and below), it is not 



too surprising that 't Hooft's and Bergknoff's results coincide. For the 't Hooft 
model, all masses are expressed in units of the basic scale /Zq = g 2 Nc/2n 0. It is 
obvious from ( |8.20| ) that M 2 vanishes in the chiral limit, m — > 0. As explained in 
161 , |162|| , this is not in contradiction with Coleman's theorem [|179|| as the 'pion' 



is not a Goldstone boson and decouples from the S-matrix. 

The condensate of the 't Hooft model has first been calculated by Zhitnitsky 
162|| using basically the identity 



{0\W\0)/N< 



c 



4tt 



M 1 



-0.28868 . 



(8.21) 



As expected, the condensate is proportional to Nq and can again be written in 
terms of Mi. The result has been confirmed analytically [ |180| , 57] and numerically 
181[ . Via (|8.21|) , the numerical value [|181|1 for the condensate leads to a numerical 



estimate for Mi, 



Mi = 3.64 ± 0.05 



(8.22) 



Though the numerical calculation did not use light-cone methods, it seems to 
favor the 't Hooft-Bergknoff value (|8.14|) , which is the standard value for the 
't Hooft model. Higher order corrections to (|8.20|) have been discussed in ||182 
without explicit calculation of the expansion coefficients. 

In recent years, both models have been serving as a testing ground for new 
techniques developed in order to solve bound state problems using light-cone 
(or, equivalently, light-front) quantization. These new methods are 'discretized 
light-cone quantization' (DLCQ) |j3], \ 



p5| , where one works in a finite volume 
leading to an equally-spaced momentum grid (cf. Chapters 4 and 7), and the 
'light-front Tamm-Dancoff approximation' | 128|| , which (drastically) truncates the 
Fock space of the theory (cf. Section 7.4), thus limiting the number of constituents 
a bound state can have. The latter approach can be viewed as a generalization 
of the techniques of 't Hooft and Bergknoff. Both methods aim at a numerical 
solution of relativistic bound state problems. DLCQ has been applied to the 
massive Schwinger model |183| and to QCD in 1+1 dimensions (for arbitrary 



Nq) ||1 78| , 184 ]. There are analogous LFTD calculations for both models, QCDi + i 



185| , |186|| and the Schwinger model ||187| , |188|| . These latter works are closer in 



spirit to ours than the DLCQ approaches. We will compare to all this recent 
work in more detail later on. 

The purpose of this chapter is to obtain the 'pion' mass squared of both the 
Schwinger and 't Hooft model to high accuracy including the third order in m. 
In addition, we want to calculate the light-cone wave function of the 'pion' with 
high precision. We will use the light-front techniques of 't Hooft and Bergknoff 



2 Some authors (including 't Hooft) use a different convention for the coupling, which amounts 
to the replacement g 2 — > 2g 2 . 
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and extensions thereof. Our starting point is a Tamm-Dancoff truncation to the 
two-particle (valence) sector. The low order calculations will be done analytically. 
To go beyond, computer algebraic and numerical methods will be applied and 
their convergence tested. The final goal is to shed some light on the advantages 
and limitations of this particular approach to bound-state equations. 

we 



The remainder of this chapter is organized as follows. In Section 8.3 



review 't Hooft's ansatz for the wave function yielding the lowest order solution 
of the bound-state equation. We compare the exact endpoint analysis with a 
variational procedure which is introduced at this point. In Section |8.4j we extend 
't Hooft's ansatz by adding more variational parameters in such a way that an 
analytic solution can still be obtained. To this end we employ computer algebraic 
methods which allow to treat up to five variational parameters. This is sufficient 
to achieve rather good convergence. In Section BJ3, these results are compared 
with the outcome of purely numerical calculations. We conclude in Section 
with some discussion of the presented as well as related work. 



8.3 't Hooft's Ansatz 



Our starting point is the bound state equation of the 't Hooft and Schwinger 
model in the two-particle sector, which, in a unified way, can be written as 



3S, 178 



M^{x) = (m 2 -l) f ix \-P CdyJ^k 



x(l - x) J (x - y) 2 

+ a dx(f)(x) . (8.23) 
Jo 

This expression corresponds to the first line of Q6.18 ) with only the matrix ele- 



ment W 2 = (qq\W\qq) retained. We will refer to ( |3.23|) as the 't Hooft-Bergknoff 



equation in what follows. 4>(x) denotes the valence part of the 'pion' wave func- 
tion, x and y the momentum fraction of one of the two fermions in the meson. 



The symbol V indicates that the integral is defined as a principal value ||127 



It regularizes the Coulomb singularity l/(x — y) in the matrix element 



2- 



The parameter a measures the strength of the anomaly. In the 't Hooft 
model, a = (no anomaly), and in the Schwinger model a = 1 (representing the 
usual chiral anomaly). The scale parameters, as already mentioned, are given by 
Ho = g 2 Nc/2ir and /i 2 , = e 2 /7r, respectively. M denotes the mass of the lowest 
lying bound state (the 'pion'). Our objective is to obtain (approximate, but 
accurate) solutions for M and 0. 

Upon multiplying ( |3.23|) with 4>(x) and integrating over x we obtain for the 



eigenvalue 

M 2 (m\ = (m 2 - r, 

-^0 -^0 -^0 



M 2 (m) = (m 2 -l)^-^ + a$ , (8.24) 
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where we have defined the integrals 



J = / dx(tf(x) , (8.25) 



o 



1 A.2i 



/, = / dx-j-^—r , (8.26) 
x 1 — X) 





I 2 = dxdy J — w , (8.27) 

dx(f>{x) . (8.28) 

o 

Jo is the norm of the wave function, Ii and I 2 are matrix elements of the mass 
and interaction term in the light-cone Hamiltonian ||129| , |187|| in the state \<f>). 
The integral I3 is the wave function at the origin. An independent formula for 
the 'pion' mass-squared is obtained by integrating (|8.23|) over x. This results in 
the simple expression 



M 2 (m) = m 2 — + a , (8.29) 
h 



with the additional integral, 



dx , . (8.30) 

x(l — x) 



Of course, for the exact wave functions, the right-hand sides of ( |S.24 ) and ( |3.29 ) 



have to coincide. As the wave functions cannot be obtained exactly, we will later 
use the agreement between both values for the mass-squared as a measure for the 
accuracy of our wave functions. To determine the latter we will use a particular 
class of variational ansatze. 



In his original work on the subject, ||127| , |158|| , 't Hooft used the following 
ansatz for the wave function 

<p(x) = x 13 (I - xf . (8.31) 

This ansatz is symmetric in x «-> 1—x (charge conjugation odd), and (3 is supposed 
to lie between zero and one so that the endpoint behavior is non-analytic. As a 
non-trivial boundary condition one has the exact solution of the massless case, 

M 2 = a, and <j>{x) = 1 , i.e. (3 = 0. (8.32) 

The main effect of having a non- vanishing fermion mass is the vanishing of the 
wave functions at the endpoints implying a non- zero (3. This suggests the follow- 
ing series expansion for (3, 

(3{m) = p 1 m + p 2 m 2 + p 3 m 3 + 0{m A ) , (8.33) 

and for the 'pion' mass squared, 

M 2 = a + M x m + M 2 m 2 + M 3 m 3 + 0(m 4 ) , (8.34) 

in accordance with ( |S.12j ). 
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8.3.1 Exact Endpoint Behavior 



The exponent (3 in ( P-31|) can actually be determined exactly by studying the 
small-x behavior of the bound state equation Q8.23Q . To this end we evaluate the 
principal value integral for x — > 0, 



Jo (x - y) 2 



x 



V 



dz- 



y/3 



—x 



(3-1 



7T/3 cot 71/3 + 0(x 



+ 0(x) 



(8.35) 



Plugging this into (|8.23| ) and demanding the coefficient of x 13 1 to vanish yields 



the transcendental equation ||127 



m 2 - 1 + up cot up = 



(8.36) 



which is independent of the anomaly a. Using this expression we can determine 
P either numerically for arbitrary m or analytically for small m by expanding 



, \/3 A 1 2 
p = — m 1 m 

H 71 \ 10 



0{m A 



(8.37) 



Note that the second order coefficient, P2, is vanishing. Furthermore, for the 
exact P one has P1/P3 = —1/10, which we will use as a check for our numerical 
results later on. 

Our task is now to determine the coefficients Mj, % = 1,2,3 in flB.34j ). The 
ansatz (|8.31| ) leads to the following results for the for the 'pion' mass squared 



M 2 (m,p) 



[m 



+ 



+ P 



B 2 (P,P) 
B(2P,2P) 



1 + a 



(2/3 +1)^ 



(8.38) 



In the above, B(zi,z 2 ) denotes the Beta function. The relevant formulae for 
double integrals like I2 (|8.27|) can be found in ||189|| and ||190|| , Appendix C. For 
a = 0, (|8.38|) has also been obtained in [|186|1 . Let us emphasize that this result, 
which expresses the 'pion' mass squared as a function of the (exactly known) 
endpoint exponent P, is only approximate as 't Hooft's ansatz ( |S.31|) for the wave 
function does not represent an exact solution of the bound-state equation. It is 
only the endpoint behavior that is known exactly; in the intermediate- a; region 
't Hooft's ansatz is only an approximation that presumably becomes rather bad 
for large masses (non-relativistic limit) where the wave function is strongly peaked 
at x = 1/2. The accuracy of the result will be discussed extensively later on. 
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In order to find M 2 to order m 3 we need to expand ( |8.38| ) to order /3 3 , as (3 
itself is of order m. The result is 



M 2 (m,/?) = — + 4m 2 + a + —f3 



+ 4 

r3 



tt 2 /3 - 3C(3) + a(n 2 /12 - 1) 



2 



+ 



-vr 4 -48C(3) + 4a 4-3C(3) p 



+ 0(/3 4 ) 
Inserting (3 from ( 8.371) one finds 



(8.39) 



M 2 {m) 



+ 



+ 



7T 

a + 2—=m 
4(2-»C(3)) + «(l-H) 



3v/3 



7T" 



-7T 

5 



48C(3))+4a(4-3C(3) 



+ 0(m 4 ) . 



(8.40) 



Let us give the explicit results for the 't Hooft and Schwinger model, i.e. for a 
and a =1, respectively. For a = 0, we have 



Mi 
M 2 
M 3 



2— = 3.627599 
\/3 

9 



3.615422 



4(2 -C(3)- 

!^( r ! 7r 4_48C(3) > ) = 0.125139 



and for a 



Mi 
M 2 
M 3 



7T 3 \5 



2^ = 3.627599 

<2-c(3)^; 

3^ 



1 - 4 = 3.399568 

7T 2 



7T' ; 



-7T 



48C(3))+4(4-3C(3) 



0.389137 . 



(8.41) 
(8.42) 
(8.43) 

(8.44) 
(8.45) 

(8.46) 



We note that the anomaly does not contribute to the first order term which 
therefore is the same in both the Schwinger and 't Hooft model. This confirms 
the observation made in Section 8.2, the coincidence of ( ft.14 ) and ( ft- 20] ). As 
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Mi is proportional to the condensate, cf. ( |3.17| ) and ( |B.21| ), the latter is also 
independent of the anomaly 



For the Schwinger model (a = 1), the second order result can be compared 
with the analytical calculation of Adam, ( |8.19| ). Astonishingly, the relative dif- 
ference is smaller than for the first order, namely 0.36 %. This is accidental, as 
we shall see below. 

As a cross check, we determine the 'pion' mass squared using the alternative 
formula ( |8.29| ) . Plugging in the ansatz flS.31| ) , we find the fairly simple expression 



2m 2 

M 2 (m, (3) = Am 2 + a + — — . (8.47) 

P 



Inserting the expansion ( |S.37| ) this becomes 



M 2 = a + M x m + M 2 m 2 + M 3 m 3 + 0(m 4 ) , (8.48) 

with the coefficients Mj explicitly given by 

Mi = 2^/^ = 3.627599, (8.49) 
M 2 = 4, (8.50) 
M 3 = Mi/10 = 0.03627599 . (8.51) 

In ( |8.48j ) the anomaly contributes only in zeroth order (as (3 is independent of 
a), so that the Mj are the same for a = and a = 1. Comparing with Q8.41 ) 



and ( |8.44| ), we see that the first order coefficients coincide for the two alternative 



mass formulae. In addition, one finds that M 2 roughly coincides with the M 2 of 
both the 't Hooft and Schwinger model to within 10-20 %. The value for M 3 is 
smaller by approximately an order of magnitude. 

We conclude that 't Hooft's ansatz works well to lowest non-trivial order in 
m but needs improvement if one wants to go further. As a preparative step for 
such a development we now introduce a variational method that can easily be 
extended to accurately calculate higher orders in m. 



8.3.2 Variational Approach 

At variance with the above we are now going to regard (3 as a variational param- 
eter to be determined by minimizing the function M 2 (m,(3) of (|8.39| ). To this 
end we insert the expansion (|8.33|) into (|8.39|) and obtain 



M 2 (m) 



a + 



n 2 1 



III 



+ 



4 + &( 



7T 

y 



i 
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+ 4 (t 



3C(3) + a 



7T 

12 



+ 
+ 



I lift 2 
3 



m 



8/?i/3 2 (y - 3C(3)) + ^(-tt 4 - 48C(3)) 



+ 0z 



71 



Pi 

0f ' ^V3 A 2 



+ a(2A/? 2 (y-4) +4/3 1 3 (4-3C(3) 
+ 0(m 4 ) . 



(8.52) 



Note that to order m only the leading coefficient 0\ contributes. We will fur- 
thermore shortly see that the dependence of M 2 on 2 and of M 3 on 3 is only 
apparent. 

Solving the minimization equation, dM 2 /d/3 = 0, for the coefficients 0i, leads 

to 



0i 

02 
03 



V3/tt = 0.55133 , 

0.11690 + 0.065612a , 

0.0049077 - 0.050811a + 0.019521a 5 



(8.53) 
(8.54) 
(8.55) 



Comparing with ( |8.37| ) we note that the coefficient 0\ is exact! Plugging it 
into (|8.52j ) we verify the statement above that M 2 is independent of 02 and M 3 
independent of 0^. Therefore, the expressions ( |8.52| ) and ( |8.40| ) coincide up to 
and including order m 2 . Mi and M 2 are thus the same, irrespective of whether 
one uses the exact endpoint exponent of ( |S . 3 7| ) or its variational estimate. The 
estimates ( |3.53^|S35| ) for the coefficients 0i differ in at least three respects from the 
exact values of ( ^.37] ): (i) 2 and 3 depend on a, (ii) 2 ^ 0, (hi) 3 ^ —0i/lQ. 
However, all these shortcomings affect at most the third order coefficient M3, and 
thus become negligible for small m. Only if one wants to have reliable numbers 
for M3, (which we will be going to produce), these effects have to be taken into 
account. The present values of 0i lead to 



M 3 (a = 0) = 0.043597 



M,(a = r 



0.190372 



.56) 



which should be compared with (|8.43|) and ( |3.46|) , respectively. The differences 
are large so that the agreement is not particularly good. These discrepancies, 
however, are not disturbing at this point, as we are calculating a third order effect 
with just one variational parameter. We will do better later on by enlarging the 
number of these parameters until we see satisfactory convergence of the results. 

The variational estimate for can also be plugged into the alternative mass 
formula, M 2 (m,0), (|8.47|) , which by itself does not constitute a variational prob- 
lem. We do not give the analytical results for the Mi here but simply refer 
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Figure 8.1: Comparison of the two alternative mass formulae, ( 8.3$) for M 2 (f3) 
and ( |8.47p for M 2 {(3), and the second order bosonization results \173[ , \170^ for 
m = 0.1. The vertical line marks the minimum of the curve M 2 (/3) yielding 
the variational estimate for f3 using 't Hooft's ansatz. For this value of (3 one 
finds M 2 (m = 0.1) = 1.3969 and M 2 {m = 0.1) = 1.3913. The bosonization 
result (short-dashed horizontal line) is M^ os (m = 0.1) = 1.390 ± 1 • 10~ 3 , where 
the error is basically an estimate of the unknown third order coefficient. Thus, 
within 't Hooft's ansatz, the alternative mass formula yields a somewhat better 
result at m = 0.1. In the next sections we will refine our methods so that the 
variational estimates for M 2 and M 2 will converge towards each other. 




to Fig. |3.1| for a qualitative comparison of the two alternative formulae, and to 
Section Tables |873| and pT7| , for the actual numbers. 

As stated repeatedly, 't Hooft's ansatz (|8.31| ) used so far has to be extended 
if one wants to accurately determine the second and third order coefficients of 
M 2 . This is our next issue. 



8.4 Extension of 't Hooft's Ansatz 



For the numerical calculations of ||186| , |190| , |188|| the following set of trial functions 
has been used: 

N 

tf>(x) = ^u k [x(l-x)f +k . (8.57) 

k=0 
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Obviously, the term with k = and normalization uq = 1 corresponds to 
't Hooft's original ansatz ( |3.31 ). The coefficients Uk are treated as additional 



variational parameters, so that, according to the variational principle, ( p. 57 ) 



must yield a better result than flS.31|) . The question is, how big the improve- 



ment will be. To this end we will calculate the coefficients Mj by adding more 
and more basis functions to 't Hooft's original ansatz ( |8.31|) , thus enlarging our 



space of variational parameters. We will follow the two different approaches men- 
tioned above, namely, (i) use the exact (3 from (|8.37|), or, (ii) treat (3 as one of 



the variational parameters. For both approaches we will compare the calculated 
coefficients Mj with the Mj obtained from the alternative mass formula (|S.29 ). 



We will continue adding basis functions until we see our results converge. The 
maximum number N of basis functions in this section will be five, i.e. 't Hooft's 
original wave function plus four corrections. 

A particular benefit of the ansatz ( |8.57| ) is the fact that the integrals I± to 
I4 can still be evaluated analytically though the formulae become rather lengthy. 
For this reason we will make heavy use of the program package MAPLE in what 
follows. As a result we have analytical expressions for all the quantities we calcu- 
late. As these expressions cover pages and pages without being very instructive, 
we do not display them but rather the evaluated numbers. In this sense, our 
treatment might be called 'semi-analytical'. A major advantage, however, of us- 
ing computer algebraic manipulations is the high calculational accuracy which is 
only limited by the maximum number of digits the machine can handle. 

8.4.1 Purely Variational Approach 

In this subsection, (3 will always be treated as a variational parameter and thus 
be obtained by minimizing the 'pion' mass squared, M 2 . In the ansatz ( |3.57| ) we 



use up to four additional basis functions, so that our maximum N is four. To 
avoid an inflation of indices we rename the coefficients k — 1, . . . , 4 by a, b, c 
and d, respectively. Their expansion coefficients, defined via 

a = aim + a 2 m 2 , 

b = b\m + &2^ 2 j 

c = C\m + C2m 2 , 

d = dim + d^rn 1 , (8.58) 

are determined (recursively) by minimizing M 2 with respect to them. Let us 
start with the 't Hooft model (a = 0). 

't Hooft Model 



In order to save space we do not display all the values for the coefficients in ( |S.58| ). 
The best values obtained are given in the last section when we discuss the quality 
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Table 8.1: Expansion coefficients of the end point exponent (3 for the 't Hooft 
model obtained by successively varying with respect to (3 (first line) , P and a 
(second line), etc. The hrst column, (3\, coincides with the result from the exact 
endpoint analysis, which in addition yields p 2 = 0, (3$ = — Pi/ 10. 

Ansatz Pi = ^/3/tt p 2 P3 

't Hooft 0.55132889 0.11689763 0.00490773 

a 0.55132889 0.00976951 -0.04010634 

b 0.55132889 0.00256081 -0.04889172 

c 0.55132889 0.00102086 -0.05209341 

d 0.55132889 0.00050738 -0.05343097 



of the wave function. Here, we rather display the results for the coefficients Pi of 
P (see Table |8~T| ), as these can be compared with the exact values of (|8.37f) . 

As an important finding we note that Pi remains unchanged at the exact 
value v^3/?r (a 'variational invariant'). P2 tends to zero, as it should. The non- 
vanishing of P2 only affects the coefficient M 3 , as Mi and M 2 do not depend on 
P 2 . The convergence of P3 towards — Pi/10 seems somewhat slow; but as all Mi, 
i = 1, 2, 3, are independent of P3 this has no observable effect. 

In Table |8.2| we list the expansion coefficients of the mass squared, M 2 . One 
notes that the convergence of the results for M 2 and M 3 is rather good. For M 2 we 
finally have a relative accuracy of 8- 10~ 7 , and for M3 of 4- 10 -5 . Furthermore, the 
coefficients are getting smaller if one adds more basis functions, in accordance 
with the variational principle. As the coefficient Mi is entirely determined by 
the 'variational invariant' Pi it remains unchanged at 27r/y / 3. This is another 
important result: Mi stays fixed at the standard 't Hooft-Bergknoff value (|8.14|) . 



If we use the alternative mass formula (|8.29|) and evaluate it using the wave 
function calculated above we find the values for M 2 listed in Table IO. 



We note that the values of the Mj are somewhat more sensitive to the values of 
the expansion coefficients Pi, as M 2 does depend on p 2 , and M 3 on p 2 and P%. To 
check the quality of our wave functions one should compare the results for Mj and 
Mi which are listed in Tables |87| and |3.3| . Again, it is gratifying to note that we 
are improving our results step by step in the variational procedure. The values 
of Mi converge towards those of Mj. The relative accuracy is approximately 
10~ 3 . This error is entirely due to the difference between the calculated and 
the real wave function. It is clear that the accuracy in the single eigenvalue M 2 
(Table |S.2| ) is higher than the one for the wave function, where in principle an 
infinite number of points has to be calculated. This represents an indirect proof 
that we are indeed improving our wave functions, and not just the eigenvalues. 
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Table 8.2: Expansion coefficients of M 2 for the 't Hooft model obtained by 
successively enlarging the space of variational parameters. M x is the standard 
't Hooft-Bergknoff result. Note the good convergence towards the bottom of the 
table. 



Ansatz Mi = 2n/V5 M 2 M 3 

't Hooft 3.62759873 3.61542218 0.043597197 

a 3.62759873 3.58136872 0.061736701 

b 3.62759873 3.58107780 0.061805257 

c 3.62759873 3.58105821 0.061795547 

d 3.62759873 3.58105532 0.061793082 



Table 8.3: Expansion coefficients of the alternative mass squared M 2 for the 
't Hooft model obtained by successively enlarging the space of variational param- 
eters. Again, the fixed value for Mi is the standard 't Hooft-Bergknoff result. 



Ansatz M 1 = 2tx/^> M 2 M 3 

't Hooft 3.62759873 3.23084437 0.130791594 

a 3.62759873 3.54922830 0.066592425 

b 3.62759873 3.57265307 0.059652468 

c 3.62759873 3.57769969 0.060282959 

d 3.62759873 3.57938609 0.060854772 
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Table 8.4: Expansion coefficients of the endpoint exponent (3 for the Schwinger 
model obtained by successively enlarging the space of variational parameters. 
The first column, f}\, coincides with the result from the exact endpoint analysis, 
which in addition yields fa = 0, P3 = — Pi/ 10. 



Ansatz Pi P2 Ps 

't Hooft 0.55132889 0.18250945 -0.026382222 

a 0.55132889 0.01177681 -0.040372413 

b 0.55132889 0.00317437 -0.048081756 

c 0.55132889 0.00127507 -0.051602713 

d 0.55132889 0.00063609 -0.053130398 



Table 8.5: Expansion coefficients of M 2 for the Schwinger model obtained by 
successively enlarging the space of variational parameters. Again, the fixed value 
for Mi is the standard 't Hooft- Bergknoff result. Note the good convergence 
towards the bottom of the table. 



Ansatz M x = 2n/V3 M 2 M 3 

't Hooft 3.62759873 3.39956798 0.19037224 

a 3.62759873 3.30906326 0.34776772 

b 3.62759873 3.30864244 0.34820193 

c 3.62759873 3.30861240 0.34820513 

d 3.62759873 3.30860791 0.34820389 



Schwinger Model 

For the Schwinger model (a = 1), we perform exactly the same calculations. The 
expansion coefficients of P are listed in Table |S.4| . Although the coefficients P2, 
Ps are somewhat different from those in Table |8J] (which they must be as they 
lose their dependence on a only in the strict limit of exact evaluation) we note 
the same tendency: P2 converges to zero and P3 towards — Pi/ 10. 

The best values for the 'pion' mass squared are again provided by the varia- 



tional results listed in Table |8J| The numerical accuracy is practically the same 
as for the analogous Table |8.2| . The comparison with the alternatively calculated 
coefficients M 2 is given in Table |8]6|. Again, everything is completely analogous 
to the case of the 't Hooft model (a = 0). 
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Table 8.6: Expansion coefficients of the alternative mass squared M 2 for the 
Schwinger model obtained by successively enlarging the space of variational pa- 
rameters. Again, the fixed value for Mi is the standard 't Hooft-Bergknoff result. 



Ansatz M x = 2ir/y/3 M 2 M 3 

't Hooft 3.62759873 2.79913596 0.57111672 

a 3.62759873 3.27031909 0.36868190 

b 3.62759873 3.29819920 0.34869835 

c 3.62759873 3.30441759 0.34753576 

d 3.62759873 3.30651525 0.34762562 



8.4.2 Variational Approach using the exact (3 

In this subsection we calculate M 2 and M 2 for both values of a using the exact 
value QS.37D for (3. Thus only a, b, c, d are treated as variational parameters. 



This is the procedure employed numerically in ||188|| . From the discussion of the 
preceding subsection, in particular the values for the $ displayed in Table |8.1| , 
which differ minimally from the exact values, we expect that the results will be 
very close to those from the purely variational approach. 

't Hooft Model 

This is exactly what happens as can be seen from Tables |S.7| and 1573. Comparing 



Table |S.7| with Table |872| one finds that the coefficients M\ and M 2 of both tables 
coincide as these only depend on f3\ which is the same in both approaches. Only 
for M 3 there are slight differences, due to the dependence on (3 2 - F° r the M iy 
i = 2, 3, the discrepancies are somewhat bigger as these coefficients do depend 
on fi 2 an d 03 (Table |T3] vs. Table |3.8| ). Still the consistency is quite obvious. 



Upon comparing the second columns of Tables |U and |873| , respectively Ta- 
bles [T7] and 57S , one notes a funny coincidence. In the first case M 2 is slightly 



bigger than M 2 , and vice versa in the second case. If one denotes the purely 
variational results with a superscript V, and the results obtained with the exact 
with a superscript 'e', one finds that M 2 , which is the same in both approaches, 
is given by the arithmetic mean 



M 2 = - (M 2 e + M|J . (8.59) 

We have checked this analytically for 't Hooft's ansatz. For the higher orders this 
is difficult to do but the numerical evidence is beyond doubt. 

Comparing the third columns of Tables |8.2j and |8.7| , one finds that Mg < M|. 
We thus see a slight tendency that the results of the purely variational procedure 
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Table 8.7: Expansion coefficients of M 2 ('t Hooft model) obtained by using the 
exact endpoint exponent (3 and successively enlarging the space of variational 
parameters.. Again, the fixed value for M x is the standard 't Hooft-Bergknoff 
result. 



Ansatz M 1 = 2ir/V3 M 2 M 3 

't Hooft 3.62759873 3.61542219 0.12513878 

a 3.62759873 3.58136873 0.06230622 

b 3.62759873 3.58107781 0.06184441 

c 3.62759873 3.58105821 0.06180178 

d 3.62759873 3.58105533 0.06179462 



Table 8.8: Expansion coefficients of the alternatively defined mass squared M 2 
('t Hooft model) obtained by using the exact endpoint exponent f3 and succes- 
sively enlarging the space of variational parameters. Again, the fixed value for 
Mi is the standard 't Hooft-Bergknoff result. 



Ansatz M x = 2tt/V3 M 2 M 3 

't Hooft 3.62759873 40 0.36275987 

a 3.62759873 3.61350915 0.08659575 

b 3.62759873 3.58950254 0.07106398 

c 3.62759873 3.58441672 0.06603497 

d 3.62759873 3.58272458 0.06406179 
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Table 8.9: Expansion coefficients of M 2 (Schwinger model) obtained by using 
the exact endpoint exponent f3 and successively enlarging the space of variational 
parameters. Again, the fixed value for Mi is the standard 't Hooft-Bergknoff 
result. 



Ansatz M x = 2ir/y/3 M 2 M 3 

't Hooft 3.62759873 3.39956798 0.38913656 

a 3.62759873 3.30906326 0.34859533 

b 3.62759873 3.30864244 0.34826210 

c 3.62759873 3.30861239 0.34821485 

d 3.62759873 3.30860791 0.34820630 



Table 8.10: Expansion coefficients of the alternatively defined mass squared 
M 2 (Schwinger model) obtained by using the exact endpoint exponent (3 and 
successively enlarging the space of variational parameters. Again, the fixed value 
for Mi is the standard 't Hooft-Bergknoff result. 



Ansatz 

't Hooft 

a 

b 

c 

d 



Mi = 2tt/\/3 

3.627598730 

3.627598730 

3.627598730 

3.627598730 

3.627598730 



M 2 

4.0 

3.347807427 
3.319085690 
3.312807215 
3.310700593 



M3 

0.3627598730 
0.3637551073 
0.3566261465 
0.3523062162 
0.3504584408 



are better than those obtained using the exact (3. The same observation has been 
made by Mo and Perry [187], in particular for large values of the fermion mass m 
(where the third order coefficient M 3 becomes important). We interpret this fact 
as a hint that for larger m the behavior of the wave function in the intermediate 
region becomes more relevant compared to the endpoint behavior. 



Schwinger Model 



For the Schwinger model, the analogous results are listed in Tables |3.9| and |8. 10 



Exactly the same remarks as above apply including the size of the errors and the 
relation between M2 and Mf, M%. In the next section we will verify our results 
by purely numerical methods. 
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Table 8.11: Numerical results using lowest order (two-particle) light- front 
Tamm-Dancoff approximation for M 2 — 1 as a function of the fermion mass m, 
for a = 1 (Schwinger model). 



Ansatz m = 0.0001 


m = 0.0005 


m = 0.001 


m = 0.005 


m = 0.01 


m = 0.05 


m = 0.1 


't Hooft 0.000362794 


0.00181465 


0.00363100 


0.0182230 


0.0366163 


0.189927 


0.397130 


a 0.000362858 


0.00181465 


0.00363086 


0.0182207 


0.0366072 


0.189695 


0.396181 


b 


0.00181461 


0.00363100 


0.0182210 


0.0366071 


0.189694 


0.396176 


c — 


0.00181459 






0.0366071 


0.189694 


0.396176 


d 








0.0366073 


0.189694 


0.396176 



8.5 Comparison With Numerical Results 



For the numerical calculations the ansatz (|8.57| ) has been used. The value for (3 



has been determined numerically from flS.36|) . Again, up to four additional basis 
functions have been included (for larger masses even five). Let us begin with 



the Schwinger model (a = 1). Here we can use the code developed in [|188|| . In 



Table 18.111 we list M 2 — 1 as a function of the fermion mass m, calculated within 



two-particle light-front Tamm-Dancoff approximation. The notations a,b,c,d are 
as in the preceding section. 

One main difference in comparison with the computer algebraic treatment of 



Section |8.4j is the numerical inaccuracy for small m. This is a general disease 



of numerical treatments, and also shared e.g. by the lattice [168, 169|| or DLCQ 



approach ||183 |. In these approaches, however, the numerical errors are typically 



much larger than ours (see the next section for an explicit comparison). Due to 
the small-m instability, our numerical calculation does not converge within an 
arbitrary number of digits. As soon as a calculated value for M 2 — 1 becomes 
bigger as the preceding one (a numerical violation of the variational principle), 
we terminate the procedure and pick the smaller value as our final result. The 
difference between these last two values can be used as an estimate of the nu- 
merical inaccuracy. The errors will be further discussed below and in the next 
section when we compare our results with related work (see Tables |S.15| , |S.16| ). 

It should also be pointed out that the numerical analysis is conceptually very 
different from the computer algebraic treatment. As one cannot perform a Taylor 
expansion numerically, the expansion coefficients of M 2 have to be obtained by 
fitting polynomials to M 2 (m). Clearly, this is an additional source of errors, and 
one expects the results to be less accurate than those of the preceding sections. 
A cubic fit to the optimum values for M 2 — 1 in Table fj.llj yields the following 
expansion coefficients of M 2 , 

Mi = 3.62609 , (8.60) 
M 2 = 3.33607, (8.61) 
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M 3 = 0.22396, (8.62) 



which should be compared with the last line of Table |3.5| . To estimate the accu- 
racy of these values we also show the results from a quartic fit, 

M x = 3.62755 , (8.63) 
M 2 = 3.31029 , (8.64) 
M 3 = 0.32984. (8.65) 



Comparing ( p.60H8.62|) and (|8.63| - |8.65|) one finds that the stability of the fits is 



not too impressive in view of the accuracy we would like to achieve. In particular 
the third order coefficient is numerically difficult to determine. We estimate the 
relative accuracy as being 10~ 3 for Mi, 10~ 2 for M 2 and 10 _1 for M 3 . This is also 
confirmed by comparing with Table |8.5| . 

For the Schwinger model it is possible to check the influence of higher particle 
sectors on the 'pion' mass squared by using the machinery developed in [|188|| for 



the wave functions of the higher Fock components. These are the amplitudes of 
finding not only two, but four, six, ... (anti-)fermions in the 'pion'. In Table |8.12| , 
we list the best values for M 2 — 1 including up to six-body wave functions. We 
also show the 2-, 4-, and 6-body content of the total wave function. It is known 



that the 'pion' is entirely 2-particle in the chiral limit ||187|| . For small mass, one 
therefore expects only small contributions from the higher Fock sectors. This 
is confirmed by the numerical results. Astonishingly, this feature persists up to 
values of at least m = 0.5 for the fermion mass. This fact is in agreement with 
the observation of Mo and Perry that the four-particle component of the wave 



function is less than 0.4 % for all values of the fermion mass [T%7\. A similar 



result has been found in the DLCQ calculations of [ |183|| . In addition we note 
that there seems to be some interesting kind of hierarchy between the relative 
strengths of the contributions from different particle sectors. If we denote the 2k- 
particle amplitude in the wave function by f 2 k, we find that |/ 2 | 2 3> I/4I 2 2> |/e| 2 5 
the individual proportions being several orders of magnitude (see Table |B.12| ). 
Here, we explicitly see the magic of light-front field theory at work: high Fock 
components in bound states tend to be largely suppressed, at variance with the 
situation encountered in field theory quantized the standard way. 

From the point of view of the variational principle, the results shown in Ta- 
ble 8.12 are a bit better (i.e. smaller) than those of Table 8.11| (apart from the 



value for m = 0.0005), but the improvement is rather small. A cubic fit yields 

Mi = 3.62667, (8.66) 
M 2 = 3.23696 , (8.67) 
M 3 = 0.36235 , (8.68) 



and a quartic fit 



Mi = 3.62747, (8.69) 
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Table 8.12: Numerical Results for M 2 — 1 as a function of the fermion mass m, 
for a = 1 (Schwinger model) including 2-, 4- and 6-body wave functions. The 
contributions of the different Fock sectors to the total wave function squared are 
given in percent. 



m = 0.0001 


m = 0.0005 


m = 0.001 


m = 0.005 


m= 0.01 


m = 0.05 


m = 0.1 


M 2 - 1 0.000362793 


0.00181481 


0.00363055 


0.0182179 


0.0365968 


0.189468 


0.395400 


% 2-body 100.00 


100.00 


100.00 


99.999990 


99.999962 


99.999213 


99.997514 


% 4-body 






0.000010 


0.000038 


0.000782 


0.002474 


% 6-body 










0.000005 


0.000013 



Table 8.13: M 2 as a function of the fermion mass m, for a = ('t Hooft model), 
obtained with Code I. 



m 


0.0001 


0.001 


0.01 


0.1 


M 2 


0.000362795 


0.00363109 


0.0366342 


0.398634 






M 2 = 3.20864, 
M 3 = 0.60365 . 




(8.70) 
(8.71) 



Compared with ( fj.60fj8.62] ) and (|8.63| - |8.65|) we do not find any absolute improve- 
ment in Mi that could be distinguished from the numerical inaccuracy. For the 
coefficient M 2 , which in bosonization schemes is 3.3874, we even get the wrong 
tendency: it becomes smaller upon including higher Fock states. The inaccuracy 
for M 3 is so large that this coefficient is only determined in its order of magnitude. 

Altogether, we arrive at the very important conclusion that the inclusion of 
higher particle sectors in the light-front bound state equation of the Schwinger 
model does not diminish the discrepancy between the results obtained via the 
't Hooft-Bergknoff method and those from bosonization techniques. 

Let us move to the 't Hooft model. If we put a = in the Schwinger model 
code above (Code I) we find the 'pion' mass squared for the 't Hooft model. The 
best values within the 2-particle sector are listed in Table 3.13 . Note that higher 
particle sectors are strictly suppressed in the large- Nc limit. Thus, there is no 
point in calculating these, unless one is interested in 1/Nq corrections to the 
't Hooft model results |178|, HM HI • 



In view of the restricted number of data points only a quadratic fit makes 
sense which yields 



Mi = 3.62754 
M 2 = 3.58796 



(8.72) 
(8.73) 
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Table 8.14: M 2 as a function of the fermion mass m, for a = ('t Hooft model), 
obtained with Code II. 



m 


0.0001 


0.001 


0.01 


0.1 


M 2 


0.000362795 


0.00363118 


0.0366341 


0.398634 



which is consistent with the results of Section |S.4| (see Table |8.2| ). A second code 
(Code II), which was independently developed for the 't Hooft model ||185|| yields 
the results displayed in Table |8.14| . From a quadratic fit we obtain 

Mi = 3.62754, (8.74) 
M 2 = 3.58806. (8.75) 



Comparing with ( |8.72j , |3.73| ) we see that both codes yield the same results within 



the numerical accuracy. This is reassuring, since, as already stated, the codes 
were developed independently from each other. 

A cubic fit to a total of 30 data points produced with Code II yields 

Mi = 3.62758, (8.76) 
M 2 = 3.58260, (8.77) 
M 3 = 0.06450, (8.78) 



in fair agreement with the variationally obtained values of Table ft.2| . To check 
the stability of this fit we can compare with the extension to fourth order, 

Mi = 3.62756, (8.79) 

M 2 = 3.58314, (8.80) 

M 3 = 0.06280, (8.81) 

M 4 = 0.001232. (8.82) 

This leads to an estimate of roughly 10~ 6 , 10~ 4 and 10~ 2 for the (absolute) numer- 
ical error in the first, second and third order coefficient, respectively. Furthermore 
it is gratifying to note that the fourth order coefficient, M 4 , is numerically small. 

As is obvious from the discussion above, the numerical errors for the 't Hooft 
model are much smaller than those for the Schwinger model. The basic reason for 
this is the Schwinger model anomaly. Note that in the bound state equation the 
anomaly factor a multiplies an integral over the wave function. To evaluate this, 
one needs the wave function as a whole. In the 't Hooft model, on the other hand, 
this term is absent and the eigenvalues are dominated entirely by the endpoint 
behavior of the wave function which is known exactly. This makes the numerical 
errors much smaller (see also Tables |8.15| , |8.16| ) . 
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8.6 Discussion 

In the preceding sections we have calculated the expansion coefficients Mj in the 
series 

M 2 = a + M x m + M 2 m 2 + M 3 m 3 (8.83) 

for the 'pion' mass squared of both the 't Hooft (a = 0) and the Schwinger model 
(a = 1). In order for the expansion flS.83 ) to make sense we have considered only 



small masses m <C 1 i.e. m <C /xq in the original units. For the 't Hooft model 
we have = g 2 Nc/2ir with Nc — > oo, g 2 Nc fixed. Thus, we are working in 
the weak coupling phase where the limit Nc — > oo (g — > 0), is taken before the 
limit m — ► 0, or, equivalently, such that one always has m ^> g ~ 1 fy/~N~c — > 
For the Schwinger model, /J 2 , = e 2 /ir, so that small fermion mass 



corresponds to strong coupling. 

We have used analytical, computer algebraic and numerical methods of dif- 
ferent accuracy. Nonetheless, the overall picture is intrinsically consistent. Our 
main findings are the following: 

(i) The first order coefficient, Mi, in the expansion of the 'pion' mass squared 
is independent of the anomaly a, i.e. the same in both the 't Hooft and Schwinger 
model. This confirms the results of 't Hooft ||127| , |158|| and Bergknoff | 129|| . 



(ii) The 't Hooft-Bergknoff value, Mi = 2n/\^3, is a 'variational invariant': it 
does not get altered by extending the space of variational parameters. This has 
been checked analytically and numerically. 

(iii) In the Schwinger model, the 't Hooft-Bergknoff value does not change 
upon inclusion of higher particle sectors. Only the second and third order coeffi- 
cients, M 2 and M 3 , are affected. 

(iv) Thus, for the Schwinger model, there remains a few percent discrepancy 
in the coefficients Mi and M 2 compared to bosonization results. 

(v) The variational calculation yields the 'pion' wave function to a high accu- 
racy. The endpoint behavior is reproduced exactly (in leading order in m). The 
behavior in the intermediate region, < x < 1, gets improved as can be seen 



from Fig. |3.3| below. 



In Table 8.15| we summarize our optimum final results for the 't Hooft model 



(a = 0) and compare with results that have been obtained previously. 



Upon inspection of Table |3.15| one finds a very good overall consistency of the 
results. All given values are in good agreement within error bars. The scale is set 
by the variational results which are the most accurate ones. Within error bars, 
they are matched by our numerical results as well as by previous analytical and 
numerical calculations. 

As far as the latter are concerned, a few remarks are in order. The coefficients 
Mi displayed in the last line of Table |S.15| are obtained by performing polynomial 



fits to the values calculated numerically in [|192|| (which basically agree with those 



from 't Hooft's original calculation displayed in Fig. 5 of |127| |). The averaged 
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Table 8.15: The expansion coefficients Mi for the 't Hooft model (a = 0). 
The errors of our results obtained within 2-particle Tamm-Dancoff (2PTD) ap- 
proximation are estimated by comparing the last two lines in Table \$.2\ (for the 
variational method) and the different polynomial fits $3.72 - ^.82j) (for the nu- 
merical results). The numerical inaccuracies given in the last line were estimated 
by us using the polynomial fit method. 



Mi 

variational 2PTD 2tt/a/3 = 3.627599 
numerical 2PTD 3.62758 ± 2 • 1(T 5 



't Hooft BBS | 



Burkardt g80] 
Li GJ1 



Li et al. p2 | 



2vr/^3 = 3.627599 
2ir/V3 = 3.627599 
3.64 ± 0.06 
3.64 ±0.03 



M 2 

3.581055 ±3- 10~ 6 
3.5829 ±3- 10~ 4 

3.5812 

3.60 ±0.06 



Ms 

0.061793 ±3-10- 
0.064 ± 1 • 10~ 3 



0.04 ±0.04 



result of the fits is shown including an estimate of the errors. It should also be 
stressed that the data in [ 192 ] are mainly obtained for m > 0.5 whereas the bulk 
of our values for the Mj are obtained for m < 0.1. Nevertheless, the agreement 
with our results is satisfactory within error bars. 

From Table |8.15| it is obvious that all numerical results favor the first order 
't Hooft-Bergknoff value of 27t/a/3. This agreement is particularly gratifying for 
the data of ||192|| which were not obtained via light-cone techniques but within 
ordinary quantization and even within a different gauge, namely axial gauge. The 



agreement is therefore highly non-trivial ||193 
For the Schwinger model 



a 



1), we summarize our findings in Table |8.16 



For comparison we have also listed analytical results obtained via bosonization 
techniques and data which we extracted from polynomial fits to numerical light- 
front Tamm-Dancoff ||187|| , DLCQ ||183|| and lattice results ||194 ]. The contents of 
Table |B.16| are graphically displayed in Fig. |8.2| , where we have chosen rescaled 
units as in ||176| . 

As already stated, due to the anomaly the numerical errors for the Schwinger 
model are at least an order of magnitude larger than for the 't Hooft model. 
In addition, the 2% discrepancy between our light-cone results and the analytic 
bosonization results does not get resolved. However, we have shown that the 
discrepancy is not due to (i) inaccuracies in the wave function or (ii) neglect 
of contributions from the next higher Fock sectors (four and six particles). It 
should be mentioned that other numerical methods (DLCQ, lattice) are by far 
too inaccurate to distinguish between the 't Hooft-Bergknoff and the bosonization 
value. We will come back to these issues in a moment. 
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Figure 8.2: The rescaled l pion' mass M = M/yl + m 2 as a function of 
the fermion mass m. The short-clashed curve represents a l phenomenologicaV 
parametrization for M with M 2 = 1 + Mim + 4m 2 , which becomes exact for very 
small and very large m and thus smoothly interpolates between the strong and 
weak coupling regions. The long-dashed curve is our second-order and the solid 
curve our third-order result. As expected, mass perturbation theory breaks down 
as m becomes of order 1. The crosses are the lattice results j \194l , the diamonds 
the lattice results J195|/ , which go down to comparatively small masses, however 
with large errors. The circles are light-front Tamm-Dancoff results /|187f , whereas 
the squares ftL83H and triangles \196j , as quoted in \1 76\] , are DLCQ results. The 
values corresponding to the triangles are not included in Table [Q6 . The 2% 
discrepancy in Mi is invisible within the resolution of the figure. 
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Table 8.16: The expansion coefficients Mi for the Schwinger model (a = 1). 
The errors for the variational 2-particle Tamm-Dancoff (2PTD) approximation 
are estimated by comparing the last two lines in Table [5]^. The errors of our 
numerical calculations are obtained from comparing polynomial fits of different 
order to the numerical results of Tables 8.11 and p.!2\ . Errors of other results are 
given where they could be estimated analogously. 



Mi 



variational 2PTD 
numerical 2PTD 
numerical 6PTD 



2tt/V3 = 3.627599 
3.6268 ±8 • 10~ 4 
3.6267 ±4 • 1(T 4 



numerical 4PTD (T87|3.62 ±0.07 
DLCQ |H 3.7 ±0.2 

lattice |T9l 3.5 ±0.2 

2e< = 3.562146 
2e 7 = 3.562146 



J19J 
Adam flTf] 
Fields et al. |7| 



M 2 

3.308608 ±4-10" 
3.32 ±0.02 
3.22 ±0.02 
3.2 ±0.3 
3.5 ±0.3 
3.7 

3.3874 
3.387399 



Ms 

0.348204 ±1 
0.28 ±0.06 
0.5 ±0.1 
0.3 ±0.2 

0.02 



10~ 



Table 8.17: Lattice results for the first order coefficient Mi and the (negative 
of the) condensate, —(0\'i/}ip\0) = Mi/4n. The results are quoted in chronological 
order. They should be compared with the bosonization results, Mi = 3.562, 
— <O|-0-0|O> = 0.283 and the 't Hooft-Bergknoff values, Mi = 3.628, -(0|^|0) = 
0.289. 



P3 



am 



[169] 



[1681 



Mi 


3.42 


3.644 


3.48 


2.97 


3.33 


3.77 


-<0|^|0> 


0.27 


0.290 


0.28 


0.24 


0.26 


0.30 



As is well known, ||1 76| , |187| , |197|| , the DLCQ data are comparatively inaccurate 
for small fermion mass m (i.e. large coupling). This is due to the fact that the 
dominating feature of the light-cone wave function, its endpoint behavior, is not 
very accurately reproduced using an equally-spaced momentum grid. The poor 
convergence of DLCQ for small m has recently been overcome by incorporating 
the exact endpoint behavior [|197|| . 

The lattice data generally suffer form the same disease of having rather large 
error bars for small m. The chiral limit, m — 0, can only be reached via ex- 
trapolation in a very inaccurate manner unless one uses the Schwinger result 
M 2 (m = 0) = 1 as a bias A representative collection of lattice results 

for Mi is given in Table |8.17| . For the convenience of the reader we also list the 
condensate (0|^|0) = -M 1 /Air. 

The values for the Schwinger model condensate in Table |8.17| should be com- 
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pared with our results. If we assume that formulae flS.17| ) and ( |3.18| ) are directly 
applicable to our calculation we get 



(0|^|0)i = -|^ = -^= = -0.28868, (8.84) 
1 / Mo 

<0|^|0> 2 = -^V^^ = -°- 28015 ' ( 8 - 85 ) 
while the standard result is 

(0|^|0) = -— = -0.28347 . (8.86) 

As already mentioned in Section 8.2, the few percent discrepancies in Mi and 
M2 immediately affect the condensate. Note that the condensate value obtained 
from the second order coefficient is closer to the bosonization result, the relative 
error being 1.2 % compared to the 1.8 % at first order. Furthermore, the exact 
result ( |3.86|) lies between our first and second order values, ( |3.84|) and (|8.85|) , 
so that their mean value, (Ol^^lO) 1/2 = —0.28442, is much closer to ( |8.86| ), the 



error being only 0.3%. A similar reduction of errors is actually at work if one 
considers the 'pion' mass-squared of the Schwinger model as a function of m. The 
first order overshoots the bosonization value while the second order contribution 
is too small. Adding both one gets closer to the exact result, at least if the 
mass m is not too tiny. For the value of Fig. |8.1| , m = 0.1, one finds to second 
order M 2 (0.1) = 1.39585, while, in this order, M| os (0.1) = 1.39009. The relative 
difference is only 0.4%. Possible resolutions of the discrepancies between (|8.84|) , 
( |8.85|) and (|8.86| ) will be discussed at the end of this section. 

Apart from the 'pion' mass squared, the eigenvalue in the 't Hooft-Bergknoff 
equation, we have also calculated the associated eigenf unction, the light-cone 
wave function. The latter has been obtained with high accuracy as can be seen 
from the good convergence of the alternatively defined mass squared coefficients, 



Mi towards the variational estimates, Mj (see Section ^4). If we denote 4>o{x 



x(l — x), our most accurate variational ansatz for the wave function can be written 

as 

A a, b, c, d] = 4 + a <^ +1 + b <^ +2 + c <^ +3 + d ■ (8-87) 

According to our mass perturbation theory, each of the variational parameters is 
expanded in powers of the fermion mass m, the coefficients being denoted fix, (32, 
03, ax, a>2, etc. From 't Hooft's endpoint analysis, j3 is known exactly. It turns 
out to be independent of the anomaly (a) and thus is the same for the 't Hooft 
and Schwinger model. In Table |8.18| we compare the exact expansion coefficients 
of j3 with their variational estimates. The best estimates for the other variational 
parameters are listed in Table for both the 't Hooft and Schwinger model. 
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Table 8.18: Comparison of the exact expansion coefficients of the endpoint 
exponent (3 with their best variational estimates for both the 't Hooft (a = 0) 
and the Schwinger model (a = 1). 



Pi & 

exact VS/n = 0.55132890 -A/10 = -0.05513289 

variational (a = 0) y/3/n = 0.55132890 0.00051 -0.053 

variational (a = 1) y/3/n = 0.55132890 0.00063 -0.053 



Table 8.19: Best estimates for the variational parameters in the light-cone wave 
function for both the 't Hooft (a = 0) and Schwinger model (a = 1). 





a i 


a 2 


bi 


b 2 


Cl 


c 2 


di 


d 2 


a = 


0.92 


1.19 


-2.2 


-3.0 


7.0 


9.2 


-9.9 


-11.8 


a = 1 


1.42 


1.03 


-2.7 


-3.1 


8.8 


9.6 


-12.5 


-12.1 



The errors in the coefficients are comparatively large and growing from the 



left to the right in Table |3.19| . Good numerical convergence is evident for the 
coefficients of a. For the other coefficients our method does not provide enough 
iteration steps to make convergence explicit. However, the sensitivity of the mass 
expansion coefficients Mj and the wave functions on the parameters b, c, d is very 
weak (see Fig. |8~3j ). 

Let us finally investigate possible sources for the few-percent discrepancies 
in the Schwinger model results like for the condensate values, (|8.84j ) and flS.85|) 



vs. ( |3.86| ). First of all, it has to be reemphasized that the light-cone calculations 
are conceptually rather different from the usual treatment based on bosonization 
within standard equal-time quantization. We are studying the lowest lying me- 
son as a relativistic fermion-anti-fermion bound state. The fermionic degrees of 
freedom are explicitly taken into account as they define our Fock basis. In the 
bosonized theory, which is a sine-Gordon model with the Lagrangian [ |199| , |164 



1 i 

£ = -d^d^cf) - -ulcj) 2 - cmu, cos V4^> , (8.88) 

_ _ 

the 'pion' is the elementary particle described by the scalar field <fi. One calcu- 
lates its mass by perturbation theory in the 'coupling' cm/io, c = e 7 /27r. In the 
bosonized theory, this 'coupling' is actually dependent on the normal-ordering 
scale used to renormalize the theory. The standard choice is the most natural 
one, namely the Schwinger boson mass Ho, which explicitly appears in the La- 
grangian QS.88p . It is, however, not very hard to determine the general scale 
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Figure 8.3: The light-cone wave function of the 't Hooft model l pion' for m 
= 0.1. The solid curve represents the result from 't Hooft 's original ansatz, the 
dashed curve our best result (with maximum number of variational parameters) . 
At the given resolution, however, the curves of all extensions of 't Hooft's ansatz 
(a, b, c, d) lie on top of each other. 




dependence of quantities calculated within perturbation theory. All one needs is 
Coleman's re- normal-ordering prescription ||199|| which for the case at hand is 



COS V^47T0 = —Np COS V^47T0 . 



5.89) 



Here, iV Mo and denote normal-ordering with respect to /io and /i, respectively. 
From ( |3.89| ) one derives the following identity for the boson mass squared (setting 
fM) = 1), 



M = 1 + Mi m + M 2 m 2 + M 3 m 3 + . . . 

= 1 + M[ fxm + M' 2 jj, 2 m 2 + n 3 m 3 + ... . 



190) 



In this expression, the Mj denote the coefficients calculated with the standard 
normal-ordering scale /io, where all tadpoles can be set to zero, while the M[ are 
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calculated with normal-ordering scale /i, where one has non-vanishing tadpole 
contributions. In order that the 'pion' mass be scale independent one has to have 
that (upon reintroducing fi ) 

M' n = M n . (8.91) 

We have checked this identity explicitly in the first two non-trivial orders of mass 
perturbation theory with the Lagrangian (|3.88| ). It is also interesting to note that 
the scale dependence of the coefficients might equivalently be associated with a 
multiplicative rescaling of the fermion mass m by defining m! = fim and using 
m! as the expansion parameter in ( |3.90 ). 



Now, the upshot of all this is the following. As we are working entirely in the 
fermionic representation, we do not a priori know, to which normal-ordering scale 
of the bosonized theory our results correspond. If we assume that it is not the 
Schwinger mass /io but a different scale fi, we can match the first order coefficients 
by choosing 

H = = , (8.92) 

where M[ now denotes our light-front result. Quite a similar point of view has 
been taken by Burkardt in his paper on light-cone quantization of the sine-Gordon 
model |[20Q| . However, as is obvious from ( |8.91|) , the change of scale ( |8.92|) prop- 



agates systematically through all coefficients. In particular, if M[ is larger than 
the corresponding bosonization result, the same has to be true for all M[. Our 
coefficient M' 2 , however, already fails this requirement. We thus conclude that a 
possible scale dependence alone cannot explain the discrepancy. 

An additional source of error could still be contributions from higher Fock 
sectors. We have shown numerically that the first few low orders have very little 
impact on the expansion coefficients Mj. The contribution from each Fock sector 
might thus be exponentially small. The whole series, however, if it could be 
summed up, might add something seizable. Let us assume, for example, that the 
contribution of higher Fock sectors to M 2 has the following formP] 

00 

5M 2 = J2 c 2n{m) , (8.93) 

n=2 

where 2n is the label of the 2n-particle sector. If we had c 2n [m) = e~ 2nm , sum- 
mation of the series would yield 5M 2 = l/2m + 0(1) and thus effectively a 
contribution to the first order in m. As we are not able to calculate the real 
high order behavior of the series in ( |8.93|) , we cannot, at the moment, make any 
definite conclusion about a possible relevance of such summation effects. All we 



3 We thank M. Burkardt for suggesting the following example. 
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can say is, that if such effects are present, they must be numerically small, i.e. at 
the few percent level. 

As a first step towards getting some better control of higher particle sectors 
one would like to develop some tools to estimate their contributions analytically. 
This amounts to determine the coefficients C2 n (m), at least for small n. So far, 
higher Fock states have only been included numerically. 

It might also be interesting to calculate excited states. For the Schwinger 
model, this again requires the inclusion of higher Fock components as pointed 
out by Mo and Perry ||187| . For the 't Hooft model, on the other hand, one can 
stay in the two-particle sector. One might then try to modify the variational 
procedure, basically by taking into account the growing number of nodes, with 
the aim to match the low energy spectrum with the high energy behavior derived 
by 't Hooft. Work in these directions is underway. 



Chapter 9 



The Pion Wave Function in the 
NJL Model 



In this last chapter we will finally study 'realistic' light-cone bound-state equa- 
tions and wave functions. To this end, we enter the world of 3+1 dimensions 
and investigate the pion wave function within a particularly suited model due to 
Nambu and Jona-Lasinio (NJL) [ [201| . 



The first derivation, analysis and solution of a light-cone bound-state equa- 
tion appeared in 't Hooft's original paper on what is now called the 't Hooft 



model ||127|| . We have discussed this model at length in the last section where we 
also rederived 't Hooft's solution. Interestingly, 't Hooft did not use the light- 
cone formalism in the manner we presented it and which nowadays might be 
called standard. This amounts to deriving the canonical light-cone Hamiltonian 
and setting up the associated bound-state equation by projecting the Schrodinger 
equation on the different sectors of Fock space (cf. Chapter 6). Instead, he started 
from covariant equations, namely the Schwinger-Dyson equations for the quark 
propagator (or self-energy), and the Bethe-Salpeter equation for the bound state 
amplitude, which needs the quark self-energy as an input. The light-cone bound 
equation was then obtained by projecting the Bethe-Salpeter equation onto hy- 
persurfaces of equal light-cone time thereby spoiling explicit covariance. In this 
way, one avoids to explicitly derive the light-cone Hamiltonian, which can be a 
tedious enterprise in view of possibly complicated constraints one may have to 
solve. This makes it worthwhile to have a closer look at this way of proceeding. 

In a condensed notation, the Schwinger-Dyson equation for the propagator 
can be written as 

S = S + S ZS , (9.1) 

where 5* is the full propagator, 
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So the free propagator, 

So(p) = , (9.3) 

p — mo 

and E the quark self-energy. The Bethe-Salpeter equation for a quark-antiquark 
bound state (a meson), on the other hand, is 

Xbs = S^Kxbs , (9.4) 

with Si, S2 the full propagators of quark and anti-quark, K the Bethe-Salpeter 
kernel, defined through the interaction one is considering, and xbs the Bethe- 
Salpeter amplitude. From the latter, one obtains the light-cone wave function 



via integration over the energy variable p [|202|, [203|, |204|, |205 



0lc(p) = J f^XBs(p) , P = P±) • (9-5) 

The program outlined above, deriving the light-cone bound-state equation from 
the covariant equations, when applied to QCD, is of course ambitious but should 
still be feasible in view of the wealth of results already obtained in studying 
the Schwinger-Dyson equations of QCD ||206| , [207|| . In a first step, which is 
understood, say, as a feasibility study, a simpler approach will be pursued, namely 
a calculation within a particularly suited framework, the NJL model. In the 
language of the Schwinger-Dyson framework, this amounts to using a model gluon 
propagator that is constant in momentum space. Obviously, such a propagator 
does not lead to confinement, which however, as pointed out in the introduction, 
does not necessarily prohibit a reasonable description of hadron structure. The 
NJL model has been intensively studied within the last decade, and a wealth of 



results has been obtained (see the recent reviews |208| , |209| , [210| , |211|| ). To the best 
of our knowledge, however, light-cone wave functions never have been calculated 
within the model. 



9.1 The Condensate 

In its simplest form, the NJL model has a chirally invariant four-fermion inter- 
action, which can be imagined as the result of 'integrating out' the gluons in the 
QCD Lagrangian. Within the instanton model of QCD, this can indeed be done 



and leads to an NJL type interaction [IB|. For two massless (m = 0) quark 
flavors (u and d), the NJL Lagrangian is 

£ NJL = i>%W> + G [{^f + (^7st^) 2 ] • (9.6) 

In what follows we will mostly suppress flavor (and color) degrees of freedom to 
keep things as simple as possible. As the coupling G has negative mass dimension 
mass -2 ), the model is not renormalizable and needs a cutoff as a parameter 
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that has to be fixed by phenomenology. Above a critical coupling, G > G c , a 
fermion mass is dynamically generated so that chiral symmetry is spontaneously 
broken corresponding to a second order phase transition. The dynamical fermion 
mass is proportional to the fermionic condensate in the vacuum as was briefly 
mentioned in the introduction, Section 1.3. 

These results are obtained via a self-consistent mean-field solution of the 
Schwinger-Dyson equation resulting in the gap equation, 

£ = m = 2iGtrS(x = 0) = — — / = -2G<0|^|0) m . (9.7) 

ZilX J K 77? 

This equation determines the dynamically generated mass m. The integral on 
the right-hand-side is of course divergent and has to be cut off (see below). 

The mean-field solution has a very intuitive explanation. One essentially 
argues that the main effect of the interaction is to generate the mass of the 
quarks which become quasi-particles that interact only weakly. Neglecting this 
interaction entirely, one can view the process of mass generation like in Section 
5.2.3 as the transition of quarks with mass m = to mass m resulting in a mass 
term mipip in the Lagrangian ( |9.6|) or in the Hamiltonian. The condensate is 
obtained via the Feynman-Hellman theorem by differentiating the energy density 
of the quasi-particle Dirac sea, 

o 

d d f , f d 2 k±m 2 + k\ 



oo 



m f dk 



8tt 3 / k+ 



16tt 3 k+ 
d 2 k ± , (9.8) 



o 



Again, as it stands, the integral is divergent and requires regularization. In the 
most straightforward manner one might choose m 2 /A < k + < A and |kj_| < A, 
so that the condensate becomes 

A A 2 

. , m f dk + f ,,,n s m . 9 , A 2 

m?/K 



Plugging this result into the gap equation (9.7) one finds for the dynamical mass 
squared, 

m 2 (G)=A 2 exp(-|^) . (9.10) 

The critical coupling is determined by the vanishing of this mass, m(G c ) = 0, 
and from ( |9.10|) we find the surprising result 



G c = 



(9.11) 
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This result, however, is wrong since one knows from the conventional treatment of 
the model that the critical coupling is finite of the order ir 2 /A 2 , both for covariant 



and non-covariant cutoff ||201|| . In addition, it is quite generally clear that in the 
free theory (G = 0) chiral symmetry is not broken (as mo = 0) and, therefore, 
this should not happen for arbitrarily small coupling, either (cf. Chapter 7). The 
remedy is once more to use an information from the ordinary calculation of the 
condensate. We translate the non-covariant, but rotationally invariant, three- 



vector cutoff, |k| < A, into light-cone coordinates ||81|| , which leads to 

2 

771 

< k 2 ± < 2Ak + - m 2 - (k + ) 2 , — < k + < 2A . (9.12) 

Note that the transverse cutoff becomes a polynomial in k + . The fc_i_-integration 
thus has to be performed first. Renaming 2 A —>■ A and introducing the longitu- 
dinal momentum fraction, x = k + /A, ( fj.l2j ) becomes 



< k\ < A 2 x(\ -x) -m 2 , m 2 /A 2 < x < 1 - m 2 / A 2 . (9.13) 
Surprisingly, this is just the often used 'invariant-mass cutoff', 

M 2 = < A 2 , (9.14) 

x{± — x) 

with M5 2 the invariant mass-squared (for a two-particle system) introduced in 
Chapter 6. Note that the cutoff is also invariant under the exchange x <-> 1 — x. 

For the condensate (|9.8|) the cutoff ( |9.12j ) yields an analytic structure different 
from rOI), 



(0|^|0) m = -^f2A 2 -m 2 ln^ , (9.15) 
Sir 2 \ m 2 J 

where we have neglected sub- leading terms in the cutoff A. From (|9.15|) , one 
infers the correct cutoff dependence of the critical coupling, 

G c = ^ . (9.16) 

The moral of this calculation is that even in a non-renormalizable theory like the 
NJL model, the light-cone regularization prescription is a subtle issue. In order 
to get a physically sensible result the transverse cutoff has to be /c + -dependent. 
Clearly, this dependence cannot be arbitrary but should be constrained from 
dimensional and symmetry considerations. For renormalizable theories, such ar- 
guments have been given by Perry and Wilson |[212|| . In the example above, it 



was rotational invariance that solved the problem. 

The condensate ( |9. 15[ ) was originally obtained in [51], where, however, a 



slightly more complicated cut-off was used. In that work, the condensate was 
defined covariantly in terms of the fermion propagator at the origin, S(x = 0) 
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like in ( |9.7| ). Performing the energy integration over k~ with the appropriate 
cutoff again leads to (|9.15|) . 

In the NJL model with its second-order phase transition of mean-field type, 
the usual analogy with magnetic systems can be made. Chiral symmetry cor- 
responds to rotational symmetry, the vacuum energy density to the Gibbs free 
energy, and the mass m to an external magnetic field. The order parameter 
measuring the rotational symmetry breaking is the magnetization. It is obtained 
by differentiating the free energy with respect to the external field. This is the 
analogue of expression ( |9.8| ) as derived from the Feynman-Hellmann theorem. 



9.2 The Light-Cone Bound State Equation 

Once the physical fermion mass m is known, it can be plugged into the Bethe- 
Salpeter equation (|9.4|) which in ladder approximation reads 



(*) = / ^iS(k)iS(k-P) Jj0pj ^K(k,k) X Bs(k) , (9.17) 



with P denoting the bound-state four-momentum. On the left-hand-side, the 
projection onto x + = (i.e. the ^"-integration) has already been carried out. 
On the right-hand-side, the two integrations over k~ and k~ still have to be 
performed. Whether this can easily be done depends of course crucially on the 
kernel K, which in principle is a function both energy variables. For the NJL 
model, however, one has (after a Fierz transformation), 

K(k, k) = 2 75 ® 75 - 7^75 ® 7 M 75 , (9.18) 

i.e. the kernel is momentum independent due to the four-point contact interaction! 
Thus, the ^"-integration immediately yields 4>lc, and the ^"-integration can be 
performed via residue technique and is completely determined by the poles of the 
propagators, S(k) and S(k — P). As a result, one finds a non-vanishing result 
only if < k + < P + , and one of the two particles is put on-shell, e.g. k 2 = m 2 , 
as already observed by Gross |[213 |. 



The upshot of all this is nothing but the light-cone bound-state equation, 
which reads explicitly 



'x, ki 



2G (# + m ) 75 (#-/>+ m ) 
x(l-x) M 2 -M 2 
1 <■ d 2 k 



x I dy I — tr 



+ 



750LC(2/, k|_) 0A(2/,kj 



x(l -x) M 2 — M 2 
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x 



dy 



d 2 k< 



-tr 



7 M 750Lc(y,k±) 0A(y,k±) • 



(9.19) 



Here we have defined the longitudinal momentum fractions x = k + /P + , y = 
k + /P + , the on-shell momentum k = (k~ , kj_, k + ) with k = {k\ + m 2 )/k + , and 
the bound state mass squared, M 2 = P 2 , which is the eigenvalue to be solved 
for. 6 A (x,]<.±) denotes the invariant mass cutoff |yj). 

Now, while ( |9.19|) may appear somewhat complicated it is actually very sim- 
ple; indeed, it is basically already the solution of the problem. The crucial ob- 
servation is to note that the two integral expressions are mere normalization 
constants, 



C A = 



D*P» = 



dy 



dy 



d 2 k A 

d 2 ~k ± 
"8^ 



■tr 



-tr 



750lcO, k ± ) A (y,k_L) , 
7 M 750Lc(y,k_L) 9 A (y, k±; 



(9.20) 
(9.21) 



Thus, up to normalization, the solution of the light-cone bound-state equation 

2GC A (j£ + m ) 75 (j£_/>+ m ) 



U9|) is 



x, k i 



M 2 - M 2 



+ 



x(l — X) 

GD A $ + m)f*y s (t ~ f + m) 



x(l — x) 



M 2 - M 2 



(9.22) 



As a first check of our bound state wave function ( |9.22| ) we look for a massless 
pion in the chiral limit. To this end we decompose the light-cone wave function 
according to |J, 

0LC = 03 + 0p75 + 0^7^75 + + 0f <V • ( 9 - 23 ) 

Multiplying ( |9.22|) with 75, taking the trace and integrating over k we find 

GC A 



° A = dX J dk± x(l-x)(M 2 -M 2 ) 6A{x ^ 

M 



+ 



2tt 3 



dx / d k 1 0a fx, ki 



(9.24) 



In the chiral limit one expects a solution for M = 0, the Goldstone pion. In this 
case one obtains 

G 



2tt 3 



dx J d 2 k±9\(x, kj 



(9.25) 
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This is exactly the gap equation (|9.7| ) upon performing the ^"-integration there 
(with an invariant mass cutoff ( j9.13|) understood in both identities). Note once 
more the light-cone peculiarity that the (Fock) measure in ( |9.25| ) is entirely mass 
independent. All the mass dependence, therefore, has to come from the (invariant 
mass) cutoff. Otherwise one will get a wrong behavior of the dynamical mass m 
as a function of the coupling G, as was the case in ( p.lOj ) . 

With this in mind, we have found that the Goldstone pion is a solution of the 
light-cone bound-state equation exactly if the gap equation holds. This provides 
additional evidence for the self-consistency of the procedure. The same results 
are well known from the covariant treatment ||201|| , and it is gratifying to get 
them also within the light-cone formalism. The deeper reason for the fact that 
the quark self-energy and the bound state amplitude satisfy essentially the same 
equation [[214 ], is the chiral Ward identity relating the quark propagator and the 
pseudo-scalar vertex | |215[ . 

Our next task is to actually evaluate the solution ( |9.22j ) of the bound-state 
equation. 0lc is a Dirac matrix and therefore is not yet a light-cone wave function 
as defined in Chapter 6. The relation between the two quantities has been given 
by Liu and Soper ||204| . 

2P + ^/x~x~ ip{x, k ± , A, A') = u(xP + , k ± , \) 1 + <pi jC (k)-i + v(xP + , -k ± , A') , (9.26) 

where we have denoted x = 1 — x to save space and included a factor (xx) 1//2 to 
match with the conventions of Brodsky and Lepage [^, [216j , |217|| . A somewhat 
lengthy calculation, using the spinor identities of App. A, yields the result 



ip(x, k ± , A, A') 



2GP' 



M 2 



Ml 



2CU 
M 



(9.27) 



with the arguments of the spinors u and v suppressed. At this point we have 
to invoke another symmetry principle. Ji et al. have pointed out p3f that the 



spin structure (uTv) should be consistent with the one obtained form the instant 
form spinors via a subsequent application of a Melosh transformation | |218| and a 
boost. Using this recipe, one obtains the following relation between the constants 
C A and D A , 

2C A /M = -D A = N A /2GP + . (9.28) 

As a result, the spin structure in (|9.27| ) coincides with the standard one used in a 
number of works [[JT], |219j , |9"2"| , |220| , |9"B| ] . The NJL wave function of the pion finally 
is 

N(A) 



ip(x,kj_, A, A') 



M 2 



M 2 



u x (M + P) l5 v y 6(A 2 -M 2 



(9.29) 



N(A) is a cutoff dependent normalization defined in (|9.28|) , and the spin (or 
helicity) structure is given by 

u(xP + , k ± , A) (M + f) 75 v(xP\ -k ± , A') = 
1 



' XX 



A 



mM + m - k 1 + M z xx 



5x,-x> ~ k- X (M + 2m)8 xx , (9.30) 
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where we have defined k\ = k 1 + 2A&; 2 . The first term with spins anti-parallel 
corresponds to L z = 0, the second one (with spins parallel) to L z = ±1. It has al- 
ready been pointed out by Leutwyler that both spin alignments should contribute 
to the pion wave function ||221| , |222|| . Note that the latter is a cutoff dependent 



quantity. This is necessary in order to make the wave function normalizable. In 
other words, it guarantees the boundary conditions ( |6.23| ) which state that the 
wave function drops off sufficiently fast in x and k±. 

As we are interested in analyzing the quality of a constituent picture, we 
approximate the pion by its valence state, 

f 1 dr f d 2 k i 

\k:P} = J2 m V / ^4^,k J ,,A,A / )|gg:g,k ±> A,A / ) (9.31) 

XX \/%(*- X) J -LOT 



which should be compared with the general expression (|6.13j ). The normalization 
of this state is given by ( 6.10|) , 



(tt : P'\n : P) = 16tt 3 P + 5 3 (P - P') . (9.32) 

As usual we work in a frame in which the total transverse momentum vanishes, 
i.e. P = (P + ,P± = 0). Expression ( |9.32| ) yields the normalization ( |6.19| ) of the 
wave function, 

Y,f dx I ^l^(x,k ± ,A,A')| 2 = l. (9.33) 
xx' ^° ^ 71 

It is of course a critical assumption that the probability to find the pion in its 
valence state is one. In this way we enforce a constituent picture by fiat, and 
it is clear that such an assumption has to be explicitly checked. An alternative 
way of proceeding would be not to fix the normalization but to use -AT (A) as an 
additional parameter that has to fixed by phenomenology. We will, however, use 
( |9.33| ) and check whether a valence (or constituent) picture makes sense. 



9.3 Observables 

With the light-cone wave function at hand, we are in the position to calculate 
observables. As a technical prelude, however, we first have to determine its 
normalization. To keep things simple we consider the chiral limit, M — 0. We 
write the pion wave function as a matrix in helicity space, 

(9.34) 

The diagonal terms correspond to parallel spins, the off-diagonal ones to anti- 
parallel spins. 
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The normalization ( |9.33| ) determines N as a function of A. Using the explicit 
form of the invariant-mass cutoff, (|9.13| ), one finds 



in' 



l-e z 



dxxil 



X) 



A 2 (x) 



where we have defined 
,2 . 



A 2 Or) 



m 2 /A 2 
A 2 x(l 



x) 



m 



A 2 



x(l 



x 



(9.35) 

(9.36) 
(9.37) 



As an additional approximation we will neglect terms of order e 2 , i.e. we assume 
that the cutoff is large compared to the constituent mass. From the standard 
values, A ~ 1 GeV, m ~ 300 MeV, we expect that this assumption should induce 
an error of the order of 10%. The advantage is that it allows for a simple analytic 
evaluation of all the integrals we will encounter. Furthermore, the leading order 
will be independent of the actual value of the constituent mass. From ( |9.35|) 
thus finds the normalization, 



one 



\N\ l 



240vr 2 



[l + 5e 2 + 0(e z 



Choosing N real we obtain to leading order, 

N = 2Vl5n/A 



(9.38) 



(9.39) 



This fixes our normalization constant. We proceed by calculating the pion electro- 
magnetic formfactor. It is defined by the matrix element of the electromagnetic 
current J£ m between pion states, 



(tt : P|J£Jtt : P) = 2{P + P')"F(Q 2 



Q 2 



(P-P 



l\2 



(9.40) 



Considering fi = + in a frame where P = (P + , 0) and P = (P + , q ± ) one is led 
to the the Drell-Yan formula |2|, ||, 



^(qi) 



AA' 



dx 



d 2 k A 

16-7T 



■r(x,k' ± ,X,X')^(x,k ± ,X,X') . 



(9.41) 



The transverse momentum of the struck quark is = kj_ + xq ± . Note that 
the formula ( |9.41|) with its overlap of two wave functions on the right-hand-side 
is rather similar to the non-relativistic result ||82|| . If one sets the momentum 
transfer q ± equal to zero, one is left with the normalization integral (|9.33|) , so 
that the form factor is automatically normalized to one. Plugging the pion wave 
function (|9.34|) into ( |9.41|) , one finds to first order in q^_, 



l-o 2 ^ 



20 



e 2 + 0(e 4 



(9.42) 
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Inserting the normalization ( |9.38| ), the formfactor becomes 



F(qi) = l-^ql 



1- |e 2 + 0(e 4 ) 



(9.43) 



This leading order result already determines a particular observable, namely the 
charge radius of the pion, which is defined via the slope of the formfactor at zero 
momentum transfer, 



(r|) 
6 



ql + 0(qi). 



To leading order in e 2 the pion charge radius thus becomes 

(rl) = 18/A 2 , 



(9.44) 



(9.45) 



and is thus entirely determined by the cutoff A. We might fix the latter by using 
the experimental result |224|1 



(rl) = 0.44 fm 2 . 

The cutoff thus becomes (to leading order in e 2 ) 

18 



(rl) 



1.3 GeVy . 



(9.46) 



(9.47) 



Our preferred way, however, to fix the cutoff is to use the pion decay constant 
= 92.4 MeV [[225 ], which is known to a higher accuracy than the charge radius. 
Though this does not matter in view of our present approximations, it seems to 
be advantageous at least in principle. f n is given by the 'wave function at the 
origin', ip , defined as 



dx 



d 2 k A 

16/T : 



^(x,k ± ,T,|)^(A 2 



Ml) 



3NA 2 
2048vr 



0(e 2 



(9.48) 



Note that only the L z = component contributes. Using the normalization ( |9.38j ) 
this becomes 

ip = — — A ~ 0.022 A . 9.49 
512 



The pion decay constant has been obtained from weak w decay |82[] where only 
the pion valence wave functions contributes to the decay matrix elements. There- 



fore, the following result [fB2| , pl7]| provides an exact constraint on the pion wave 
function, 

ifj = ^= . (9.50) 



c 
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Figure 9.1: The pion formfactor squared vs. momentum transfer q 2 = q\. The 
full line is the monopole fit of f?2§, \F\ 2 = n/(l + gi(r 2 )/6) 2 with n = 0.991, 



(r 2 ) = 0.431 fm 2 ; the dashed line is the same fit with our values, n = 1, (r 2 ) 
0.50 fm 2 . The agreement is consistent with the expected accuracy of 10%. 
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Comparing with ( |9.49| ) we find that the cutoff is fixed through 

256 



A 



9^ 



f n ~ 12.7 f n - 1-2 GeV . 



(9.51) 



Having fixed the cutoff, we are in the position to actually predict the pion charge 
radius via (|9.45[) . The result is 



(r 2 ) = = 0.50 fm 2 = (0.71fm) 2 



(9.52) 



which differs from the experimental value ( |9.46| ) by 13% and therefore is consistent 
with the expected accuracy 



This is also evident from Fig. |9.1| where we compare with the experimental 
results of |[224j|. The full and dashed line correspond to monopole fits of the form 



n 



;i + qi(r 2 )/6) 2 ' 



(9.53) 



with n = 0.991, (r 2 ) = 0.431 fm 2 dH and our result fl9T52D , n = 1, (r 2 ) = 0.50 
fm 2 , respectively. 
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If we compare our result fl9.51| ) for the cutoff with the value for the cutoff A x 
below which chiral effective theories make sense, A x = 4nf n ~ 12. 6f n [[226j| , we 
find that the agreement is surprisingly good. Clearly, this is somewhat accidental, 
as we do not expect an accuracy better than 10%. This can be seen explicitly 
again by alternatively using the pion charge radius to fix /„.. Plugging ( |9.47| ) into 
the general relation ( |9.51| ) between cutoff and and solving for we predict 

Q^/5 

/tt = ——A ~ 0.077 A ~ lOOMeV , (9.54) 
256 

which indeed is accurate to within 10%. 

We proceed by calculating the mean transverse momentum, 



^l) = Y.f dx Ijh^ l^k ± ,A,A')| 2 #(A 2 -M ( 



(9.55) 



A A 

,2 



To leading order in e this becomes 

2240tt 2 



(kl) = , (9.56) 



or, upon inserting the normalization 



(kl) = 4 A " - (°- 43 GeV ) 2 • ( 9 - 57 ) 
28 

The r.m.s. transverse momentum is thus of the order of 400 MeV, which is a 



reasonable value ||221|| . In particular, this value shows that the pion is a highly 
relativistic system as the momenta of its (valence) constituents exceed the pion 
mass. 

The inverse of the r.m.s. momentum can be identified with the pion core 
radius, for which one finds 

R w = (k 2 ± }- 1/2 = v/28/3 A" 1 ~ 0.47 fm . (9.58) 

As expected, it is smaller than the charge radius which we attribute to the fact 
that the charge distribution measured by the charge radius does not coincide with 
the distribution of baryon density. The core radius is sometimes related to 



the decay constant /„- via the dimensionless quantity [227 



C = UR n . (9.59) 
In constituent quark models one typically gets 

/. = ^ , (9-60) 



146 



CHAPTER 9. THE PION WAVE FUNCTION IN THE NJL MODEL 



and if the quarks are confined by a harmonic oscillator potential, their mass 
is m = 3/2R n , so that C ~ 0.4. This implies the fairly large value R n ~ 
0.8 fm. Using standard many-body techniques, Bernard et al. have calculated 
this quantity in a model treating the pion as a collective excitation of the QCD 
vacuum, and find C ~ 0.2 ||228[| . This result is close to what we get from fl9.51| ) 
and (|938l) , 




C = — — \ — ~ 0.24 . (9.61) 
256 V 3 v ; 

Thus, though we work within a constituent picture, we do get a reasonable value. 
We believe that this is due to the intrinsic consistency of the light-cone framework 
with the requirements of relativity. 
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9.4 The Pion Distribution Amplitude 



We follow Brodsky and Lepage [82, 216] and define the pion distribution ampli- 



tude by integrating the wave function over k±. Due to the cutoff dependence of 
the wave function, the integration extends effectively only up to the cutoff scale 
A(x), 

0bl(x) = J k J_, A, A') . (9.62) 

As usual, we evaluate this to leading order in e so that the result becomes inde- 
pendent of the constituent mass, 

0bl(x) = [x(l - x)f 2 = [x(l - x)f 2 . (9.63) 

47T 97T 

The normalization of <fi is obtained by integrating this over x, which simply 



yields the 'wave function at the origin', ip , defined in (|9.48| ). We have used 
the subscript 'BL' (for Brodsky-Lepage), because in what follows we will use a 
distribution amplitude which is normalized to one, that is, 

^)-^0bl = ^[x(1-x)] 3 / 2 . (9.64) 

It has been shown by Brodsky and Lepage, that the distribution amplitude obeys 
an evolution equation in momentum scale Q 2 which can be solved exactly in the 
limit Q 2 — > oo. The result is called the asymptotic distribution amplitude, 

as (x) = 6x{l - x) . (9.65) 

The (non-asymptotic) pion distribution amplitude has been a rather controversial 
object. For a while people have tended to believe in a 'double-humped' shape of 
the amplitude (due to a factor (1 — 2x) 2 ), which was originally found by Chernyak 



and Zhitnitsky via QCD sum rules ||229|| . Recent CLEO data, however, ||230|| seem 



to support an amplitude that is not too different from the asymptotic one ||231|| . 



It turns out, that our result ( 9.64j) is qualitatively consistent with these findings 



232, 233, 234 . Belyaev and Johnson have recently reported two constraints 



which should be satisfied by the distribution amplitude [p3 

0( X = O.3) = 1 ± 0.2 , 

0( X = O.5) = 1.25 ±0.25. (9.66) 

In Figure |T2] we have displayed our distribution amplitude in comparison with the 
asymptotic one. The vertical lines represent the constraints (|9.66j ). Obviously, 
these are not satisfied, as we have 

(j)(x = 0.3) = 1.31 , (9.67) 
0( x = 0.5) = 1.70 . (9.68) 
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As we already overshoot the asymptotic amplitude (for 0.3 < x < 0.7) this does 
not come as a surprise. In order to improve on that we presumably have to 
increase our accuracy by taking into account finite pion and quark masses. Still, 
we get an amplitude which is at least qualitatively correct. 



Figure 9.2: The pion distribution amplitude. The full curve represents the asymp- 
totic amplitude, the broken one the amplitude calculated within the NJL model. 
The vertical lines mark the constraints (|9.6'6|) of Belyaev and Johnson. 



(x) 1 




It should also be mentioned that it is not obvious what the actual momentum 
scale is to which our result corresponds. Our transverse-momentum cutoff is x- 
dependent, A 2 (x) ~ A 2 x(l — x), so, if we choose an average x of (x) ~ 1/2, a 
natural scale seems to beQ 



1/2 



(x)(l-(x)) A ~ A/2 ~ 600 MeV . 



(9.69) 



This is somewhat lower than fx ~ 1 GeV, which has been assumed by Belyaev 
and Johnson in their analysis of the distribution amplitude in terms of light-cone 
quark models ||233|| . 



We thank W. Schweiger for discussions on this point. 



Conclusions 



In this report we have discussed an alternative approach to relativistic (quantum) 
physics based on light-cone dynamics. It makes use of the fact that for relativistic 
systems the choice of the time parameter is not unique. Our particular choice 
of relativistic dynamics singles out Dirac's front form so that null-planes tangent 
to the light-cone become hypersurfaces of equal time. This apparently trivial 
change of coordinates has far-reaching consequences: 

• The number of kinematical (i.e. interaction independent) Poincare gener- 
ators becomes maximal; there are seven of them instead of the usual six, 
among them the boosts. 

• Lorentz boosts in z-direction become diagonal; the light-cone time and 
space coordinates, x + and x~, respectively, do not get mixed but rather get 
rescaled. 

• As a consequence, for many-particle systems one can introduce frame- 
independent relative coordinates, the longitudinal momentum fractions, Xj, 
and the relative transverse momenta, k_u. 

• Because of a two-dimensional Galilei invariance, relative and center-of- 
mass motion separate. As a result, many formulae are reminiscent of non- 
relativistic physics and thus very intuitive. 

• This is particularly true for light-cone wave functions which, because of 
the last two properties, are boost invariant and do not depend on the total 
momentum of the bound state. They are therefore ideal tools to study 
relativistic particle systems. 

• The last statement even holds for relativistic quantum field theory where 
one combines the unique properties of light-cone quantization with a Fock 
space picture. The central feature making this a reasonable idea is the 
triviality of the light-cone vacuum which accordingly is an eigenstate of the 
fully interacting Hamiltonian. This implies that the Fock operators create 
the physical particles from the ground state. 
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As expected, however, the principle of 'conservation of difficulties' is at work 
so that there are problems to overcome. We have focused on the issue of vacuum 
properties. In instant-form quantum field theory, in particular in QCD, many 
non-perturbative phenomena are attributed to the non-triviality of the vacuum 
which shows up via the appearance of condensates. These are non- vanishing vac- 
uum expectation values which support the idea that the instant-form vacuum is a 
complicated many-body state (like, e.g. the BCS ground state). Of further phys- 
ical relevance is the fact that many of these condensates signal the spontaneous 
breakdown of a symmetry. The conceptual problem which arises at this point is 
how to reconcile the existence of condensates with the triviality of the light-cone 
vacuum. 

We have shown that for scalar field theory in 1+1 dimensions, the field mode 
with longitudinal momentum k + is the carrier of non-trivial vacuum properties. 
Its c-number part coincides with the vacuum expectation value of the field. The 
reflection symmetry, mapping the field onto its negative, is thus seen to be spon- 
taneously broken. 

For fermionic fields with a bilinear condensate reflecting the spontaneous 
breakdown of chiral symmetry, we had to proceed differently. We calculated 
the bound state spectrum of model field theories and reconstructed ground state 
properties like the condensates from the particle spectrum. The paradigm iden- 
tity was the Gell-Mann-Oakes-Renner relation which relates the masses of the 
bound states and the fermions to the condensates. Alternatively, one can try to 
derive gap equations like in the NJL model, which determine the mass of the 
dressed, effective fermions in terms of the condensate. 

It is crucial for the latter approach to solve the associated bound-state prob- 
lem, or, in other words, the light-cone Schrodinger equation. For a relativistic 
quantum field theory, this, in principle, amounts to solving an infinite system of 
coupled integral equations for the amplitudes to find an ever increasing number 
of constituents in the bound state. Experience, however, shows that the light- 
cone amplitudes to find more than the valence quanta in the bound state tend 
to become rather small. We have confirmed this fact for the theories we have 
studied. Note that the same is not true within ordinary, that is, instant-form 
quantization. The light-cone approach, however, provides some confidence that 
also in real QCD a constituent picture may emerge, so that hadrons are consis- 
tently described as light-cone bound states of a minimal number of two or three 
constituents. 

A first step in this direction was performed using an effective field theory, the 
NJL model. We have seen that, though we made a number of approximations, 
in particular by enforcing a constituent picture, a number of pionic observables 
are predicted with reasonable accuracy. 

At this point one has to ask whether the ideas just outlined can really be used 
in QCD. This is a difficult question. The reason why light-cone QCD has not 
been presented in this report is that, in our opinion, there are still technical prob- 
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lems to overcome. First of all, the quantization procedure itself is not entirely 
clear beyond perturbation theory There are the usual gauge redundancies to be 
taken into account together with the peculiar first-order nature of the light-cone 
Lagrangian. We have reported on this elsewhere |[235|| . One generally uses the 
light-cone gauge A + = 0, a variant of an axial gauge which has its own well-known 
problems, among them residual gauge invariance and infrared divergences. How- 
ever, it is only in this gauge that the canonical quantization is straightforward. 
The latter has been performed quite some time ago ||236| , |237| , |217|| and success- 
fully been used as the basis of light-cone perturbation theory. Some impressive 
achievements of the latter method can be found in the recent review |[238|| which 
focuses on the light-cone formulation of gauge theories. 

Still, even if one wants to utilize perturbation theory for bound-state cal- 
culations (in case this is justified) one is confronted with difficulties. One big 
problem is to find a systematic renormalization program in order to get rid of the 
ubiquitous infinities. As one does neither have explicit covariance nor rotational 
invariance there are very few guiding principles and restrictions. Accordingly, 
there is a rather large variety of divergences which can even become nonlocal, 
power counting gets more complicated, and the number of counter-terms tends 
to become large ||137| . As a result, to our knowledge, the light-cone renormaliza- 



tion program has not been pushed beyond one loop so far. 

We therefore consider it worthwhile to first pursue a less ambitious program, 
the first steps of which have been undertaken in the last chapter. One should 
make use of effective field theories, that are obtained from QCD in a well-defined 
239], [240|| , and in general are not renormalizable. The cutoff is a pa- 



manner 



rameter of the theory which sets the energy scale below which the theory makes 
sense. Thus, one only has to deal with the issue of regularization to render the 
appearing integrals finite. Still, for the case of the NJL model, we have seen 
that already this can be a nontrivial problem. An alternative candidate suited 
for quark bound-state studies is the effective quark theory derived within the 
QCD instanton model [IE|. It is also of the NJL type, however, with a more 
complicated (and presumably more realistic) interaction. 

Further improvements can be achieved by using more elaborate kernels in the 
covariant Bethe-Salpeter equation. These can be found in the literature on the 
Schwinger-Dyson approach to QCD ||206| , [20711 and basically amount to choosing 
a model for the gluon propagator. The NJL model corresponds to an infrared- 
finite propagator which is constant in momentum space. Therefore, one does not 
have confinement. The latter can be implemented by using a gluon propagator 
which diverges for small momenta. Other constraints follow from multiplicative 
renormalizability. Clearly, the procedure of deriving the light-cone bound-state 
equation will no longer be as simple as in the last chapter but should still be 
feasible. Work in this direction is underway. 
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Appendix A 
Light-Cone Spinors 



In this appendix we collect some basic formulae representing our conventions for 
light-cone spinors in d = 3+1, together with some frequently occurring identities 
involving these. Our conventions for Dirac matrices are those of Bjorken and 
Drell 0. 

Let ip{x) be a solution of the free Dirac equation 

(0 - m)ip(x) = . (A.l) 
The Fock expansion of if) reads 

#s,0) = ^ fjrj^ [b(k,XHk,X)e- ik - x + d\k,X)v(k,X)e ik - x 

(A.2) 

where 

k = (k + , kj_) , x = (aT, xj_) , k ■ x = \k + x~ — kj_ ■ xj_ . (A. 3) 

The projections ip± are obtained through the action of the projectors A±, ip± = 
A±tp, with (in a two-by-two block notation) 

1 T ± _ 1 ( 1 ±<r 3 



» ]■ (kA) 

As A + and A_ are projectors they obey 

A| = A±, Ai = A ± , (A.5) 

and 

A + + A_ = 1 , A + A_ = A_A+ = . (A.6) 
The following identities are frequently used (i = 1,2), 

(A.7) 
(A.8) 
(A.9) 



A±/3 


= 


A ± f 


= f A± 


A±a i 


= a { A T 
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APPENDIX A. LIGHT-CONE SPINORS 



The spinors u and v in ( |A.2| ) satisfy the equations 



— m) w(fc, A) 
+ m) v(k, A) 








and are explicitly given by [[24 1| , 216 , |86 

1 



u(k, A) 
u(fc, A) 



1 



— [3m + CKj^) X_ A . 



The X\ are eigenspinors of the projector A + with eigenvalue one, 



A + X\ — X\ 



X X X X , — 5\\> 



or explicitly, 



X x =X 



1 

7! 



i 



1 





4 



1 

7! 



o 
i 



-1 



If we introduce the two-spinors x A , the usual eigenspinors of cr 3 , 



1 





we can write for X\, 



X— i — X I 



Xa 





1 



Using the identities 



A, 







A, 







we can rewrite the elementary spinors ( |A.12j , |A.13| ) according to 



u(k, A) 
v(k,\) 




— [mA_ + (k + + a i k i )K + ] 







- [mA„ + (k + + c/F)A+] 







(A.10) 
(A.ll) 



(A.12) 
(A.13) 



(A. 14) 



(A.15) 



(A.16) 



(A.17) 



(A.18) 



(A.19) 
(A.20) 
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The projections of the spinors u and v are 

u+(k,X) = Vk+Xx = u + (k + , A) , (A.21) 
v+(k,X) = Vk+X^ x = v + (k + ,X) , (A.22) 

w_(ife,A) = -= (m/3 + Oik) X x , (A.23) 
V k + 

u_(fe,A) = -7= (-m/3 + X_ A . (A.24) 

v fc + 

One has the completeness relation in the '+ '-sector, 

u +(k, X)ul(k, X) = J2 A )4( fe > A ) = fc+A + • ( A - 25 ) 

A A 

The following provides a (very incomplete) list of matrix elements (for further 
matrix elements, see ||216|| ), 



75 = i7 Q 7V7 3 = ( ? J ) , (A.30) 



?4(A; + ,A)« + (A; + ,A0 = w 1 ( .(A; + ,A)t; + (A; + ,A / ) = A;+5 AV , (A.26) 

w^(fe,A)^« T (fe,A') = rndxy =F A£;_a<5a,-a' , (A.27) 

*4(fe,A)/?%(M') = -m^±AM-A,A', (A.28) 

where we have defined 

k x = k 1 + i\k 2 , so that k x k_ x = k\ . (A.29) 

Introducing 

■ 1 2 "? / 

<Vr = 7.0/ 0/ 0/ 0/ = I 

x 1 

one finds that X A is also an eigenspinor of 75, 

l5 X x = XX X . (A.31) 

Furthermore, one has 

k'-y'Xx = XkxX-x ■ (A.32) 
The last two identities lead to the matrix elements 

u ] + {k\X) 1 ,v + {p\X') = V^Y^va', (A.33) 

mU^AHYTsM^A') = ^/t^k_ x 5xx>, (A.34) 

which have been used in Chapter 9 to determine the spin structure of the pion 
wave function. 
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